ESD  TDR  64-394 


* 


>* 

a, 

8 

8 

i — i 

u_ 


■  £ 
ULJ 


ESD-TDR-64-394 


EST1  PRC 

ESO  RCCORO  COPY  □  DOC  T AH  □  FROJ  OFFICEH 

Cr.r.T  RETURN.T'0  O  ACCESSION  M*5TE»  FILE 

SCIENTIFIC  4  TECHNICAL  »R!MIAT?W  mflll 

(tsri),  suiioifia  un  " 


□ 


COPY  NR 


u 


OF 


+ 


CO*£S 


DATE, 


CV9 


ESPOD  OPERATING  INSTRUCTIONS""""^^^ 

CY  NIY _ L -  OF - J- 

AND  CARD  FORMATS 

TECHNICAL  DOCUMENTARY  REPORT  NO.  ESD-TDR-64-394 


JUNE  24,  1964 


496L  Systems  Program  Office 

ELECTRONIC  SYSTEMS  DIVISION 
AIR  FORCE  SYSTEMS  COMMAND 
UNITED  STATES  AIR  FORCE 

L,  G.  Hanscom  Field,  Bedford,  Massachusetts 


PROJECT  ES-3-496L-3627 

Prepared  under  Contract  AF  19(628)-594 


STL  NO.  8497-6066- RU -000 


PfoDloDkOW- 


ESD-TDR-64-394 


ESPOD  OPERATING  INSTRUCTIONS 
AND  CARD  FORMATS 

TECHNICAL  DOCUMENTARY  REPORT  NO.  ESD-TDR-64-394 

June  24,  1964 


496L  Systems  Program  Office 
Electronic  Systems  Division 
Air  Force  Systems  Command 

United  States  Air  Force 
L.  G.  Hanscom  Field,  Bedford,  Massachusetts 


Project  ES-3-496L-3627 
Prepared  Under  Contract  AF  19(628)-594 


STL  No.  8497- 6066- RU- 000 


FOREWORD 


This  document  is  one  of  three  reports  which  describe  ESPOD,  a  general 
satellite  orbit  determination  program  prepared  for  the  Air  Force  Elec¬ 
tronic  Systems  Division  for  use  in  the  SPACETR ACK/ SPAD ATS  Center 
at  Ent  Air  Force  Base,  Colorado  Springs,  Colorado. 

This  report  is 


ESD-TDR-64-394 


ESPOD  Operating  Instructions 
and  Card  Formats 


(STL  No.  849 7- 6066 -RUO 00) 


The  companion  reports  are 


ESD-TDR-64-393 


ESPOD  Functional  Description 
(STL  No.  8497-606 7- RU000) 


ESD-TDR-64-395 


ESPOD  Mathematical  and 
Subroutine  Description 
(STL  No.  849 7- 6065- RU 000) 


The  ESPOD  Program  was  prepared  by  TRW  Space  Technology  Labora¬ 
tories  under  Air  Force  Contract  Number  AF  19(628)-594. 
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1.  INTRODUCTION 


This  report  provides  operating  instructions  and  other  aids  pertinent 
to  the  ESPOD  program.  The  sections  are  arranged  in  the  order  that  they 
will  be  used  to  set  up  and  run  the  ESPOD  program. 

a)  Input  Control  Card  Descriptions  and  Deck  Setup 

b)  Operational  Instructions 

c)  Output  Descriptions 

d)  Program  Constants 

The  operator  or  analyst,  making  frequent  use  of  this  report,  will 
find  the  following  features  useful. 

a)  Card  formats  are  exact  full  scale  size  with  unused  or 
blank  areas  cross  hatched. 

b)  A  separate  listing  of  titles  for  preliminary  data  cards. 

c)  Tables  providing  diagnostic  assistance  in  card  usage 
and  program  options. 

For  an  abridged  technical  description  of  the  ESPOD  Program, 
ESD-TDR-64-393,  "ESPOD  Functional  Description,"  may  be  consulted. 

The  mathematical  and  computer  processes  for  the  ESPOD  Program 
are  described  in  ESD-TDR-64-395,  "ESPOD  Mathematical  and  Subroutine 
Description. " 
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2.  INPUT  INSTRUCTIONS  AND  CARD  FORMAT 


This  section  is  concerned  with  the  detailed  description  of  the  ESPOD 
input  quantities.  The  first  part  of  the  section  describes  the  card  format, 
content,  and  groupings  of  the  B2  systems  cards,  the  ESPOD  program  con¬ 
trol  cards,  and  the  program  option  cards.  The  standard  SPADATS  system 
observations  cards  and  sensor  cards  which  ESPOD  uses  are  also  described. 

The  second  part  of  this  section  is  concerned  with  the  input  sequence. 
Sample  ESPOD  load  sheets  with  the  corresponding  input  decks  are  de¬ 
scribed.  The  various  categories  of  the  input  cards  are  functionally 
grouped  and  tabulated  as  to  input  mode  to  aid  the  analyst  in  preparing 
an  ESPOD  run.  There  is  also  a  complete  tabulation  of  the  program 
response  to  omitted  or  included  cards,  or  both. 

2.  1  INPUT  CONTROL  CARD  DESCRIPTION 

Cards  for  input  to  ESPOD  fall  into  six  classes  and  are  listed  below 
in  their  normal  input  sequence: 

1)  B2  System  Job  Start  Cards 

2)  Job  Description  Card,  JDC 

3)  Preliminary  Data 

4)  Observation  Cards 

5)  Sensor  Cards 

6)  B2  System  Job  Finish  Cards 

All  classes  are  not  necessarily  needed  for  a  particular  run.  This 
will  be  discussed  at  length  in  Section  2.  2. 

2.  1.  1  B2  System  Inputs 

There  are  two  types  of  input  cards  which  must  be  used  with  ESPOD 
inputs  when  running  in  the  scheduled  mode.  The  Job  Start  and  Job  Finish 
input  cards  are  described  on  the  following  pages. 
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70  SCHTP 


70  SCHTP 


Field 

Column 

Content 

Description 

1 

1-8 

70  A  SCHTP 

2 

9 

11-8-2 

11-8-2  multipunch 

3 

10-80 

Not  used 

Figure  2-1.  70  SCHTP  Card  Format 


2-3 


SPADATS  JOB 


SPADATS  JOB 


000000000000000 

111111111111111 

2222222222222222 

3333333333333353 

4444444444444444 

5555555555555555 

6666666666666666 

7777777777777777 

8888888888888388 

9999999999999999 

2j  2i  27  ?»  21  M  J1  r  13  W  »  *  37  U  39  4C 


Field 

Column 

Content 

De  scr  iption 

1 

1-16 

Not  used 

Z 

17-19 

JOB 

3 

20-24 

Not  used 

4 

25-40 

NNNN 

Job  number 

5 

41-79 

Not  used 

6 

80 

J 

Figure  2-2  SPADATS  Job  Card  Format 


SPADATS  REMARKS 


SPADATS  REMARKS 


Field 

Column 

Content 

Desc  r  iption 

1 

1-16 

Not  used 

2 

17-19 

REM 

3 

20-24 

Not  used 

4 

25-40 

xxxx 

Arbitrary  remark 

5 

41-80 

Not  used 

Figure  2-3.  SPADATS  Remarks  Card  Format 
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RUN 


RUN 


Field 

Column 

Content 

Description 

1 

1  - 1 6 

Not  used 

2 

17-19 

RUN 

3 

20-24 

Not  used 

4 

25-34 

ESP0D, 

DATA 

5 

35-80 

Not  used 

Figure  2-4.  RUN  Card  Format 
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ENDDATA 


ENDDATA 


Field 

Column 

Content 

Description 

1 

1-16 

Not  used 

2 

17-23 

ENDDATA 

3 

24-80 

Not  used 

Figure  2-5.  ENDDATA  Card  Format 
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END0FJ0B 


END0FJ(Z>B 


Field 

Column 

Content 

De  scr  lption 

1 

1-8 

END0FJ0B 

2 

9 

11-8-2 

11-8-2  Multipunch 

3 

10-80 

Not  used 

Figure  2-6. 


End  of  Schedule  Mode  Card  Format 
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ENDSCHED 


ENDSCHED 


Field 


Columns 


Content 


1 


1-8 


ENDSCHED 


2 


9 


11-8-2  multipunch 


3 


10-80 


Not  used 


Figure  2-7.  Schedule  Tape  END  Card 


2-15 


2.  1.  2  JDC  —  Job  Description  Card 


The  JDC  serves  as  a  control  card  for  the  flow  of  information  through 
ESPOD.  This  card  is  always  the  first  card  of  an  input  deck.  It  defines 
the  program  sections  to  be  used,  and  selects  certain  options.  It  identifies 
the  vehicle  number  and  name  and  permits  a  short  arbitrary  remark. 
Column  description  of  the  JDC  card  follows.  A  pictorial  description  of 
the  JDC  card  appears  at  the  top  of  the  load  sheet  (see  Section  2.  2.  3  for 
illustration) . 

2.1.3  Preliminary  Data 

The  preliminary  data  input  cards  control  the  detailed  options  avail¬ 
able  to  ESP0D  and  provide  numerical  information  specific  to  the  case. 
Preliminary  data  must  be  provided  for  all  cold- start  runs  except  for 
those  which  only  sort  and  list  the  observations.  An  ENDPR  card  must 
follow  the  last  of  the  preliminary  data  cards.  The  full  set  of  fields  re¬ 
quired  to  set  up  the  preliminary  data  cards  is  illustrated  in  Figure  2-9. 

It  should  be  noted  that  they  have  a  fixed  format.  An  Identifier  Field 
(I-Field)  (3,  5,  7  or  9,  Figure  2-9)  is  five  columns  wide.  A  Variable 
Field  (V-Field)  (4,  6,  8  or  10,  Figure  2-9)  is  14  columns  wide.  All 
Identifiers  are  left  adjusted  in  the  I-Field.  All  numerical  values  appear¬ 
ing  alone  in  a  V-Field  must  have  a  decimal  point  indicated  somewhere 
and  if  the  numerical  value  is  expressed  in  floating  decimal  format,  the 
exponent  must  be  right  adjusted  in  the  V-Field.  A  value  in  a  V-Field 
which  is  not  preceded  by  a  name  in  the  associated  I-Field  will  be  consid¬ 
ered  as  a  value  for  an  element  in  the  array  identified  in  the  first  I-Field- 1 
of  the  card.  Subsequent  values  after  the  first  are  put  into  consecutive 
locations  in  the  array.  If  a  V-Field  is  left  blank,  the  field  is  ignored. 

The  following  is  an  example  of  the  same  value  entered  into  its  V-Field 
in  two  different  forms,  fixed  decimal  and  floating  decimal  formats. 


Column  number 

10 

1  1 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

Fixed  decimal 

- 

1 

4 

5 

3 

2 

6 

7 

Floating  decimal 

- 

1 

4 

5 

3 

2 

6 

7 

+ 

3 

- 

0 

♦ 

1 

4 

5 

3 

2 

6 

7 

+ 

0 

3 

Floating  decimal 

- 

0 

1 

4 

5 

3 

2 

6 

7 

+ 

0 

2 

Floating  decimal 

- 

0 

1 

4 

5 

3 

2 

6 

7 

+ 

0 

0 
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JDC 


JDC 


Column 

Content 

Description 

1-3 

JDC 

Identifies  JDC  card 

4-7 

XXXX 

Vehicle  number- right  adjusted 

8-17 

XXXX 

Vehicle  name 

18-29 

XXXX 

User’s  header 

30 

0 

Cold  start 

1 

Conditioned  start 

2 

Conditional  start 

31 

0 

Sensor  and  observation  data 
not  to  be  processed 

1 

Sensor  and  observation  data 
to  be  processed 

32 

0 

Do  not  print  sensor  data 

1 

Print  sensor  data 

33 

0 

Do  not  print  observations 

1 

Print  observations 

34 

0 

Sensors  on  SEAI  tape 

1 

Sensors  not  on  SEAI  tape 

35 

0 

Observations  on  SRADU  tape 

1 

Observations  not  on  SRADU  tape 

Figure  2-8.  Job  Description  Card  Format 
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JDC 


JDC 


Column 

36 

37-40 

41 


42 


43 

44 

45 


46 

47-50 

51 

52 


53 


54 

55 


56 


57-80 


Content 

0 

1 

Not  used 

0 

1 

0 

1 

0 

1 

0 

1 

0 

1 

2 

0 

1 

Not  used 
0 
1 

0 

1 

0 

1 

Not  used 

0 

1 

0 

1 

Not  used 


Description 

Do  not  print  program  constants 
Print  program  constants 

ESP0DDC  not  needed 

ESP0DDC  needed 

A  priori  S  matrix  not  input 
this  run 

A  priori  S  matrix  input  this  run 
T  -  1 

Do  not  punch  (A  A)  this  run 
T  -  1 

Punch  (A  A)  this  run 
Do  not  punch  A^A  this  run 
Punch  A  1  A  this  run 
UVW  residuals  print 
STW  residuals  print 
LLH  residuals  print 
Use  proved  elements 
Use  new  elements 

ESP0DEPH  not  needed  this  run 

ESP0DEPI1  needed  this  run 

x,  v,  z,  T^,  not  present 

this  run 

x,  y,  z,  T\  ,  Trn  present  this  run 
L)  r  i) 

Do  not  generate  special  binary 
empheris  tape  on  this  run 

Generate  special  binary 
ephemeris  tape  on  this  run 

Do  not  perform  matrix  update 
this  run 

Perform  matrix  update  this  run 

Do  not  punch  inverse  of  updated 
covariance  matrix 

Punch  inverse  of  covariance 
matrix 


Figure  2-8.  Job  Description  Card  Format  (Continued) 


2-19 


PRELIM  DATA 


PRELIM  DATA 


Field 

No. 

Field 

Name 

Column 

Content 

Description 

1 

1-2 

NN 

Sequence  number  identifying 
successive  cards  starting  with 
first  and  having  the  same  identi¬ 
fier  in  Field  3  or  special  identifier 
99 — identifies  99  card 

DAC  Cards 

SIGMA  Cards 


2 

3-4 

NN 

Iteration  number;  identifie s  punched 
output  from  successive  iterations 

3 

I-Field  1 

5-9 

xxxxx 

Principal  variable  Identifier 

4 

V -Field  1 

10-23 

xxxxx 

Value 

of  first  variable 

5 

I-Field  2 

24-28 

xxxxx 

Name 

of  second  variable 

6 

V  -  Field  2 

29-42 

xxxxx 

Value 

of  second  variable 

7 

I-Field  3 

43-47 

xxxxx 

Name 

of  third  variable 

8 

V -Field  3 

48-61 

xxxxx 

V  alue 

of  third  variable 

9 

I-Field  4 

62-66 

xxxxx 

Name 

of  fourth  variable 

10 

V -Field  4 

67-80 

xxxxx 

Value 

of  fourth  variable 

F  igure 

2-9.  Field 

F  ormat 

of  Preliminary  Data  Card 
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To  readily  locate 
of  the  card  formats  by 
below. 

Card  Title 

ICTYP 

1C0ND  No.  1 

IC0ND  No.  2 

ITIME  No.  1 

1TIME  No.  2 

SPADATS-7  No. 

SPADATS- 7  No. 

SPADATS-7  No. 

SPADATS-7  No. 

SPADATS-7  No. 

SPADATS-7  No.  i 

SPADATS-7  No. 

DNREV 

CAT  1 

CAT2 

BNDS 

BISES 

NITER 

TSTEP 

DELET 

STYPE 

T  MAX 

SIGMA 

DRAG 

APF10 

Z0NAL 

SECTR 

TESSR 

Z0NAL-M0D2 

99-M0D2 


any  particular  preliminary  data  input  card  a  list 
title  in  order  of  apprarnncr  in  this  section  is  given 


Function  in  ESPOD 

Page 

Initial  conditions 

2-25 

Initial  conditions 

2-27 

Initial  conditions 

2-29 

Initial  conditions 

2-  11 

Initial  conditions 

2-3  3 

Initial  conditions 

2-35 

Initial  conditions 

2-  17 

Initial  conditions 

2-39 

Initial  conditions 

2-41 

Initial  conditions 

2-4  1 

Initial  conditions 

2-45 

Initial  conditions 

2-47 

Initial  conditions 

2-49 

Solution  variables 

2-51 

Solution  variables 

2-55 

Solution  variables 

2-57 

Solution  variables 

2-61 

Convergence  control 

2-63 

Integration  control 

2-65 

Data  editing 

2-67 

Data  editing 

2-69 

Data  editing 

2-71 

Data  editing 

2-73 

Drag 

2-75 

Drag 

2-77 

Earth  potential 

2-83 

Earth  potential 

2-85 

Earth  potential 

2-86 

Earth  potential 

2-89 

Earth  potential 

2-91 

2-22 


Card  Title 

Function  in  ESPOD 

Page 

SECTR-M0D3 

Earth  potential 

2-93 

SECTR-M0D4 

Earth  potential 

2-95 

TESSR-M0D4 

Earth  potential 

2-97 

RADPR 

Radiation  pressure 

2-99 

SMAT 

Matrix  inputs 

2-101 

UPMAT 

Matrix  inputs 

2-105 

DELTT 

ESP0DEPH  options 

2-111 

DAC  No.  1 

ESP0DEPH  options 

2-113 

DAC  No.  2 

ESP0DEPH  options 

2-115 

DAC  No.  3 

ESPCJDEPH  options 

2-117 

PRDCT 

ESP0DEPH  options 

2-119 

99 

Constants 

2-121 

ENDPR 

End  preliminary  data 

2-123 

The  following  sections 

give  the  details  of  the  various 

inputs  on  the 

preliminary  data  cards  and  how  to  specify  the  available  options. 


2 .  1.3.  1  Specification  of  Initial  Conditions 


The  initial  conditions  can  be  specified  in  any  of  four  different  types 
of  input  coordinates  defined  in  sections  which  follow.  The  input  cards  can 
appear  in  three  different  forms. 

•  The  ICTYP  card  specifying  the  coordinate  system;  the 
IC0ND  cards  providing  the  initial  conditions;  and  the 
ITIME  card  giving  the  epoch  for  the  input  solution  vector 
completely  define  the  initial  conditions. 

•  The  ICTYP  card  specifying  SPADATS  7- card  elements; 
and  an  optional  DNREV  card  (to  update  to  new  epoch)  define 
a  set  of  initial  conditions. 


•  A  third  form  is  obtaining  initial  conditions  from  the  SEAI 
tape  and  an  optional  DNREV  card. 

2.  1.  3.  1.  1  ICTYP  Card  Description.  The  variable  field  associated  with 
this  name,  derived  from  Initial  Condition  Type,  is  an  integer  between 
1  and  4  which  indicates  the  type  of  initial  conditions,  input  by  the  IC0ND 
set  or  by  the  7- card  elements.  The  coordinate  systems  and  type  are 
listed  below: 
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1.  ECI  Polar  Spherical,  Right  Ascension  Reference  (ADBARV) 


Q 

Right  ascension 

(degrees) 

6 

Declination 

(degrees) 

p 

Flight  path  angle 

(degrees) 

A 

A  zimuth 

(degrees ) 

R 

Radius 

(km) 

V 

Velocity 

(km/ sec) 

ECI  Polar  Spherical,  Longitude  Referen 

V 

Longitude 

(degrees) 

6 

Declination 

(degrees) 

P 

Flight  path  angle 

(degrees) 

A 

Azimuth 

(degrees) 

R 

Radius 

(km) 

V 

Velocity 

(km/  sec) 

ECI 

Carte  sian 

x  1 

I 

(km) 

>  Position  coordinates 

(km) 

.  J 

i 

(km) 

*  1 

i 

(km/  sec) 

y 

\  Velocity  coordinates 

(km/  sec) 

Z  _> 

l 

(km/ sec) 

SPADATS -7-Card  Mean  Elements 

a 

Semimajor  axis 

(km) 

e 

Eccentricity 

(degrees) 

i 

Inclination 

(degrees) 

P 

n 

Nodal  period 

(days / rev. ) 
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ICTYP 


ICTYP 


Field 

Column 

Conte  nt 

1 

1-2 

01 

2 

3-4 

Blank 

3 

5-9 

ICTYP 

4 

10-11 

1  Z..  3.. 

or  4 . 

4 

12-23 

Not  used 

4-10 

24-80 

Not  used 

Description 
Sequence  number 
Iteration  number 
Principal  variable  identifier 


F igur  e  Z  -  I  0. 


ICTYP  Card  Format 


Z-Z5 


ft  Right  ascension  of 

ascending  node  (degrees) 

°°0  Argument  of  perigee  (degrees) 

t  Time  of  epoch  (days) 

2.  1.3.  1.2  ICCftND.  This  name,  derived  from  Initial  Conditions,  identi¬ 
fies  the  initial  position  and  velocity  for  the  following  coordinates: 

1.  (ICTYP  No.  1)  ECI  polar  spherical  —  a,  6,  (3,  A,  R,  v 

2.  (ICTYP  No.  2)  ECI  polar  spherical  —  6,  (3,  A,  R,  v 

3.  (ICTYP  No.  3)  ECI  Cartesian  —  x,  y,  z,  x,  y,  z 

The  IC0ND  cards  are  only  for  these  first  three  types  of  initial  conditions. 
The  SPA  DATS  7- card  element  set  is  discussed  in  Section  2.  1.3.  1.4. 

Since  six  variables  specify  the  initial  conditions,  two  cards  constitute  an 
IC0ND  set.  The  IC0ND  card  format  is  delineated  by  field  and  content 
(see  Figures  2-11  and  2-12).  The  field  for  mat  cor  responds  to  the  standard 
preliminary  data  card  and  the  Load  Sheet  description  (Section  2.  2). 
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IC0ND  No.  1 


IC0ND  No.  1 


ojo  0000000000  0  OOj 

It  11  1J  13  1«  IS  II  17  1|  II  n  21  n  2 
111111111111111 

2222222222222^ 

3|3  333333333333  3c 

4444444444444  4^ 

5|5  5  5  5  5  5  5  5  5  5  5  5  5  5§ 

66666666666666 6§ 

7777777777777  7b 

B  8  8  8  8  8  8  8  8  8  8  8  8  8  8f 

999999999999999 

id  11  1?  13  14  15  II  17  II  II  20  J1  22  23 


00000000000000 

21  10  31  13  13  M  JS  II  37  »  M  40  41  *T 

11111111111111 

2  2  2  2  2  2  2  2  2  2  2  2  2  2t 
33333333333333| 
44444444444444; 
55555555555555J 
66666666668666; 
7777777777777?! 
88888888886888] 
9999999999999  9p 

21  »  31  32  U  34  35  »  37  31  3»  41 4142 


00000000000000 

41  41  so  si  S?  S3  54  55  U  J7  M  St  M  II 

11111111111111 
12  2222222222222 
]  3  3333333333333 
44444444444444 

[55555555555555 
(66666666666666 
7  7  7  7  7  7  7  7  7  7  7  7  7  7; 


9999999999999$ 

[  41  41  SI  SI  S?  S3  54  55  %  57  SI  bl  60  41 


10 

0  0  0  0  0  0  0  0  0  0  0  0  I  0 

lunnunnnnnnnnm 

11111111111111 

22222222222222 

33333333333333 

44444444444444 

55555555555555 

66666666666666 

77777777777777 

88888888888888 

99999999999999 

17  M  H  70  71  n  >3  74  75  71  V  l{  75  H 


Content  of  Field 


Field 

Columns 

ICTYP 

1 

ICTYP 

2 

ICTYP 

3 

1 

1-2 

01 

01 

01 

2 

3-4 

Blank 

Blank 

Blank 

3 

5-9 

ICQ)ND 

IC0ND 

ICCjND 

4 

10-23 

a 

X 

X 

5 

24-28 

Blank 

Blank 

Blank 

6 

29-42 

6 

6 

y 

7 

43-47 

Blank 

Blank 

Blank 

8 

48-61 

p 

P 

z 

9 

62-66 

Blank 

Blank 

Blank 

10 

67-80 

A 

A 

X 

Figure  2-11.  IC0ND  Card  Format  No.  1 
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ICOND  No.  2 


ICOND  No.  2 


Content  of  Field 


Field 

Columns 

ICTYP 

1 

ICTYP 

2 

ICTYP 

3 

1 

1-2 

02 

02 

02 

2 

3-4 

Blank 

Blank 

Blank 

3 

5-9 

IC0ND 

IC$ND 

IC0ND 

4 

10-23 

R 

R 

y 

5 

24-28 

Blank 

Blank 

Blank 

6 

29-42 

V 

V 

z 

7 

43-47 

Blank 

Blank 

Blank 

8 

48-61 

Blank 

Blank 

Blank 

9 

62  -66 

Blank 

Blank 

Blank 

10 

67-80 

Blank 

Blank 

Blank 

Figure  2-12.  IC0ND  Card  Format  No.  2 
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Z.  1.  3.  1.  3  ITIME.  This  name,  derived  from  Initial  Time,  identifies  the 


instant  in  time  when  the  initial  conditions  are  valid,  that  is,  the  epoch. 
Epoch  is  specified  by  year,  month,  day,  hours,  minutes,  and  seconds. 

Two  ITIME  cards  are  needed.  They  are  meaningless  alone  and  are  used 
only  to  accompany  an  ICTYP  1,  Z,  or  3  card  and  the  corresponding  IC0ND 
cards  to  complete  the  specification  of  the  input  vector.  The  format  of  the 
card  is  listed  below.  See  also  Section  Z.  Z  (Load  Sheets). 


Z-  30 


ITIME  No.  1 


ITIME  No.  1 


Field 

Column 

Content 

Description 

i 

1  -2 

01 

Sequence  number 

2 

3-4 

Blank 

Iteration  number 

3 

5-9 

ITIME 

Principal  variable  identifier 

4 

10-23 

Y 

Y  ear 

5 

24-28 

Blank 

6 

29-42 

M 

Month  number 

7 

43-47 

Blank 

8 

48-61 

D 

Day  number 

9 

62-66 

Blank 

10 

67-80 

H 

Hour  of  day 

Figure  2-13.  ITIME  Card  Format  Card  No.  1 
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ITIME  No.  Z 


ITIME  No.  2 


1 


hhs 

3  C  l  I  1 


8 


10 


0  0  0  0  0  0  0  0  0  0  0  0  0  op 

II  11  12  13  14  13  II  12  11  112121  22  I 

2222222222222: 

3  3  3  3  3  3  3  3  3  3  3  3  3  3R 
4444444444444 
5555555555555 
6  B  6  6  6  6  6  6  6  6  6  6  6  6  6g 
777777777777 
8  8  8  8  8  8  8  8  8  8  8  8  8  8  I 


9999993999999! 

10  11  12  13  14  13  it  II  il  II  20  21  22  2 


ro  0  0  0  0  0  0  0  0  0  0  0  0  c 

I2I  30  31  32  J3M  35  3I37  3I3I  40  41  4; 
111111111111111 

22222222222222 

3  3  3  3  3  3  3  3  3  3  3  3  3  3! 

4  4  4  4  4  4  4  4  4  4  4  4  4  4! 

55555555555555 

66666666666666 

77777777777777 


199999999999999 

12130)1  32  33  34  33  X  37  31  31  40  41  42 


>\\  W.  v\K\\Vis \n  A*  V  . . » 


Field 

Column 

Content 

1 

1-2 

02 

2 

3-4 

Blank 

3 

5-9 

ITIME 

4 

10-23 

Min 

5 

24-28 

Blank 

6 

29-42 

Sec 

7-10 

43-80 

Not  used 

Description 
Sequence  number 
Iteration  number 
Principal  variable  identifier 
Minute  s 

Seconds 


Figure  2-14.  ITIME  Card  Format  Card  No.  2 
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2 .  1.3.  1.4  SPADATS  7  -Card  Mean  Elements .  These  cards  are  standard 


to  the  SPADATS  system.  ESPOD  accepts  the  standard  7 -card  element  set 
as  input  (ICTYP  4).  The  shaded  columns  are  not  used.  The  formats  for 
these  cards  appear  on  the  following  pages. 
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SPADATS  -7 ,  No.  1 


SPADATS  -7 ,  No.  1 


/ 


1 


00  0 

1  2  3 

I  1  1 


4  4 


2  2 

3  3 

4  4 

5  5 
66 

7  7  7 

8  8  8 

9  9 
1  2  1 


I  I  I 


& 

0«0 


0  0  0_ 

.  5 

"I 

2  2  2  2^2 

3  3  3  3^3 

‘f 

5|5  5  5^ 

S  6  6  6^ 

I 

1 

8  8  8^1 
9j  9  9  9^9 


000000000  0^ 

I  ft  II  t2  13  14  IS  II  17  1|^ 

222222222  2$ 

> 

333333333  3$ 

5 

444444444  4% 

5555555555  5$ 

6666666666  6^ 

7  7  7  7  7  7  7  7  7  7  7$ 
88888880881 
9  9  9  9  9  9  9  9  9  9§$fc 

J  10  11  12  13  14  IS  It  17  UsV|kjr 


00000000000000 

24  25  X  27  7t  7t  30  31  J2  U  34  35  *| 
11111111111111 

J2  2222222222222 

J3  3333333333333 

^4  4444444444444 

$5  555555555555  5 

$6  6  6  5  6  6  t>  u  b  6  6  6  *  6 

\ 

S, 

$7  7  7  7  7  7  7  7  7  M  7  7  7 
8888888888868 
^99  =  9  9  9  9  9  9  8939, 

^3  24  2S  »  27  21  21  30  31  32  33  34  35  X 


m 


00000000000000 
37  31  3*  40  41  4?  43  44  45  4«  47  41  49  SO 
11111111111111 

22222222222222 

33333333333333 

44444444444444 

55555555555555 

66666666666666 

11111111111111 


99953999939999 

V  31  39  40  41  42  4  3  44  <tS  41  47  44  49  SO 


0000000000000 

si  S?  S3  S4  ss  SC  ST  SI  S9  U  11  12  13 

11111111111111 


10  11 


2222222222222 

3333333333333 

44444444444444 

5  5  5  5  5  5  5  5  5  5  5  5  5  5: 
66666666666666 
77777777777777 


99999999999999 

51  SJ  S3  ‘A  55  K  57  51  59  U  tl  V  U  (4 


Field  Column 

1  1-3 

2  4-6 


3  7 

4  8 

5  9-18 


6 

7 

8 

9 

10 

11 


19-22 

23-36 

37-50 

51-64 

65-79 

80 


Description 

Satellite  number — justified  right 

Element  set  number — justified 
right 

Not  used 

Card  number  (card  #=  -  1) 

Satellite  name  for  Element 
File  Update 

Not  used 

Nq  —  Epoch  revolution 
e  —  eccentricity 
i — inclination  (degrees) 

Not  used 

Card  type,  E  =  Nodal  Elements 


Figure  2-15.  SPADATS -7  Element  Card  Format,  Card  No.  1 
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SPADATS  -7 ,  No.  2 


SPADATS -7,  No.  2 


Field 

Column 

Description 

1 

1-3 

Satellite  number — justified  right 

2 

4-6 

Element  set  number 

3 

7 

Not  used 

4 

8 

Card  number  (Card  #  =  2) 

5 

9-12 

Year  of  T 

o 

6 

13-22 

Not  used 

7 

23-36 

T  — Time  of  Epoch  (day  of  year 
and  fraction  of  day) 

8 

37-40 

Not  used 

9 

41-50 

Not  used 

10 

51-64 

L^ —  Mean  Longitude — degrees 

11 

65-79 

Not  used 

12 

80 

Card  Type,  E  =  Nodal  Elements 

Figure  2-16.  SPADATS-7  Element  Card  Format,  Card  No.  1 
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SPADATS  -7 ,  No.  3 


SPADATS  -7,  No.  3 


/ 


ID 

i  i 

I  1 
3  3 

3  3 

4  4 

5  5 

6  § 

7  7 

8  8 
9  9 

i  i 


- ITT 

2  ,34 


Jon 

4  J  I  1 

mi 

3|3  2  2 

3  3  3  3 

4  4  4  4 

5  5  5  5 


I  0  0  11  0  0  0  0  0  0  0  0  0  0 

i  n  it  <i  u  m  13  it  n  ii  it  a  21  n 
1111111111111 


I” 

is  s 


6(8  6  6  ft 

J 

1?  7  7®  7 


| 

HU! 


9  9  9  9; 

4  S  *%  I 


0  0  0  0  0  0  0  0  0  0  0  0  0  0  3  0  0  0  0  0  0  0  0  0  0  0  0  0) 

23242321  27  2121)031  32  S3  >4  AN  DM3I4Q41  *2  4  j  44  IS  44  4  ’  41  *4  30 


72222222222222 
B3P  1113333333333 
M  4  4  4  4  4  4  4  4  4  4  4  4 
5  5  5  5  5  5  5  5  5  5  5  5  5 
S  6  6  6  6  6  6  6  S  6  6  b 
77777777777777 
S  8  S  8  e  8  8  8  8  8  8  8  8  8  8 


9  9 


9999999999993 

t  ID  I*  Lj  1)  14  It  II  u  II  It  70  21  23 


I  1 


111111111111111111111111111 
2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2fc 


0000000000000 

|1  57  S3  S4  »  M  1)  U  M  M  II  V  II  M 

I  1  1  1  1  1  I  I  1  1  I  1  I 


33333333333333333333333333333 


6  6 


44444444444444 

5  5  5  5  5  5  5  5  5  5  5  5  5  9 

b6bo666B6f'libL 
7  7  7  7  7  7  7  7  7  7  ’??? 

88888888888888 

9999999999^999 

:j  21  »  »  2WI  21  M  JI  K  U  M  a  * 


444444444444444 

5  5  5  5  5  5  5  5  5  5  5  5  5  5  5 
b  6  6  6  6  6  6  6  6  6  6  6  6  6  8 
7777777777/777 
888888888868888 


999999999999999 

3T  a  fe  40  41  42  43  44  45  4*  «2  *1  4t  X  51 


8 


’10 


2222222222222 

3333333333333 

4444444444444 

5555555555555 

66666666666666 

7/777777777777 

8888888888888 


9999999999999 

j2  U  S4  a  %  5?  si  M  <4  It  l?  13  |4 


Field  Colum  n 

1  1  -3 

2  4-6 


3  7 

4  8 

5  9-22 


6  23-36 


7  37-50 

8  51-64 

9  65-79 

10  80 


Des  c  ription 

Satellite  number — justified  right 

Element  set  numbe r —justified 
right 

Not  used 

Card  number  (Card  #  =  3) 

P  — Anomalistic  period  at  Epoch  — 
day s / rev 

$2  Q  —  Right  ascension  of  ascending 
node  degrees 

to0  —  Argument  of  perigee — degrees 
qo — Perigee  distance — Earth  radii 
Not  used 

Card  type,  E  =  Nodal  Elements 


Figure  2-17.  SPADATS-7  Element  Card  Format,  Card  No.  3 
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SPADATS  -7,  No.  4 


SPADATS  -7,  No.  4 


/ 


Tio 


9  9  9 

i  i  i 


2  34 


oooo 

1  3  3 
1  1  1 


2  2  2  2 

3  3  3 

4  4  4  4 

5  5  5 
666666 
1  1  1 
8  8  8 


IDOOOOOOOOOOOOOOO 

I  11  11  13  13  14  13  II  17  11  I)  W  31  22 
11111111111111 


22222222222222 
333333333333333 
444444444444444 
555555555555555 
616  6  6  6  6  6  6  6  6  6  6  6  6  6 
7777777777777 
888888888888888 


9  9  9 

4  S  I 


39  9  9  9  9  9  9  9  9  9  9  9  9  9  9 

3  10  11  13  13  14  13  II  17  II  11  30  31 


000000000000 
Z3  34  733l  77  7f?VV  31  33  33  14  35]l 
11111111111111 


0(0  0  0  0  0  0  0  0  0  0  0  0  0 
17  31  H  40  41  4?  43  44  44  41  47  41  41 
11111111111111 


22222222222222 

33333333333333 

4444444444444 

55555555555555 

66666666666666 

77777777777777 


1 


22222222222222 

33333333333333 

444444444444444 

55555555555555 

66666666666666 

77777777777777 


9999999999999999999999999999 

24  ft  21  ?7  ?l  21  B  11  32*1  M  »  X  r  31  }«  40  41  4?  O  «4  4S  41  4T  41  «  U 


v\  ^v\V\V'\wV\Ksv\\\  A  v\\\VvvW\\v,\V\v 


Field 

1 

1 

3 

4 


Column 

1-3 

4-6 

7 

8 


Description 

Satellite  number — justified  right 

Element  set  number — justified 
right 

Not  used 

Card  number  (Card  #  =  4) 


5 

6 

7 

8 


9-22  C  —  Rate  of  change  of  period  — 

days/ (rev)^ 

23-36  £2  — Time  derivative  of  Right 

ascension  of  ascending  node  — 
degrees/ day 

37-50  coG  —  Time  derivative  of  argument 

of  perigee — ( degree s / day ) 

5  1  -64  Not  used 


9  65-79 


Not  used 


10  80 


Card  type;  E  =  Nodal  Elements 


Figure  2-18.  SPADATS-7  Element  Card  Format,  Card  No.  4 
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SPADATS  -7 ,  No.  5 


SPADATS -7,  No.  5 


0  g  0  0  0  0  D  0  0  R  A  D  D  0  0 

♦  >i  ii  it  i)  it  is  it  >j  ii  it  » ii  27 

mninimii 
ill  222222  2222222 
333333333333333 

mUHUU444 
5ISI555S5SSSI55  5 
fij  6  b  6  fi  6  fi  £  6  B  B  &  B  S  5 
777777777777; 
dfi  B  B  fl  &  &  &  &  B  0  0  &  &  B 
9|S  9  9  9  9  9  9  9  $  9  9  9  9  9: 

9  1#  II  11  I J  1*  IS  II  17  Ifl  11  2V  21  21' 


iasaSS®® 


Field 

Column 

Des  c ription 

1 

1-3 

Satellite  number — justified  right 

2 

4  -6 

Element  set  number — justified 
right 

3 

7 

Not  used 

4 

8 

Card  number  (Card  #-  ~  5) 

5 

9-22 

d  —  decay  acceleration 

6 

23-79 

Not  used 

7 

80 

Card  type;  E  =  Nodal  Elements 

Figure  2-19.  SPADATS-7  Element  Card  Format,  Card  No.  5 


SPADATS- 7 ,  NO.  6 


SPADATS- 7 ,  NO.  6 


/ 


000 

1  2  3 
1 1 1 

2  2 

3  3  3 

4  4  4 

5  5  5 

6  6  6 
7  7 
888 


9  9  9 

1  I  3 


34 


oooSo 

4  3  1^1 


0(0  0  0  0  0  0  0  0  0  0  0  0  0  0  0 

I  II  11  12  13  14  IS  1«  17  II  II  2«  21  22 


3  3  3t 


7  7  7 


0000000000000 

23242S2I?7  2I213I313?  33  34  3S]| 
11111111111111 


2222222222222 
3333333333333333 
44444444444444 
555555555555555 
666666666666666 
7  7  7  7  7  7  7  7  7  7  7  7  7 
888888888888888 


99999999999999 

J  10  U  12  13  14  IS  II  17  II  II  20  21  22fc: 


2222222222222 

333333333333 

44444444444444 

55555555555555 

66666666666666 

777777777777? 

8888888888892 


99999999999999 

3  24  25  21  27  21  21  31  3«  12  33  >4  3S  31 


00000000000000 
37  31  31  41  41  42  43  44  45  4|  47  4|  4#  SO 
11111111111111 


222222222222222 

333333333333333 

44444444444444 

55555555555555 
6  S  6  6  6  6  6  6  6  6  6  6  6  6 
77777777777777 

988888888888888 


99999999999999 

3?  31  J!  4|  41  42  43  44  4S  4|  47  41  41  Ml 


wv\w  \S\\v\  VvVXvVVVVVXvXvANvvtt'wuv 


Field 

Column 

Description 

i 

1-3 

Satellite  number  -  justified  right 

z 

4-6 

Element  set  number  -  justified  right 

3 

7 

Not  used 

4 

8 

Card  number  (Card  #  6) 

5 

9-ZZ 

a  -  semimajor  axis  -  Earth  radii 

6 

Z3-36 

-  Nodal  period  -  days/ rev 

7 

37-50 

CNT  -  Rate  of  change  of  nodal  period 
day s / ( rev)^ 

8 

51-79 

Not  used 

9 

80 

Card  type;  E  =  Nodal  Elements 

Figure  Z-ZO. 


SPADATS-7  Element  Card  Format,  Card  No.  6 


SPADATS-7 ,  NO.  7 


SPADATS- 7  ,  NO.  7 


Field  Column  Description 


1 

1-3 

Satellite  number  -  justified  right 

2 

4-6 

Element  set  number  -  justified  right 

3 

7 

Not  used 

4 

8 

Card  number  (Card  $  -  7) 

5 

9-22 

Not  used 

6 

23-29 

Initial  revolution,  decimal  may  be  punched  in  column  29 

7 

30-36 

Final  revolution,  decimal  may  be  punched  in  column  36 

8 

37-50 

Expiration  date  of  Bulletin,  in  format;  YMMDDHHMMSS. 
decimal  punched  in  column  48 

9 

51-58 

RMS,  in  format  XXXXX.  XX;  decimal  punched  in  56 

10 

59-66 

Number  of  observations  used  in  obtaining  RMS 

11 

67 

IST0P 

Blank  or  0  =  correct  the  inclination  element 

1  =  do  not  correct  the  inclination 

2  =  do  not  correct  the  drag  parameter 

4  =  correct  time  equation  only 

12 

68 

Blank  or  0  =  mean  elements 

1  =  K-25  elements 

2  =  osculating  elements 

13 

69-79 

Not  used 

14 

80 

Card  type;  E  =  Nodal  Elements 

Figure  2-21.  SPADATS-7  Element  Card  Format,  Card  No.  7 
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2.  1.  3.  1.  5  DNREV.  This  function,  from  Delta  Number  of  Revolutions, 
controls  the  update  of  SPADATS  7  -  card  elements  from  the  epoch  given  with 
them  to  the  epoch  required  for  ESP0D.  The  DNREV  card  may  be  omitted 
if  automatic  program  logic  will  establish  the  desired  epoch.  It  also  con¬ 
trols  the  update  when  file  elements  (SEAI  tape)  are  used.  When  the  input 
vector  is  obtained  from  the  SEAI  tape,  the  epoch  is  automatically  chosen 
as  the  last  time  of  observation.  If  desired,  an  alternate  epoch  may  be 
selected  with  the  DNREV  card. 

The  first  variable  is  assigned  one  of  the  following  four  values:  "l.", 
n2,H,  n3.M  ,  or  n-l.n.  These  numbers  indicate  the  manner  in  which  the 
second  variable  is  interpreted. 
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DNREV 


DNREV 


Field  Content 


Desc  ription 


1 

Z 

3 

4 

5 

6 


01 

Blank 
DNREV 
X.  0 
Blank 

y 


x  (first  variable) 

1. 0 


Z.  0 


3.  0 

-  1.  0 


Sequence  Number 
Iteration  Number 
Principal  variable  identifier 
X  =  -1.0,  1.0,  Z.O  or  3.0 


Specification  of  new  epoch,  with  units  and  reference 
update  as  specified  in  tabulation  below. 

y  (second  variable) 

Time  in  days  (D)  and  fractions  of  days  (DF)  from 
first  of  year.  1.0  =  00  hrs,  universal  time,  on 
1  January  of  year. 

Number  of  revolutions  (REVS)  and  fractions  of  revo¬ 
lutions  (REVS  F)  from  epoch  given  with  SPADATS 
mean  elements  or  on  the  SEAI  tape,  to  epoch  required 
for  ESP0D. 

Revolutions  (REVS)  and  fraction  of  revolutions 
(REVS  F)  from  launch  (first  ascending  node  =  REV  1.0) 
to  epoch  desired  for  ESP 

Variable  fields  are  left  blank.  This  specifies  the  up¬ 
dated  epoch  when  SPADATS,  mean  elements  are  used 
to  the  time  of  last  observations. 


If  no  automatic  update  to  the  time  of  last  observation  is  desired  when 
using  file  elements  (SEAI  tape),  a  value  of  X  =  Z.  0  and  y  =  0.  0  is  selected. 


FIELDS  7,  8,  9  and  10  are  not  used. 


Figure  Z-ZZ.  DNREV  Card  Format 
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2.  1.3.  2  Specification  of  Solution  Variables 


Solution  variables  are  specified  with  two  different  cards  having 
different  formats.  These  cards  set  up  solutions  for  the  orbit  elements 
and  the  sensor  parameters  respectively. 

2.  1 .  3.  2.  1  CAT1 — Category  1  Variables 


Sy  mbol 

Name 

Cor re  s ponding 
Card  Column 

a 

o 

Right  ascension 

10 

6 

Declination 

11 

o 

Po 

Flight  path  angle 

12 

A 

Azimuth 

13 

o 

R 

Radius 

14 

o 

V 

o 

Velocity 

15 

C^A/ 2m 

Drag  parameter 

16 

K 

Drag  parameter  variation 

17 

The  first  eight  columns  of  the  V-Field-1  in  the  CAT  1  card  corre¬ 
spond  in  the  order  listed  above  to  the  Category  1  variables  to  be  solved 
for.  The  V-Field  columns  must  contain  either  "ones"  or  "zeros."  A  11 1 n 
punched  in  a  column  of  the  V-Field  indicates  that  the  corresponding  varia¬ 
ble  is  to  be  solved  for.  A  n0n  punched  in  the  same  column  indicates  that 
the  corresponding  variable  is  not  to  be  solved  for.  This  convention  holds 
true  for  eight  Category  1  variables.  For  example,  a  card  punched  as 
follows 


0J1 

L 

C  |A|  T,  1  | 

1|1|0|1|1|°|0|0,  I  1  I  I  1 

- V-Field-1 - - 

indicates  that  on  this  run  the  variables  a  ,  6  ,  A  ,  and  R  are  to  be  solved 

o  o  o  o 

and  that  6  ,  V  ,  C^A/Zm,  and  K  will  be  held  constant.  Note  that  the  six 
ro  o  D 

righthand  columns  of  the  V-Field-1  are  left  blank. 
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CAT1 


CAT1 


Field 

Column 

Content 

Description 

1 

1-2 

01 

Sequence  number 

2 

3-4 

Blank 

Iteration  number 

3 

5-9 

CAT1 

Principal  variable  identifier 

4 

10-17 

XXXXXXXX 

X  =  0  variable  not  to  be  solved 

X  =  1  variable  to  be  solved 

4 

18-23 

Blank 

5-10 

24-80 

Not  used 

Figure  2-23.  CAT1  Card  Format 
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2.  1.3.  2. 2  CAT2  —  Category  Z  Variables 


Symbol 

Name 

Cor  re  s ponding 
Card  Column 

Rb 

Range  bias 

14 

33 

52 

71 

Ab 

Azimuth  bias 

15 

34 

53 

72 

Eb 

Elevation  bias 

16 

35 

54 

73 

Rb 

Range  Rate  bias 

17 

36 

55 

74 

RA, 

b 

Right  Ascension  bias 

18 

37 

56 

75 

Db 

Declination  bias 

19 

38 

57 

76 

Time  bias 

20 

39 

58 

77 

*b 

Latitude  bias 

21 

40 

59 

78 

Longitude  bias 

22 

41 

60 

79 

hb 

Height  bias 

23 

42 

61 

80 

Category  Z  variables  are  solved  for  selected  sensors.  Any  combina¬ 
tion  of  CATZ  variables  may  be  solved  for  the  same  sensor.  The  convention 
for  specifying  CATZ  variables  to  be  solved  is  similar  to  that  used  in  CAT1 
cards. 

The  14  columns  of  the  V  -  fields  are  partitioned  into  two  subfields  of 
four  columns  and  ten  columns  respectively. 

•  The  first  4  columns  contain  the  sensor  number  right  adjusted 

•  The  last  10  columns  correspond  in  the  order  listed  above  to 
the  Category  Z  variables  to  be  solved  for. 

For  example,  the  card  punched  as  follows  indicates  by  means  of  the 

contents  of  the  V-Field-1  that  R.  ,  A,  ,  E.,  and  h.  of  sensor  3Z4  are  to  be 

b  b  b  b 

solved  for. 


Vi  i 


CAT  2 

0,3,2,41 ,1 .10.0,0,0,0,0,1 

-1,11  1 

1  1  1  1  1  1  1  1  1  1  1  1  1 

1  1  1  1 

I  I  1 1  I  1  I 
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If  a  CAT2  card  is  omitted,  no  Category  2  variables  are  solved  for. 
Unlike  the  Category  1  variables,  no  initial  estimates  of  the  Category  Z 
variables  are  needed,  although  initial  bias  estimates  can  be  specified  (see 
Section  2.  1.  3.2.4). 

If  additional  station  biases  (CAT2  variables)  are  to  be  solved,  more 
CAT2  cards  must  be  added.  The  only  format  difference  from  card  to  card 
is  the  first  field  location,  the  sequence  number. 
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CATZ 


CATZ 


2  2  2  2 


9  9  9  9 


0000000000 

15  16  17  IH  19  20  21  22  23 

1111111111 


222222222 

3333333333 

4444444444 

5555555555 

6666666666 

7777777777 


9999999999 


0  11  12  13  H  15  16  17  II  19  21  21  22  K 


_8_ 

000000000 

>2  53  54  65  56  57  54  59  60  61 

111111111 

222222222 

333333333 

444444444 

555555555 

6666666666 

777777777 


999999999 

52  53  54  55  56  57  51  59  60  61 


1C 

0000000000 

71  7?  73  74  75  71  77  71  79  10 

1111111111 

222222222 

333333333 

444444444 

555555555 

666666666 

7777777777 


99999999999999 


71  7?  73  74  75  76  77  76  13  Li 


Field 

Column 

Content 

Description 

1 

l-Z 

01 

Sequence  number 

2 

3-4 

Blank 

Iteration  number 

3 

5-9 

CAT2 

Principal  variable  identifier 

4 

10-13 

XXXX 

First  sensor  number  (right  adjusted) 

4 

14-23 

XXXX 

X  =  0  variable  not  to  be  solved  for 

X  =  1  variable  to  be  solved  for 

5 

24-28 

Blank 

6 

29-32 

XXXX 

Second  sensor  number 

6 

33-42 

XXXX 

X  =  0  variable  not  to  be  solved  for 

X  =  1  variable  to  be  solved  for 

7 

43-47 

Blank 

8 

48-51 

XXXX 

Third  sensor  number 

8 

52-61 

XXXX 

X  =  0  variable  not  to  be  solved  for 

X  =  1  variable  to  be  solved  for 

9 

62-66 

Blank 

10 

67-70 

XXXX 

Fourth  sensor  number 

10 

71-80 

XXXX 

X  =  0  variable  not  to  be  solved  for 

X  =  1  variable  to  be  solved  for 

Figure  Z-Z4.  CATZ  Card  Format 


Z  -55 


2.  1  .  3.  Z.  3  Bounds  on  Variable  Changes.  BNDS  (abbreviation  from  bounds) 
identifies  the  start  of  an  array  which  contains  the  bounds  to  be  used  in  the 
differential  correction  process.  Nominal  values  for  bounds  are  stored  in 
the  program,  but  if  any  of  these  is  to  be  modified,  a  new  set  of  bounds 
must  be  provided.  For  each  variable  selected  for  solution  by  either  a 
Category  1  or  Category  2  card,  a  bound  on  the  size  of  correction  is  speci¬ 
fied  to  control  convergence.  The  variables  selected  for  solution  are  placed 
in  the  same  order  as  the  columns  of  the  Category  1  and  Category  2  cards. 
Multiple  Category  2  cards  are  taken  in  their  sequential  order.  The  total 
count  of  variables  tagged  with  a  11 1 11  is  the  dimension  of  the  solution  vector. 
A  set  of  bounds,  corresponding  one-to-one  with  the  elements  of  the  solution 
vector,  is  listed  in  the  same  order  starting  in  the  first  variable  field  on  the 
bounds  card  and  continuing  with  no  intervening  blank  variable  fields  across 
the  first  card  and  through  any  necessary  following  cards  until  the  proper 
number  of  bounds  is  specified. 

Figures  2-25  and  2-26  show  variable  field  assignments  for  the  bounds 
which  correspond  to  the  solution  vector  of  the  C AT  1  and  CAT2  variables 
given  in  the  previous  section. 
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BNDS  NO.  1 


BNDS  NO.  1 


Field 

Column 

Content 

Description 

1 

1-2 

01 

Sequence  number 

2 

3-4 

Blank 

Iteration  number 

3 

5-9 

BNDS 

Principal  variable  identifier 

4 

10-23 

B 

a 

Bound  on  right  ascension 

5 

24-28 

Blank 

6 

29-42 

B6 

Bound  on  declination 

7 

43-47 

Blank 

8 

48-61 

ba 

Bound  on  azimuth 

9 

62-66 

Blank 

10 

67-80 

br 

Bound  on  radius 

Figure  2-25,  BNDS  Card  Format,  Card  No.  1 
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BNDS  NO.  2 


BNDS  NO.  2 


10 

i  i 
1  1 

1  1 

J  3 

4  4 

*5 
6  6 
7  7 
B  B 
9  9 

t  2 


9  9  9  9 

i  m 


JO  0  Q  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0 

II  I  I  t  II  11  12  13  U  13  II  17  II  II  JO  71  n 
1111111111111111111 

222222222222222222 

.3  333133333333333333 

4444444444444444444 

595555555555555555 

6666666666666666666 

7777777777777777777 

1866588888888888888 


99999999999999 

10  11  12  13  14  13  II  17  II  II  20  21  22 


10 


0000000000000 
17WII70  71  77  73  74  7I  70  77  7I7III 

1111111111111 

2222222222222 

3333333333333 

4444444444444 

55555555555555 

666666666666S 

777/7777777777 

88888888888888 

99999999999999 

H  It  70  71  77  >]  74  73  71  77  7|  7|  M 


Field 

Column 

Content 

i 

1-2 

02 

2 

3-4 

Blank 

3 

5-9 

BNDS 

4 

10-23 

br 

Kb 

5 

24-28 

Blank 

6 

29-42 

B\ 

7 

43-47 

Blank 

8 

48-61 

bf 

Eb 

9 

62-66 

Blank 

10 

67-80 

X 

Desc  ription 
Sequence  number 
Iteration  number 
Principal  variable  identifier 
Bound  on  range  bias 

Bound  on  azimuth  bias 

Bound  on  elevation  bias 

Bound  on  height  bias 


Figure  2-26.  BNDS  Card  Format,  Card  No.  2 
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Z.  1.  3.  Z.4  BISES.  This  name,  derived  from  BIAS  ESTIMATE,  identifies 
the  start  of  an  array  which  contains  the  initial  bias  estimates  of  the  Cate¬ 
gory  Z  variables  in  the  differential  correction  solution  vector.  They  are 
specified  in  the  same  order  as  that  for  bounds .  (See  Section  Z .  1 .  3 .  Z .  3.  ) 

If  any  Category  Z  initial  estimate  is  specified,  all  must  be  specified.  If 
a  particular  initial  estimate  is  not  known,  zero  (0.0)  may  be  entered.  The 
first  four  initial  estimates  of  the  set  are  entered  in  V-Field  1-4  of  the  first 
BISES  card.  Additional  bias  estimates,  if  any,  are  arranged  in  successive 
cards  with  no  intervening  blank  V-Fields. 

The  following  example  illustrates  the  format  of  the  BISES  card.  The 
bias  estimates  correspond  to  the  Category  Z  variables  discussed  in  a  pre¬ 
vious  section. 

Category  Z  initial  estimates  (and  later  improved  values)  are  alge¬ 
braically  subtracted  from  the  measured  observations  in  the  bias  correction 
portion  of  ESP0D.  If  it  is  desired  to  remove  a  known  bias,  and  not  to 
allow  it  to  change,  the  following  procedure  can  be  used: 

•  Specify  a  bias  estimate  with  BISES 

•  Select  the  appropriate  Category  Z  variable  in  the  solution 
vector  with  CATZ 

•  Set  the  corresponding  bound  to  zero  (0.  0)  with  the  BNDS 
input. 
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BISES 


BISES 


t9  9  9  9  9  9 

fcs  I  7  •  I  II 


10000000000001 

I  11  1?  13  14  IS  H  13  II  It  »  ?1  22  2 

1 1  M  1 1 1 1 1 1 1 1  1 


222222222222  2S 
3333333333333  3g 
4444444444444 
555555555555  5§ 
66  6^6  666666666  6q 
7777777777777  7^ 
8888888888888 8| 
999999999999  9R 


^0  0000000000000] 

(2*  JO  31  U  O  34  »  31  37  a  »  40  41  42  j 
\\  11  11  1  11  1  1  11  1  1 

b 2222222222222; 

S3  3333333333333; 
^4  4444444444444; 
:5  5555555555555) 
b6  6666666666666) 
?7  7777777777777) 
Jo  8888888888888) 
19  9  9  9  9  9  9  9  9  9  I  9  9  9! 

30  31  32  U  34  33  31  37  31  31  40  41  42  | 


so  0  0  0  0  0  0  0  0  0  0  0  0  of 

f40  41  U  si  H  S3  »4  »  SO  19  *  Si  M  111 

Si  1  1  11  1  1  11  11  1  1  if 
[22222222222222 

S3  3333333333333 
S  4  4444444444444 
555555555555555 
^6  6666666666666 
‘77777777777777 
;  8  8  8  8  8  8  8  8  8  8  8  8  8  8 
1999999999999 

41  41  50  51  S2  S3  54  55  SI  57  SI  »  M  II 


goooooooooooooo 

|l7HM70  7inn74TS7l717l71N 
^11111111111111 

[22222222222222 

f33333333333333 

^44444444444444 

[55555555555555 

[66666666666666 

77777777777777 

88888888888888 

99999999999999 

|ir  M  II  70  71  7?  ^3  74  75  71  77  7|  71  I 


F  ield 

Column 

Content 

De  sc  ription 

i 

1  -2 

01 

Sequence  number 

2 

3-4 

Blank 

Iteration  number 

3 

5-9 

BISES 

Principal  variable  identifier 

4 

10-23 

Estimated  R, 
b 

Range  bias  estimate 

5 

24-28 

Blank 

6 

29-42 

Estimated  A, 
b 

Azimuth  bias  estimate 

7 

43-47 

Blank 

8 

48-61 

Estimated  E. 

b 

Elevation  bias  estimate 

9 

62-66 

Blank 

10 

67-80 

Estimated  h, 
b 

Height  bias  estimate 

Figure  2-27.  BISES  Card  Format 
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2.  1.  3.  3  Convergence  Control  —  NITER 


ESPOD  will  automatically  come  to  a  stop  if  the  solution  either  con¬ 
verges  or  diverges  according  to  internal  criteria.  Under  certain  cases  of 
slow  convergence  or  instability  ESPOD  could  run  an  undesirably  long  time. 
For  this  reason,  it  is  necessary  to  limit  the  number  of  iterations  which 
the  program  will  be  permitted  to  try.  This  is  accomplished  automatically 
with  an  internally  contained  value  (5),  but  the  value  may  be  changed  with  a 
NITER  (N  iterations)  card  (see  Figure  2-28). 
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NITER 


NITER 


000000000000 

it  ii  i2  i)  14  ii  it  ii  it  i»  n  2i  n 

11111111111111 


22222222222222 
3333333333333 
4  4  4  4  4  4  4  4  4  4  4  4  4  4j 
5  5  5  5  5  5  5  5  5  5  5  5  5  5 
6666666666666 
77777777777777 
8888868888888 
99999999999999 

10  ll  12  I)  14  II  II  II  II  II  71  21  22  2: 


F  leld 

Column 

Content 

i 

1  -2 

01 

2 

3-4 

Blank 

3 

5-9 

NITER 

4 

10-23 

N. 

5-10 

24-80 

Not  used 

Description 
Sequence  number 
Iteration  number 
Principal  variable  identifier 
Maximum  iterations  permitted 


Figure  2-28.  NITER  Card  Format 
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2.  1.3.4  Integration  Control 


ESP0D  assumes  an  initial  step  size  of  numerical  integration  of  one 
minute  unless  a  TSTEP  (Time  STEP)  card  is  provided  (see  Figure  2-29). 
Any  other  value  of  initial  step  size  equal  to  (P)(2)  (in  minutes)  where  P 
and  m  are  integers  which  may  be  specified.  The  maximum  initial  step 
size  allowed  is  128  minutes.  The  initial  step  size  specifies  the  integration 
step  at  the  beginning  of  the  integration.  ESP0D  internal  logic  varies  the 
step  size  to  optimize  the  integration  process.  The  maximum  and  minimum 
integration  time  step  limits,  and  other  constants  in  the  integration  process 
can  be  changed  with  a  99-card  input  (see  Section  5.  7  of  ESD- TDR  -  64- 395) . 
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TSTEP 


TSTEP 


Field 

Column 

Content 

1 

1-2 

01 

2 

3-4 

Blank 

3 

5-9 

TSTEP 

4 

10-23 

X.  0 

5-10 

24-80 

Not  used 

Desc  r  iption 
Sequence  number 
Iteration  number 
Principal  variable  identifier 
Nominal  step  size,  in  minutes 


I  igure  2-29.  TSTEP  Card  Format 


2.  1.  3.  5  Data  Editing 


This  section  is  concerned  with  editing  residuals  and  rejecting  those 
which  are  not  to  be  included  in  the  differential  correction.  Items  may  be 
manually  identified  for  rejection,  and/or  they  may  be  rejected  by  com¬ 
parison  with  a  gross  outlier  test,  and/or  they  may  be  rejected  by  compari¬ 
son  with  some  multiple  K  of  the  root  mean  square  of  the  weighted  residuals 
from  the  previous  iteration,  and/or  they  may  be  rejected  if  the  correspond¬ 
ing  observations  are  removed  more  than  some  arbitrary  time  from  epoch. 
The  following  inputs  are  provided  to  correspond  to  the  above  options: 

a.  DELET  card 

b.  STYPE  card  to  change  gross  outlier  criterion  and/or 
r  ef  ractivity 

c.  99- card  item  452  to  change  K 

d.  TMAX  to  limit  time  from  epoch 

Residuals  are  weighted  by  the  reciprocal  of  the  standard  deviation 
associated  with  the  observation  type  and  the  sensor  class.  The  sensor 
class  is  specified  in  a  table  stored  in  the  program  which  can  be  changed 
with  a  STYPE  card,  item  b.  The  standard  deviations  are  also  in  a  table 
stored  in  the  program  which  can  be  changed  with  a  SIGMA  card  (see  Section 
5.  3  of  ESD- TDR  -  64  -  3 95,  f,ESPOD  Mathmatical  and  Subroutine  Description11 ) . 

2.  1.  3.  5.  1  DELET.  DELET  marks  the  beginning  of  a  table  of  identifica¬ 
tion  numbers  of  residuals  to  be  deleted.  The  identification  numbers  are 
taken  from  the  residuals  listing  output  after  running  the  computer  at  least 
once  with  the  case  in  hand  (see  Section  4.  2.  5).  V-Fields  in  DELET  cards 
are  specified  in  pairs  (a,  b)  which  mark  respectively  the  first  and  last 
residual  of  a  list  of  residuals  to  be  deleted.  The  program  deletes  all 
residuals  whose  identification  number  lies  between  a  and  b  inclusive.  If 
a  =  b,  a  single  residual  is  deleted.  Up  to  100  lists  may  be  deleted. 
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DELET 


DELET 


Field 

Column 

Content 

Description 

i 

1-2 

01 

Sequence  number  (max  =  25) 

2 

3-4 

Blank 

Iteration  number 

3 

5-9 

DELET 

Principal  variable  identifier 

4 

10-16 

al 

First  residual  of  first 
list  to  be  deleted 

4 

17-23 

bl 

Bast  residual  of  first 
list  to  be  deleted 

5 

24-28 

Blank 

6 

29-35 

a2 

First  residual  of  second 
list  to  be  deleted 

6 

36-42 

b2 

Last  residual  of  second 
list  to  be  deleted 

7 

43-47 

Blank 

8 

48-54 

a3 

First  residual  of  third 
list  to  be  deleted 

8 

55-61 

b3 

Last  residual  of  third 
list  to  be  deleted 

9 

62-66 

Blank 

10 

67-73 

a4 

First  residual  of  fourth 
list  to  be  deleted 

10 

74-80 

b4 

Last  residual  of  fourth 
list  to  be  deleted 

Figure  2-30. 

DELET  Card  Format 
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2.  1.  3.  5.  2  STYPE.  Sensor  characteristics  contained  in  the  sensor  table 


can  be  changed  with  an  STYPE  card.  The  following  parameters  are  con 
tained  in  the  sensor  table  identified  by  sensor  numbers  (corresponding 


card): 

Symbol 

Name 

cr 

t 

Sigma  type 

G 

s 

Gross  outliers  criterion 

R 

Refraction  flag 

N 

s 

Mean  surface  value  of  refractivity 

Format  of  the  STYPE  card  is  shown  in  Figure  2-31 


2-68 


STYPE 


STYPE 


F  ield 

Column 

Content 

Description 

i 

1  -2 

01 

Sequence  number 

2 

3-4 

Blank 

Iteration  number 

3 

5-9 

STYPE 

Principal  variable  identifier 

4 

10-13 

No. 

Sensor  number 

4 

14-1  5 

aT 

Sigma  type 

4 

16-19 

G 

s 

Gross  outliers  criterion 

4 

20 

RF 

Refraction  flag 

0  =  no  refraction  correction 

1  =  refraction  correction 

4 

21  -23 

N 

s 

Mean  surface  value  of 
refractivity 

5-10 

24-80 

Not  used 

Figure  2-31.  STYPE  Card  Format 
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2.  1.  3.  5.  3  TMAX.  ESPOD  will  automatically  reject  observations  which 
are  more  than  ten  days  from  epoch,  This  criterion  interval  can  be  changed 
with  a  TMAX  (TIME  Maximum)  card. 
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TMAX 


TMAX 


1 


10 


0000000000001 

11  i]  ii  m  is  ii  ii  it  it  m  it  a  a 
11111111111111 

222222222222 

3  3  3  3  3  3  3  3  3  3  3  3  3; 

44444444444444 

5555555555555S 

6  6  6  6  6  6  6  6  6  6  6  6  6  6! 

77777777777777 

688888888888888 

9  9  9  9  9  9  9  9  9  9  9  9  9  9jj 

<0  11  1?  13  14  IS  II  17  II  II  21  J1  22  nk 


Field 

Column 

Content 

1 

1-2 

01 

2 

3-4 

Blank 

3 

5-9 

TMAX 

4 

10-23 

N.  XX 

5-10 

24-80 

Not  used 

Description 

Sequence  number 

Iteration  number 

Principal  variable  identifier 

Integration  limit  from  epoch 
in  days 


Figure  2-32.  TMAX  Card  Format 
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2.  1.  3,  5.  4  SIGMA.  ESPOD  always  performs  a  weighted  differential  cor¬ 
rection.  Weights  are  assigned  to  observations  from  a  row  of  the  standard 
deviation  (SIGMA)  table.  The  row  is  selected  appropriate  to  a  particular 
sensor.  The  row  number  is  the  "sigma  type"  (see  field  4  of  STYPE  card, 
Section  2.  1.  3.  5.  2)  recorded  in  the  sensor  table.  There  are  60  rows  in  the 
SIGMA  table.  Each  row  contains  four  values  which  are  interpreted  dif¬ 
ferently  for  radar  measurements  than  for  camera  measurements. 


Table  2-1.  SIGMA  Table 


Value 

1 


V  alue 

2 


V  alue 
3 


V  alue 
4 


RADAR  <r  cr^ 

Banke  r 

Nunn  o'  (FR)  <rA  (FR) 

^  CL  O 

Game  ra 

(FR)  Field  Reduced  Data 
(PR)  Precision  Reduced  Data 


c r 


R 


o-  (PR)  <r  (PR) 

a  o 


Any  line  entry  of  the  SIGMA  table  can  be  changed  with  a  SIGMA  card.  As 
many  may  be  changed  as  necessary.  Nominal  values  stored  in  the  program 
may  be  printed  for  review  on  option.  See  Figure  2-33  for  field  layout. 


2-72 


SIGMA 


SIGMA 


F  ield 

Column 

Content 

Description 

i 

1-2 

NN 

Sensor  type 

2 

3-4 

Blank 

Iteration  number 

3 

5-9 

SIGMA 

Principal  variable  identifier 

4 

10-23 

aR 

Radar:  Standard  deviation  in  range 

<ra  (FR) 

Camera:  Standard  deviation  in  right 
ascension,  field  reduced 

5 

24-28 

Blank 

6 

29-42 

Radar:  Standard  deviation  in  azimuth 

a6  (FR) 

Camera:  Standard  deviation  in 

declination,  field  reduced 

7 

43-47 

Blank 

8 

48-61 

"e 

Radar:  Standard  deviation  in  elevation 

<Ta  (PR) 

Camera:  Standard  deviation  in  right 

ascension,  precision  reduced 

9 

62-66 

Blank 

10 

67-80 

Radar:  Standard  deviation  in  range 
rate 

cr6  (PR) 

Camera:  Standard  deviation  in  decli¬ 
nation,  precision  reduced 

Figure  2-33. 

SIGMA  Card  Format 
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2.  1.  3.  6  Atmospheric  Drag 


ESPOD  provides  on  option  three  drag  models  and  four  atmosphere 
models.  Two  of  the  atmosphere  models  require  the  specification  of  the 
parameters  and  F^.  Refer  to  Section  5.  10  of  ESD-TDR-64-395, 

"ESPOD  Mathematical  and  Subroutine  Description,"  for  full  information 
regarding  the  use  of  the  atmospheric  drag  options. 

2.  1,  3.  6.  1  Drag.  To  include  atmospheric  drag  in  the  Category  1  variables 
a  DRAG  card  must  be  input.  The  following  items  are  specified  on  the  DRAG 
card: 

Symbol  Description 

C^A/2m  (=£  0)  Drag  parameter  estimate 

K(^  0)  Drag  variation  options: 

No  variation 

1.  Secular  variation 

2.  Periodic  variation 

Model  Atmosphere  Options: 

Blank"] 

0.  >COESA  62  Static 

1.  ARDC  59 

2.  Paetzold  dynamic 
4.  COESA  dynamic 


Blank^ 

°-  J 
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DRAG 


DRAG 


Field 

Column 

Content 

Description 

1 

1  -2 

01 

Sequence  number 

2 

3-4 

Blank 

Iteration  number 

3 

5-9 

DRAG 

Principal  variable  identifier 

4 

10-23 

C pA/  Zm 

Drag  parameter  estimate  0) 

5 

24-28 

Blank 

6 

29-42 

K 

Drag  parameter  variation 
estimate  0) 

7 

43-47 

Blank 

8 

48-61 

1.  or  Z. 

Drag  model  option 

9 

62-66 

Blank 

10 

67-80 

1.,  2.,  3,, 
or  4. 

Atmospheric  model  options 

Figure  2-34.  DRAG  Card  Format 
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2  .  1 . 3 . 6.  2  APF 1  0 .  When  a  dynamic  atmosphere  is  chosen  on  the  DRAG 

card,  an  APF10  card(s)  must  be  supplied.  A^  is  the  daily  geomagnetic 

planetary  amplitude  or  an  approximation  to  it.  Fjq  is  the  10.7  centimeter 

solar  flux.  See  Section  5.  10.4  of  ESD- TDR- 64- 39 5,  "ESPOD  Mathematical 

and  Subroutine  Desc ription,  "  for  sources  of  A  and  Fir,. 

r  p  1 0 


A  single  value  of  A^  and  F^q  specified  only  for  the  day  of  epoch  indi¬ 
cates  that  it  is  intended  as  an  average  value,  and  ESPOD  will  use  it  for 
each  day  covered  by  the  curve  fit  or  trajectory.  Values  of  A^  and  P^q  may 
be  specified  separately  for  each  day  of  integration. 
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APF10,  No.  1 


APF10,  No.  1 


Field 

Column 

Content 

Description 

i 

1-2 

01 

Sequence  number 

2 

3-4 

Blank 

Iteration  number 

3 

5-9 

APF10 

Principal  variable  identifier 

4 

10-13 

NNN. 

Day  of  year  corresponding 
to  values  of  Ap,  Flo. 

4 

14-18 

A 

P 

Geomagnetic  planetary 
amplitude 

4 

19-23 

F10 

10.  7  cm  solar  flux 

5 

24-28 

Blank 

6 

29-32 

NNN. 

Day  of  year 

6 

33-37 

A 

P 

6 

38-42 

F10 

7 

43-47 

Blank 

8 

48-51 

NNN. 

8 

52-56 

A 

P 

8 

57-61 

F10 

9 

62-66 

Blank 

10 

67-70 

NNN. 

10 

71  -75 

A 

P 

10 

76-80 

F1  0 

Figure  2-35.  APF1 0  Card  Format,  Card  No.  1 
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APF10,  No.  2 


APF10,  No.  2 


l  9  9  9  9  9  9  9  9  9  9  9  9  9  9  91 

I  ’0  II  12  IJ  U  IS  II  IT  II  II  n  21  22  n\ 


/s%y/ss/s///ssA'//*///*s/ss%Ss//s/s/y.'/**SS/*s/*/s*/si*s/*s''//*/s**.'/*s/sS'/JtosysS/>//ss*s'' 


Field 

Column 

Content 

Description 

1 

1-2 

01 

Sequence  number 

2 

3-4 

Blank 

Iteration  number 

3 

5-9 

APF10 

Principal  variable 
identifier 

4 

10-13 

NNN. 

Day  of  epoch 

4 

14-18 

A 

P 

Geomagnetic  planetary 
amplitude, average  value 

4 

19-23 

F10 

10.  7  cm  solar  flux, 
average  value 

5-10 

24-80 

Not 

used 

Figure  2-36.  APF10  Card  Format,  Card  No.  2 
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2 .  1.3.7  Earth  Potential  Model 


ESPOD  automatically  provides  a  simple  Earth  potential  model  in  the 
absence  of  any  specified  input  (see  Model  1  in  Table  2-II).  For  different 
or  more  complicated  models,  modifications  to  the  stored  parameters  and 
callouts  for  higher  potential  harmonics  are  required.  See  Section  5.8  of 
ESD- TDR-64-395,  "ESPOD  Mathematical  and  Subroutine  Description,"  for 
detailed  discussion  of  the  Earth's  gravitational  potential  model. 

Modifications  to  the  stored  constants  are  made  by  means  of  "99  " 
cards  (see  Section  2.  1.3.  11).  Callout  of  higher  potential  harmonics  are 
made  by  means  of  the  ZONAL,  SECTR,  and  TESSR  cards.  For  the  ZONAL 
and  SECTR  cards,  a  "1"  in  a  particular  column  indicates  that  the  corre¬ 
sponding  coefficient  of  the  spherical  harmonic  is  to  be  included.  See 
Table  2-II  for  card  column  assignments.  In  the  TESSR  card,  the  order  of 
callout  is  not  important.  A  "42"  in  an  allowed  subfield  indicates  that  the 
^  tesseral  harmonic  is  to  be  included.  The  range  of  coefficients  of  the 
tesseral  harmonics  permissible  (14  in  all)  is  given  in  Table  2-IL 

Table  2-II  is  a  summary  of  the  potential  models  and  the  respective 
harmonic  coefficients.  Model  1  is  the  nominal  model;  that  is,  it  is  the 
model  which  ESPOD  will  use  when  none  other  is  specified.  Note  that  the 
coefficients  of  the  zonal  harmonics  ^h>  anc^  ^ \z  are  nom^na^y  zero 

in  the  program.  The  same  is  true  for  the  coefficients  of  the  sectorial  and 
tesseral  harmonics  of  order  5  and  6.  Values  for  these  coefficients,  whose 
core  locations  are  indicated,  may  be  input  on  99-cards. 
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Table  2 -II.  Earth  Potential  Model  Options  and  Nominal  Values  of  Harmonics 
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Table  2-II.  Earth  Potential  Model  Options  and  Nominal  Values  of  Harmonics  (Continued) 
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These  values  must  be  specially  inserted  by  99-card  input  when  using  Models  2,  3,  and  4. 
See  Section  5.0,  ESPOD  Operating  Instructions  and  Card  Formats. 


Z0NAL 


Z0NAL 


Field 

Column 

Content 

Desc  nption 

i 

1-2 

01 

Sequence  number 

2 

3-4 

Blank 

Iteration  number 

3 

5-9 

Z0NAL 

Principal  variable  identifier 

4 

10 

0  or  1 

J2 

0  = 

delete 

L. 

1  = 

include 

4 

1 1 

0  or  1 

^3 

0  = 

delete 

1  = 

include 

4 

12 

0  or  1 

J4 

0  = 

delete 

1  = 

include 

4 

13 

0  or  1 

J5 

0  = 

delete 

D 

1  = 

include 

4 

14 

0  or  1 

J6 

0  = 

delete 

1  = 

include 

4 

15 

0  or  1 

J7 

0  = 

delete 

1  = 

include 

4 

16 

0  or  1 

J8 

0  = 

delete 

1  = 

include 

4 

17 

0  or  1 

J9 

0  = 

delete 

1  = 

include 

4 

18 

0  or  1 

J10 

0  = 

delete 

1  = 

include 

4 

19 

0  or  1 

J11 

0  = 

delete 

1  = 

include 

4 

20 

0  or  1 

J 1  2 

0  = 

delete 

1  = 

include 

4 

21-23 

Not  used 

5-10 

24-80 

Not  used 

F  igure 

2-37,  Z0NAL  Card  Format 
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SECTR 


SECTR 


Field 

Column 

Content 

Desc  r  iption 

i 

1-2 

01 

Sequence  number 

2 

3-4 

Blank 

Iteration  number 

3 

5-9 

SECTR 

Principal  variable  identifier 

4 

10 

0,  1 

J2,  2 

0  =  delete 

1  =  include 

4 

1 1 

0,  1 

J3,  3 

0  =  delete 

1  =  include 

4 

12 

0,  1 

J4,  4 

0  =  delete 

1  =  include 

4 

13 

0,  1 

J5,  5 

0  =  delete 

1  =  include 

4 

14 

0,  1 

J6,  6 

0  =  delete 

1  =  include 

4 

15-23 

Not 

used 

5-10 

24-80 

Not 

used 

Figure  2-38.  SECTR  Card  Format 
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TESSR 


TESSR 


Field 

Column 

Content 

Description 

1 

1-2 

01 

Sequence  number 

Z 

3-4 

Blank 

Iteration  number 

3 

5-9 

TESSR 

Principal  variable  identifier 

4 

10,  11 

n,  m 

J 

4 

12 

Blank 

n,  m 

4 

13,  14 

n,  m 

J 

4 

15 

Blank 

n,  m 

4 

16,  17 

n,  m 

J 

4 

18 

Blank 

n,  m 

4 

19,  20 

n,  m 

J 

4 

21 

Blank 

n,  m 

4 

22,  23 

n,  m 

J 

5 

24-28 

Blank 

n,  m 

6 

29,  30 

n,  m 

J 

6 

6 

31 

32,  33 

n,  m 

n,  m 

J 

n,  m 

Figure  2-39.  TESSR  Card  Format 
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TESSR 


TESSR 


Field 

Column 

Content 

Description 

6 

34 

Blank 

6 

35,  36 

n,  m 

J 

n,  m 

6 

37 

Blank 

6 

38,  39 

n,  m 

J 

n,  m 

6 

40 

Blank 

6 

41,  42 

n,  m 

J 

n,  m 

7 

43-47 

Blank 

8 

48,  49 

n,  m 

J 

n,  m 

8 

50 

Blank 

8 

51,  52 

n,  m 

J 

n,  m 

8 

53 

Blank 

8 

54,  55 

n,  m 

J 

n,  m 

8 

56 

Blank 

8 

57,  58 

n,  m 

J 

n,  m 

8 

59-61 

Not  used 

9-10 

62-80 

Not  used 

For  J 

:  n  =  3,  •  • 

•  6,  m  =  1 ,  •  • 

•  5,  and  n  >  m.  Fourteen  values  of 

n,  m 


m  are  Poss^kle  and  anY  number  of  values  from  one  through  14  may  be 

used.  They  need  not  be  in  order  for  n  or  m  but  two  or  more  identical 

values  of  J  on  any  one  card  are  not  permitted, 

n,  m  7 


Figure  2-39,  TESSR  Card  Format  (Continued) 
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2.  1.  3.  7.  1  Deck  Setup  for  Potential  Model  2,  This  model  includes  the 
coefficients  of  the  zonal  harmonics  through  J^.  A  99-card  is  requir 
to  change  the  zonal  constants  J^,  J ^  and  from  those  stored  in  progr 


Figure  2-40.  Model  2  Card  Setup 


Z^NAL  MOD  2 


7.0NA  L  MOD  2 


Field 

Column 

Content 

1 

1-2 

01 

2 

3-4 

Blank 

3 

5-9 

Z0NAL 

4 

10-17 

nl,r  in  every  column 

4 

18-20 

Blank 

4 

21-23 

Not  used 

5-10 

24-80 

Not  used 

Figure  2-41.  Z0NAL  Card  Format,  Potential  Model  2 


2-89 


99  MOD 


99  MOD  Z 


8 


0000000000000000000 

11  11  17  11  14  1)  It  11  II  It  M  21  23  23  71  »  »  27  71  7 

1  11  11  11  11  11  11  1  11  1  1  1  J  1 


0  0  0  0  0  0  0  0  0  0  0  0  (  0  0  0  0  C  0  0  0  0  0  0  0  0  0  0  0  0  0  0 


X  31  r  11  X  35  3t  37  M  M  40  41  * 

111111111111 


22222222222222 
333333333333333 
i  444444444444444 

555555555555555 
4666666666666666 
7  7  7  7  7  7  7  7  7  7  7  7  7  7 
88888888888888) 


2222222222222222222 

33333333333333333 

44444444444444444 

55555555555555555 

66666666666666666 

777777777777777777 

888888888888888888 


43  44  *5  a  4  «4  «•  50  5!  52  S3  54  55  H  5?  *  51  «  II 
1111111111111111111 

2  2  2  2  21  2  2  2  2  2  2  2  2  2  2  2  2  2  2 

43333333333333333333 

4(4  444444444444444444 

5|5  5  5  5  5|  5  5  5  5  5  5  5  5  5  5  5  5  5  5 
66666666666666666666 
7777777777777777777 


9  9  9  9  9  9  9  9  9  9  9  9  9 

10  H  1?  13  14  ii  It  if  II  It  a  21  22  23 


9  6  9  9 
u  n  a  27  a 


9999999999999 

2tM3l&»)4)S3l3?X3l4l4t 


9  9  9  9  S 

44  45  *4  4 


699999999999999 

7  4l4«SaSlUUS4»MSTMMMI 


F  ield 

Column 

Content 

i 

1  -z 

99 

z 

3-4 

Blank 

3 

5-9 

00159 

4 

10-Z3 

108Z.  48  AAAA 

-06 

5 

Z4-Z8 

00160 

6 

Z9-4Z 

-Z.56Z  AAAAA 

-06 

7 

43-47 

00161 

8 

48-61 

-1.84  AAAAAA 

-06 

9  and  1  0 

6Z-80 

Blank 

Figure  Z-4Z.  99-Card  Format,  Potential  Model  Z 
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2 .  1 . 3 . 7 . 2  Deck  Setup  for  Potential  Model  3  .  Model  3  consists  of  the 
coefficients  of  the  zonal  harmonics  through  J^.  (J^,  J^,  are  the  alter¬ 

nate  values)  and  the  second  coefficient  of  the  sectorial  harmonic, 


/ 


SECTR 


/ 


99 


Z0NAL  Card  format  shown  in  Figure  2-41. 
99  Card  format  shown  in  Figure  2-42. 


Figure  2-43.  Model  3  Card  Set-up 
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SEC/l  E  MOD  i 


SECTR  MOD 


Field 

Column 

Content 

1 

1-2 

01 

2 

3-4 

Blank 

3 

5-9 

SECTR 

4 

10 

m  }  it 

4 

11-14 

Blank 

4 

15-23 

Not  used 

5-10 

24-80 

Not  used 

Figure  2-44.  SECTR  Card  Format,  Potential  Model  3 
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2.  1.  3.  7.  3  Deck  Setup  for  Potential  Model  4.  Model  4  includes  the 
coefficients  of  the  zonal  harmonics  through  J^,  the  second,  third,  and 
fourth  coefficients  of  the  sectorial  harmonics,  and  the  coefficients  of  the 
tesseral  harmonics  through  ^ .  Alternate  values  of  J^,  J  ^ are 
required  as  in  Model  2. 


r 


r 


TESSR 


SECTR 


99 


Z(^NAL  Card  format  shown  in  Figure  2-41. 
99  Card  format  shown  in  Figure  2-42. 


Figure  2-45.  Model  4  Card  Setup 


SECTR  MOD  4 


SECTR  MOD  4 


Field 

Column 

Content 

1 

1-2 

01 

2 

3-4 

Blank 

3 

5-9 

SECTR 

4 

10-  12 

n  1 n  in  ev 
column 

4 

13-14 

Blank 

4 

15-23 

Not  used 

5-  10 

24-80 

Not  used 

Figure  2-46.  SECTR  Card  Format,  Potential  Model  4 
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TESSR  MOD  4 


TESSR  MOD  4 


F  ield 

Column 

Content 

i 

1-2 

01 

2 

3-4 

Blank 

3 

5-9 

TESSR 

4 

10-11 

31 

4 

12 

Blank 

4 

13-14 

32 

4 

15 

Blank 

4 

16-17 

41 

4 

18 

Blank 

4 

19-20 

42 

4 

21 

Blank 

4 

22-23 

43 

5 

24-28 

Blank 

6 

29-42 

Blank 

7 

43-47 

Blank 

8 

48-58 

Blank 

8 

59-61 

Not  used 

9-10 

62-80 

Not  used 

Figure  2-47.  TESSR  Card  Format,  Potential  Model  4 
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2.  1.  3.  8  Radiation  Pressure 


ESPOD  accounts  for  radiation  pressure  as  a  function  of  spacecraft 
mass  and  effective  area.  The  effective  area  is  a  complex  function  of  the 
surface  reflection  and  the  type  of  reflection.  See  ESD- TDR -  64- 395, 
"Mathematical  and  Subroutine  Description, 11  Section  5.  11,  for  a  discussion 
of  the  radiation  pressure  model  in  ESPOD.  However,  the  effective  area 
is  at  least  equal  to  the  spacecraft  area  and  is  typically  10  percent  greater. 
The  effect  of  satellite  eclipsing  may  be  accounted  for  by  multiplying  the 
effective  area  by  (1  -  f),  where  f  is  the  fraction  representing  the  time  the 
satellite  is  in  earth  shadow. 

For  the  inclusion  of  radiation  pressure  effects,  a  RADPR  card  speci¬ 
fying  the  mass  and  area  (modified,  if  needed)  must  be  included  in  the  pre¬ 
liminary  data.  Its  format  is  shown  in  Figure  2-48. 
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RADPR 


RADPR 


ODOO  000000000000000; 

10  11  t2  I)  14  IS  If  17  If  It  n  21  22  23$ 
1111111111111 


22222222222222^ 
3333  3|3  333333333333  3] 
44  4  444444444444444  4 \ 
5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  ] 
%  S  S  6  S  6  6  6  6  6  6  6  6  6  6  6  6  6 

7777777777777/3 
8888888888888 3 3 
9  3$ffi9  999  9  9999999999999  9 E 

7  I  f  7  I  C  JO  1MJ  tj  U  IS  If  17  II  If  20  21  22  ’ll 


00000000000000 
30  31  32  JJ  >4  33  *  J7  »  3t  40  41  K\ 

11  I  M  11  11  1  M  1 
22222222222222 
33333333333333 
44444444444444 
5  5  5  5  5  5  5  5  5  5  5  5  5  51 
666666666SEE66 

7  7  7  7  7  T  7  7  7  7  7  7  7  7 

8  8  8  8  8  8  8  8  8  8  8  8  8  8 
9999999999999 

m)J32l334353»3?3»JJ«0  41O 


.  .  .  .  .  .  .  ir  ^ 


10 


F  ield 

Column 

Content 

Description 

i 

1-2 

01 

Sequence  number 

2 

3-4 

Blank 

Iteration  number 

3 

5-9 

RADPR 

Principal  variable  identifier 

4 

10-23 

M2 

Effective  area  of  the  space¬ 
craft  in  square  meters 

5 

24-28 

Blank 

6 

29-42 

Kg 

Mass  of  the  spacecraft  in 
kilograms 

7-10 

43-80 

Not  used 

Figure  2-48.  RADPR  Card  Format 
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2 .  1 .  3 .  9  Matrix  Inputs 


Two  matrix  inputs  are  used  by  ESPOD.  SMAT  conditions  the  A  A 
matrix  before  a  one-pass  differential  correction.  UPMAT  provides  an  a 
priori  covariance  matrix  for  an  ESPCjDEPH  run  only. 

2.  1.3.9.  1  SMAT.  For  a  standard  run  by  ESPOD,  all  observations  will  be 
present  in  the  input.  However,  ESPOD  can  accept  an  a  priori  A^A-matrix 
(input  flag  on  column  42  of  JDC  car  d)  which  brings  into  the  differential  cor¬ 
rection  summarized  data  from  observations  already  processed  on  previous 
runs  of  ESPOD.  The  a  priori  A^A-matrix  will  be  valid  only  if  the  initial 
conditions  for  the  new  differential  correction  are  the  same  set  which  gen¬ 
erated  the  a  priori  A^A-matrix.  Note  that  the  new  differential  correction 
will  cause  the  initial  conditions  to  vary  from  iteration  to  iteration.  Hence, 
the  a  priori  A^A-matrix  is  valid  for  only  one  iteration. 

T 

The  a  priori  A  A-matrix  (S-matrix)  is  input  as  an  upper  triangular 
matrix  on  SMAT  (S-matrix)  cards.  The  total  number  of  variables  input 
is  N(N  +  l)/2,  the  maximum  order  being  27.  They  are  ordered  by  row 
and  columns  within  row,  starting  with  the  first  element  input.  The  follow¬ 
ing  example  indicates  the  placement  of  elements  for  a  5  x  5  S-matrix: 


Card 

1: 

sir 

S12’ 

S1  3’ 

S14 

Card 

Z: 

S15’ 

S22’ 

S23’ 

S24 

Card 

3: 

S2  5’ 

S33’ 

S34< 

S35 

Card 

4: 

S44> 

s45> 

S55 
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SM A  I  Nos.  I.  2.  3 


SMA'I  Nos.  1,2,5 


F  ield 

Column 

Content 

Card  1 

Card  2 

Card  3 

i 

1-2 

01 

02 

03 

2 

3-4 

Blank 

Blank 

Blank 

3 

5-9 

SMAT 

SMAT 

SMAT 

4 

10-23 

S11 

S1  5 

S25 

5 

24-28 

Blank 

Blank 

Blank 

6 

29-42 

S12 

S22 

S33 

7 

43-47 

Blank 

Blank 

Blank 

8 

48-61 

S 1  3 

S23 

S34 

9 

62  -  66 

Blank 

Blank 

Blank 

10 

67-80 

S 14 

S24 

S35 

Figure  2-49.  SMAT  Card  Format,  Cards  1.  2,  and  3 
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SMA'l  NO  4 


SMA'I  NC)  4 


Field 

Column 

Content 

1 

1-2 

04 

2 

3-4 

Blank 

3 

5-9 

SMAT 

4 

10-23 

44 

5 

24-28 

Blank 

6 

29-42 

S  _ 

45 

7 

43-47 

Blank 

8 

48-61 

S55 

9-10 

62-80 

Blank 

Figure  2-50.  SMAT  Card  Format,  Card  4 
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2.  1.  3.  9.  2  UPMAT.  ESP0DEPH  can  be  operated  with  or  without ESPCjDDC . 
If  it  is  preceded  by  ESP0DDC,  the  covariance  matrix  resulting  from  the 
differential  correction  is  available  for  propagating  uncertainties.  If  it  is 
not  preceded  by  ESP0DDC,  UPMAT  is  available  for  inputting  an  a  priori 
covariance  matrix.  Note  that  uncertainties  in  orbit  elements  are  updated 
by  ESP0DEPH,  thus  UPMAT  is  restricted  to  being  6x6,  7x7,  or  8  x  8 
matrix  depending  on  whether  a  6  x  6,  drag  and  drag  variation  are  or  are 
not  included.  The  covariance  matrix  is  symmetric  and  only  the  lower 
triangular  portion  is  input. 

As  an  example,  a  6  x  6  matrix  is  illustrated  showing  the  order  of 
assigning  matrix  elements  to  cards,  followed  by  the  format  of  the  UPMAT 
cards . 


U11 

U21  U22 

U31  U32 

U33 

U41  U42 

U43 

U44 

U51  U52 

U53 

U54 

U55 

U6 1  U62 

U63 

U64 

U65 

U66 

UPMAT,  card 

input  order: 

Card  1: 

uir 

U2 1  ’ 

U22  ’ 

U31 

Card  2: 

U32’ 

U33' 

U4 1  ’ 

U42 

Card  3: 

U43’ 

U44’ 

U  5 1  ’ 

U 

52 

Card  4: 

U  53’ 

U54’ 

U55’ 

U6 1 

Card  5: 

U62’ 

U63’ 

O 

4^ 

U65 

Card  6: 

U66 

2-104 


UPMAT  Nos.  1.  2,  3,  4.  5 


UPMAT  Nos.  1.  2.  3.  4.  5 


/ 


0000000000000 

II  11  17  13  14  IS  l|  17  II  II  »  21  n  n 

2222222222222 
3333333333333 
44444444444444] 
555555555555555 
66666666666666] 
7777777777777  7§ 
88988888888888) 
99999999999999! 

10  11  17  13  M  15  II  17  II  II  »  21  22  23  5 


0  0  0  0  0  0  0  0  0  0  0  0  0  0! 

21  JO  31  32]J34]SJt37JI)9  40  41 
11111111111111 

22222222222222 

33333333333333 

44444444444444 

55555555555555 

66666866666666 

7777^777777777 


19  9  9  9  9  9  9  9  9  9  9  9  9  9. 

|7t  »  31  J3  33  M  W  %  V  3»  3S  40  41  4! 


10 

>0  0  0  0  0  0  0  0  0  0  0  0  0  0 

17  IS  M  7*  71  72  73  74  K  7*  77  7|  71  M 

11111111111111 
22222222222222 
33333333333333 
44444444444444 
55555555555555 
:6  6  6  6  6  6  6  6  6  6  6  6  6  6 
11111111111111 


99999999999999 

||7  60  II  70  71  72  *3  74  7S  71  77  7|  79  00 


Field 

Column 

Content 

Card  1 

Card  2 

Card  3 

Card  4 

Card  5 

1 

1-2 

01 

02 

03 

04 

05 

2 

3-4 

Blank 

Blank 

Blank 

Blank 

Blank 

3 

5-9 

UPMAT 

UPMAT 

UPMAT 

UPMAT 

UPMAT 

4 

10-23 

U11 

U32 

U43 

U53 

U62 

5 

24-28 

Blank 

Blank 

Blank 

Blank 

Blank 

6 

29-42 

U21 

U33 

U44 

U54 

U63 

7 

43-47 

Blank 

Blank 

Blank 

Blank 

Blank 

8 

48-61 

U22 

U41 

U51 

U55 

U64 

9 

62-66 

Blank 

Blank 

Blank 

Blank 

Blank 

10 

67-80 

U31 

U42 

U52 

U6 1 

U65 

Figure  2-51.  UPMAT  Card  Format,  Cards  1,  2,  3,  4,  and  5 
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UPMAT  NO  6 


UPMAT  NO  6 


0000000000000 

10  11  II  13  14  13  II  17  II  II  X  II  77  23 

1 1 1 1 1 1 1 1 1 1 1 1t  1 
2222222222222 
33333333333333 
44444444444444 
5555555555555 
66666666666666 
1111111111111 
88888888888888 
99999999999999 

10  11  II  13  14  13  11  17  II  II  20  II  22  23 


Field 

Column 

Content 

1 

1-2 

06 

2 

3-4 

Blank 

3 

5-9 

UPMAT 

4 

10-23 

U66 

5-10 

24-80 

Blank 

Figure  2-52.  UPMAT  Card  Format,  Card  6 
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2.  1.  3.  10  ESP0DEPH  Options 


ESP(^DEPH  integrates  the  orbit  from  epoch  to  some  specified  time 
or  times,  computes  the  position  and  velocity,  and  can  compute  the  uncer¬ 
tainty  in  position  and  velocity.  Specification  of  update  times  at  which 
the  satellite  position  will  be  calculated  can  be  input  on  three  types  of  cards 
Only  one  type  of  input  card  must  be  used.  They  are  discussed  separately 
below. 

2.  1.  3.  10.  1  DELTT.  Update  times  may  be  specified  at  which  the  satellite 
position  will  be  calculated  and  printed  by  the  DELTT  (AT  and  T)  table. 
Times  are  measured  from  epoch  and  integration  proceeds  monotonically 
away  from  epoch  in  either  a  positive  or  a  negative  direction.  When  the 
direction  changes,  integration  is  reinitiated  at  epoch. 

AT  specifies  the  step  size  between  successive  calculated  positions; 

T  specifies  the  last  point  at  which  the  calculation  is  desired.  Up  to  eight 
AT,  T  pairs  (four  cards)  may  be  specified. 

As  an  example,  the  positions  to  be  calculated  may  be: 

(1)  80  minutes  before  epoch  and  at  every  2  minutes  for 
20  minutes 

(2)  170  minutes  before  epoch  and  at  every  2  minutes  for 
20  minutes 

(3)  10  minutes  after  epoch  and  at  every  2  minutes  for 
20  tninutes 

(4)  100  minutes  after  epoch  and  at  every  2  minutes  for 
20  minutes 

Time  scale: 


4  “I* 

420r 

0 

420  r 

420r 

-170 

-80 

i  +  10 

+  100 

TIME-MINUTES 

EPOCH 
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To  specify  this  output,  the  following  DELTT  cards  are  input: 


Card  1  Card  2 


AT1  = 

-60 

AT  ^  =  -70 

Tl  = 

-60 

T3  =  -150 

at2  = 

-2 

AT.  =  -2 

4 

T2 

-80 

T4  =  -  170 

Card  3 

Card  4 

at5  = 

10 

AT7  =  70 

T5 

10 

T?  =  100 

AT6 

2 

AT8=  2 

T6 

30 

T8=  180 

2-109 


DELTT 


DELTT 


010000000000000000 


i  3  I  III 

Mill 


I*  11  12  13  14  13  13  17  II  IS  20  21  22 

M  M  1  1  1  1  1  1  1  1  1  1 


222222222222222222 

333333333333333333 

4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  41: 
555555555555555555 

5  6  6  G  6  6  6  6  6  6  6  6  6  6  6  6  6  6  £ 
777777777777777777  7; 
I  8  8  8  I  8  8  8  8  8  t  8  8  1  1  8  8  8 

(939999999999999999 

3  1  J  I  *  jo  11  12  13  14  13  11  17  II  13  20  21  22 


10 


00000000000000 

57  M  M  71  71  77  71  74  73  71  77  71  7|  M 

11  I  I  I  I  I  I  I  I  I  1  I  1 
22222222222222 
33333333333333 
44444444444444 
55555555555555 
66666666666666 
11111111111111 
88888888888888 
99999999999999 

C7  H  19  70  71  7?  >3  74  73  N  77  70  71  « 


Field  _ Content _  Description 


Card  1 

Card  2 

Card  3 

Card  4 

1 

01 

02 

03 

04 

Sequence  number 

2 

Blank 

Blank 

Blank 

Blank 

Iteration  number 

3 

DELTT 

DELTT 

DELTT 

DELTT 

Principal  variable 
identifie  r 

4 

AT. 

AT 

at5 

AT7 

5 

Blank 

Blank 

Blank 

Blank 

6 

T1 

T 

3 

T5 

T 

7 

7 

Blank 

Blank 

Blank 

Blank 

8 

at2 

AT4 

AT6 

AT8 

9 

Blank 

Blank 

Blank 

Blank 

10 

T 

2 

T4 

T6 

T 

1 8 

Figure  2-53.  DELTT  Card  Format,  Cards  1,  2,  3  and  4 
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2.1.3.10.2  DAC  Cards.  The  presence  or  absence  of  these  cards  is 
indicated  in  Column  52  of  the  JDC  card  (see  Section  2.1.2). 

DAC  cards  are  an  input  to  ESPOD  to  provide  an  interface  with 
special  predicting  and  comparison  requirements.  DAC  cards  typically 
come  in  a  set  of  three: 

1)  Header  card 

2)  11  08"  card  with  test  Cartesian  coordinates  and  days 

3)  ”  12”  card  with  fractions  of  days. 

The  header  card,  if  input,  is  ignored.  The  predict  time  is  taken 
from  the  08  and  12  cards  (second  and  third  cards  of  the  three  card  set). 
Cartesian  position  is  used  as  a  comparison  point  for  ESPOD1  s  predicted 
value . 

ESPOD  will  accept  up  to  a  maximum  of  four  sets  of  DAC  cards. 
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DAC,  No.  1 


DAC,  No.  1 


/ 


4 


7 


Xp  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0 

I  n  U  I  M  )  n  II  11 11  H  15  II  1 7  14  II  2»  21 
111111111111111111111 


oooooooooooooooooooooooooooooooooeoooooooooooo 

23  24  25  X  77  2l2«J0  3irU34&JI37JIJ»4l  41  42  43  44  45  4S  47  4|4S30S1S7S3  54S5SIS7SI3«SQItt2UI4lSttU 


0  0  0  0  0  0  Q00IIM 

61  70  71  72  73  74  71  *  7  7  7«  ?f  *0 

1  1  1  1  1  1  1  1  1  1  I  1  I 


3  3  3 


4  4 


2  2  2  2  2|2 

3  3  3  3 

4  4  4  4 
555555 
668666 
777777 
8  8  8  8 


8  8  8 


222222222222222 

33333333333333 

444444444444444 

555555555555555 

666666666666666 

7777777777777777 

88888888888888 


22222222222222222222222222222222222222222272 

33333333333333333333333333333333333333333333 

44444444444444444444444444444444444444444444 

55555555555555555555555555555555555555555555 

66666666666666666666666666666666666666666666 

7777777777777777777777777777777777777777777 

88888888888888888888888888888888888888888888 


2  2  2 


222222222222 
331333333333 
44444444444444 
55555555585555 
668666666666  » 
7777777777777 
888888888888 


3  3  3 

4 

5 

6  6  6 
7  7 


8  88 


9  9  9  9  9 

1  2  3  4  5 


9999999999999999 


S  7 


I  I  .0  11  12  13  14  IS  It  17  11  13  20  21 


9  9  9  3  9  9  9  9  9  9  9  9  9  9  9  8  9  1  9  9  9  9  9  8  9  8  9  9  9  )  9  9  9  9  9  8  9  9  9  9  9  9  9  9 

23  24  252l?7  3l2tM  31  32  33  >433S  37  3B3l«  41  42  43  4<4S4t  47  4MtllSlUS3M9»JI  57  5l5««0«l£2l3HI5 


999 


999999999999 

I  70  71  7?  73  74  75  71  77  71  TS  « 


Field 

Column 

Content 

De  sc  r  lpt  ion 

1 

1-5 

Same  as  Input  Header  Card 

Z 

6-7 

ss 

Source  code 

3 

8-21 

Same  as  Input  Header  Card 

4 

22 

N 

ID  number 

5 

23-64 

Same  as  Input  Header  Card 

6 

65-67 

NNN 

Sequence  number 

7 

68-80 

t  "  “ 

Filler 

This  first  card  of  the  three  card  set  (DAC)  is  not  necessary  for  input. 
However,  ESPOD  will  read  it  and  ignore  the  contents. 


Figure  2-54.  DAC  Card  Format,  Card  1  of  3 
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DAC,  No.  2 


DAC,  No.  2 


T" 

213 


rr 

3  4 


6  7 


8 


910 


11 


14 


o  o  o 

1  I 
1  1 


0000000000000 

)  4  5  I  71  •  II  It  12  I)  14  is  n 
11111111111111 


2  2  2  2  2  2  2  2  2  2  2  2  2  2  212 

3  3  3  3  3  3  3  3  3  3  3  3  3  3  3|3 

4  4  4  4  4  4  4  4  4  4  4  4  4  4  4 14 
555555555555555 

6  6  6  6  6  6  6  6  6  6  6  6  6  6  616 
7777777777777777 
8888888888888888 


9  9 
1  2 


99999999999999 

3  4  5  I  70  I  10  M  12  1)  14  IS  II 


0  0 

17  II 

1  1 

2  2 

3  3 

4  4 

5  5 

6  6 

7  7 

8  8 


0000000000000 

ll  21  II  22  21  H  n  21  27  2»  21  30  31 
1111111111111 

2222222222222 

3333333333333 

4444444444444 

5555555555555 

6666666666666 

7777777777777 

8888888888888 


999999999999999 


9  9 

1*  20  21  22  23  24  25  21  27  21  21  30  M 


0000000000000000 

U  31  37  31  31  40  41  42  43  44  45  41  47 
1111111111111 


32 

1  1  1 

222222222222222  2|2  22 

3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3|333 

4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4(4  4  4 

555555555555555555 

6666666666666666666 

7777777777777777777 

8888888888888888888 


9999999999999 

35  X  37  31  31  40  41  42  43  44  45  41  47 


99 


0000000000000  010  00 

I152U5455X  57  5OXMI1I2I3MX0I7 

11111111111111111 


22222222222222 

33333333333333 

44444444444444 

555555555555555 

66666666666666 

77777777777777 

88888888888888 


999999999999999 

51  52  U  54  55  X  57  54  55  60  Cl  I?  43 


0000000000000 

I  71  71  72  71  74  75  71  77  71  71  II 

11  11  1  11  11  11  1  1 


2222222222222222 

3333333333333333 

4444444444444444 

5555555555555555 

6666666666666666 

77777/7777777777 

8888888888888888 


9  9  9  9 

a  u  i7 


999999999999 

U  M  70  7>  7?  73  74  75  71  77  7|  71  » 


Field 

Column 

Content 

Desc  ription 

1 

1-2 

08 

Identifying  code  for  X 

2 

3-15 

±mm - mm±ee 

X  -  Earth  radii 

3 

16 

1 

Separator 

4 

17-18 

09 

Identifying  code  for  Y 

5 

19-3  1 

±mm - mm±e  e 

Y  -  Earth  radii 

6 

32 

t 

Separator 

7 

33-34 

10 

Identifying  code  for  Z 

8 

35-47 

±mm - mm±ee 

Z  -  Earth  radii 

9 

48 

1 

Separator 

10 

49-50 

11 

Identifying  code  for  time  - 
whole  days 

11 

51-63 

±mm - mm±ee 

Time  -  whole  days  from 
January  0,  Jan  1=1 

12 

64 

9 

Separator 

13 

65-67 

NNN 

Sequence  number 

14 

68-80 

f  ~  ~  ~  > 

F  iller 

Figure  2- 

55.  DAC  Card  Format,  Card  2  of  3 
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DAC,  No.  3 


DAC,  No.  3 


/ 


1 

2. 

3 

4 

5 

00 

1  2 

1 1 

2  2 

3  3 

4  4 

5  5 

6  6 

7  7 

8  8 

9  9 

1  2 

0000000000000 

1  4  S  1  7  1  1  11  It  11  IJ  14  11 

1 1 11 1 1 1 11  1 1 1 1 

2222222222222 

3333333333333 

4444444444444 

5555555555555 

6666666666666 

7777777777777 

8888888888888 

9999999999999 

)  «  i  u  n  w  ii  ii  i]  >4  it 

0000000000000000000000000000000000000000000000000 

11  1MI  11  »  71  72  a  N  21  21  71  21  2*  30  11  12  U  M  M  ft  37  31  W  40  41  4?  43  44  45  H  4)  4l  «  »  51  51  U  54  55  M  S7  U  SO  M  11  *7  «  M 

I  1  1  1  1  1  1  1  1  1  1  1  1  11  1  1  1  1  1  1  1  1  1  1  II  1  1  1  1  1  1  1  1  1  1  11  1  1  1  1  1  1  1  1  1  1 

2222222222222222222222222222222222222222222222222 

3333333333333333333333333333333333333333333333333 

4444444444444444444444444444444444444444444444444 

5555555555555555555555555555555555555555555555555 

6666666666666666666666666666666666666666666666666 

lnmmmmmmmmmnmiimnmmi 

8888888888888888888888888888888888888888888888888 

9999999999999999999999999999999999999999999999999 

II  IT  11  if  .*0  21  22  2J  24  2b  21  IT  2#  2*  30  J1  12  13  M  M  X  1MI  20  40  41  42  41  44  «S  47  4|  4J  51  U  U  S4  SS  %  57  SI  59  SI  11  17  »  M 

0  0  0 
MM  11 

1  1  1 

2  2  2 

3  3  3 

4  4  4 

5  5  5 

6  6  6 

7  7  7 

8  1  8 

999 

M  M  11 

0000000000000 
m  m  w  n  n  u  it  n  n  n  n  m 

1  1  1  1  1  1  1  1  1  1  1  1  1 

2222222222222 

3333333333333 

4444444444444 

5555555555555 

6666666666666 

7777777777777 

8  8  S  8  8  8  8  8  8  8  8  8  8 

9999999999999 

M  Vi  JO  71  7 1  >1  74  75  71  7  7  71  7)  00 

Field 

Column 

Conte  nt 

De  script  ion 

1 

1-2 

12 

Code  for  time  -  fraction 
of  day 

2 

3-15 

±mm - mm±ee 

Time  -  fraction  of  day 

3 

1 6  -  64 

»  """» 

F  iller 

4 

65-67 

NNN 

Sequence  number 

5 

68-80 

t  t 

Filler 

Note  that  the  field  designation  of  the  DAC  cards  departs  from  the  ESPOD 
preliminary  data  format.  For  DAC  cards,  any  fields  not  containing  input 
parameters  are  filled  with  commas  as  filler  characters.  This  applies 
to  field  14  of  card  2,  and  fields  3  and  5  of  card  3. 


Figure  2-56.  DAC  Card  Format,  Card  3  of  3 
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2.1.3.10.3  PRDCT.  The  PRDCT  card  specifies  from  one  to  four  separate 
output  times  at  which  to  predict.  Times  are  specified  in  days  and  frac¬ 
tions  of  days,  given  in  a  single  variable  field.  Only  one  predict  card  is 
permitted  (see  Figure  2-57). 
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PRDCT 


PRDCT 


o  o  0  0  0  0  0  0  0  0  0  0  OB 

10  11  12  13  14  13  II  IT  II  II  2«  21  22  2 

11111111111111 

22222222222222J 
3  3  3  3  3  3  3  3  3  3  3  3  3  ) 
44444444444444 
55555555555555) 
66666666666666 
77777777777777 
88888888888888 
99999999399999 

0  11  12  13  14  15  U  U  II  II  20  21  27  23 


V\\\V 


6 

00000000000000 

21  30  11  M  33  34  35  31  3J  31  31  40  41  42 

11111111111111 

22222222222222 

33333333333333 

44444444444444 

55555555555555 

66666666666666 

77777777777777 

88888888888888 

99999999999999 

21  30  31  12  3J  M  35  %  37  3*  39  4u  4’  42 


41  50  51  52  W  54  55  56  57  51  5«  M  II 


vnwVv 


10 

0000000000000 

IT  M  N  ft  71  77  TJ  ?4  71  H  77  T|  T1  M 

1111111111111 

2222222222222 

3333333333333 

4444444444444 

55555555555555 

6666666666666 

77777777777777 

88888888888888 

99999999999999 

V  n  M  70  71  72  73  74  75  ?l  V  ’I  7*  10 


Field 

Column 

Content 

1 

1-2 

01 

2 

3-4 

Blank 

3 

5-9 

PRDCT 

4 

10-23 

ll 

5 

24-28 

Blank 

6 

29-42 

l2 

7 

43-47 

Blank 

8 

48-61 

S 

9 

62-66 

Blank 

10 

67-80 

l4 

Description 
Sequence  number 
Iteration  number 
Principal  variable  identifier 
1st  update  time 

Znd  update  time 

3rd  update  time 

4th  update  time 


Figure  2-57.  PRDCT  Card  Format 
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2.1.3.11  Constant  Change  (99  Card) 


Most  of  the  constants  which  are  used  in  the  ESPOD  Program  can  be 
changed  with  a  99  card  (see  Figure  2-58).  A  list  of  the  ESPOD  constants, 
their  numerical  values,  and  relative  core  location  is  given  in  Section  5.0. 
A  more  detailed  list  is  given  in  Section  3.5.2  of  ESD- TDR- 64- 39 5  ''ESPOD 
Mathematical  and  Subroutine  Description."  Full  details  are  in  the  CODE 
EDIT  listing. 
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99 


99 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0  0 

1  1 

\  1 

2  2 

3  3 

44 

55 

1  1 

.. 

9  9 

1  I 

l 

Off  000 

MMI 

HIM 

22222 

3  3  3  3  3 

4  4  4  4  4 

5  5  5  5  5 

0  5  6  0  0 

7  7  7  7  7 

I  8  0  8  a 

9  3  9  19 

3  17  11 

0  0  0  0  0  0  0  0  0  Offff ff  o 

'0  ll  11  13  14  It  II  11  11  >1  Mil  ttfl 
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22222222222222 

33333333333333 

44444444444444 

55555555555555 

0  0  0  0  0  6  6  6  6  6  6  6  0  6 

77777777777777 

88018008008168 

99999999999999 

0  11  12  17  14  It  11  IT  11  H  XH7M3 

0  0  0  0  0 

14  7i  »  27  H 

mi  i 

2  2  2  2  2 

3  3  3  3  3 

4  4  4  4  4 

5  5  5  5  5 

6  6  0  5  5 

7  7  7  7  7 

0  0  1  0  8i 

9  9  9  9  9 

w  a  m  it  » 

00000000000000 

39  3D  31  JU  n  M  JS  *  J7  M -9  <0  41  4| 

t  m  i  m  n  n  n 

22222222222222 

33333333333333 

44444444444444 

55555555555555 

66666666666666 

77777777777777 

00808000080100 

99999999999993 

n  S  )l  U  H  U  U  )f  M  J0  4HI  41 

0  0  0  0  0 
41  u  «  4|  41 

nm 

2  2  2  2  2 

3  3  3  3  3 

4  4  4  4  4 

5  5  5  5  5 

6  6  6  6  6 

7  7  7  7  7 

8  8  0  0  0 

9  9  9  9  9 

U  44  U  «  *T 

0  0  0  0  0  0  0  0  0  0  0  0  0  0 

*»  «l  SI  51  57  5]  W  »  M  51  M  M  M  11 

It  t  1  It  t  1  11  1  11  1 

22222222222222 

3  3  3  3  3  3  3  3  3  3  3  3  3  3 

44444444444444 

55555555555355 

G  G  6  6  G  6  6  E  6  6  6  6  E  6 

77777777777777 

G  8  8  8  8  8  8  8  8  8  8  8  8  8 

99999999999999 

a  «  SO  SI  H  53  M  »  »  SI  51  M  K  SI 

0  0  0  0  0 
I2UWHN 

Mill 

2  2  2  2  2 

3  3  3  3  3 

4  4  4  4  4 

5  5  5  5  5 

6  6  6  6  6 

7  7  7  7  7 

8  8  8  8  8 

9  9  9  9  9 

1?  «  m  a  u 

00000000000000 

vumitminunitnnnu 

11111111111111 

22227222222222 

33333333333333 

44444444444444 

55555555555555 

66666666666666 

777/7777777777 

8  8  8  8  8  8  8  8  8  8  8  8  8  t 

99999999999999 

kl  Cl  19  70  71  17  73  74  75  7t  7  7  7|  71  M 

Field 

Column 

Content 

Description 

1 

1-2 

99 

Card  identifier 

2 

3-4 

Blank 

Iteration  number 

3 

5-9 

XXXXX 

Core  location  of  new  constant 

4 

10-23 

New  constant 

5 

24-28 

XXXXX 

Core  location 

6 

29-42 

New  constant 

7 

43-47 

XXXXX 

Core  location 

8 

48-61 

New  constant 

9 

62-66 

XXXXX 

Core  location 

10 

67-80 

New  constant 

Figure  2-58.  99  Card  Format 
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2.1.3.12  End  of  Preliminary  Data 


An  ENDPR  (End  of  Preliminary  Data)  card  is  mandatory  whether 
or  not  any  preliminary  cards  are  present.  If  any  are  present,  it  must 
follow  directly  behind  the  last  one  of  them.  An  ENDPR  card  is  shown 
in  F igure  2-  59  . 
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FN'D  PR 


END  PR 


Field 

Column 

Content 

1 

1-4 

Not  used 

2 

5-9 

ENDPR 

3 

10-80 

Not  used 

Figure  2-59.  ENDPR  Card  Format 
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2.  1.4  Observation  Data 


2.  1.4.  1  Observation  Cards 

ESPOD  accepts  the  standard  SPADATS  observation  cards,  as  given 
in  Figure  2-60.  If  more  than  7  04  observations,  the  maximum  number,  are 
to  be  processed  on  a  single  run,  they  must  be  presorted  in  chronological 
sequence  before  ’’read  in.”  If  more  than  7  04  unsorted  observations  are 
input,  the  run  will  exit  on  an  input  error. 

When  more  than  7  04  chronologically  sorted  observations  are  read 
in  on  cards,  and  initial  conditions  are  specified  on  cards,  the  epoch  must 
be  at  or  later  than  the  time  of  the  last  observation. 

Figure  2-60  shows  the  format  of  the  SPADATS  observations  cards. 
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OB  VS 
CARD 


OBVS 

CARD 


TT“ 


•0*0000000 


1  2  I 
1 1 1 


2222222222 

3333333333 


7  7  7 


6666S6666 


888888888 


999999999 


4  I 


17  1* 
1111 


4444444444 

5555555555 


7  7  7  7  7  7 


17  11 


00000000000000 


n  i}  is  u  is|ii  u  it  ii  n  2i  a  a  »?7?i»mhi:m)4»34ji 
1111111111111 


11111 
2  2  2  2  2 


33333333333333 


4  4  4  4  4 

5  5  5  5  5 
6  6  6  6  6  6 
777777 
888888 


111111111 

222222222 


444444444 

555555555 

666666666 

777777777 

888888888 


999999999999999 


0000000000000 


1111111 

2  2  2  2  2  2K2  2  2  2  2  2  2  2  2  2  2  2  2  2|2  2  2  2  2  2  2  2  2 

3  3  3  3  3  3|3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3 

4444444 


4444444444444 

5555555555555 


66666666666666666666 
7777777 


7777777777777 

8888888888888 


99999999999999999999 


II  II  111)  14  13  11  17  II  II  2t  21  a  23  M  »  »  27  21  21  »  >1  8  B  M  31  »  37  ■  *  40  41  4MJ  44  4S  «  47  «  «  M  3MM)  M 


0  0  0  0  0  0  0  0  0  0  0  0  0  0  I  0  0  0  0  0  0  0 

31  31  41  41  42  4)  4^4  41  4}  41  41  M  II  8  If  >4*  *  51*  * 
111111111 


8888888 


10 


444444444 


5555555555555555 


6666666666666 


777777777 

888888888 

999999999 


33333333333333333 


1  1  1 
2  2  2j 


4  4  4  4 

5  5  5  5 


6  6 

7  7 

8  8 
9  9 

»  54  5*54  54 


7  7  7 

8  8  8 


9  9  9  9 


1  3 


4  4 


0  0  0  0 


NIOGIC 


1 1 


222222 


4  4 


6  6 


8  8 


1  1 


6  6 


7  7  7  7, 


8  8 


9  9  9  9 

(4  II  13  U  M 


44444444 


555555555555 


Field 

1 


Z 

3 


Column 

Des  cr  iption 

1-3 

Satellite  number.  Column  1  contains  a  minus 
sign  if  this  is  a  classified  observation;  column  1 
cannot  be  double  punched. 

4-5 

Equipment  type 

6-9 

Station  number 

4  10  Accuracy  or  signal  strength  (see  Supplement  1  - 

OBVS  card,  Page  Z-1Z8). 

5  11-15  Date 


6 

7 

8 

9 


16-24  Time  (Z) 

25-3  0  Elevation/declination.  Column  2  5  can  be  over¬ 

punched  +  or  -. 

31-37  Azimuth/ right  ascension.  Column  3  1  can  be 

overpunched  +  or  -.* 

38-44  Slant  range  (KM) 


A  minus  overpunch  in  column  31  indicates  fields  7  and  8  are  declination 
and  right  ascension  respectively. 


Figure  2-60.  Observation  Card  Format 
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OBVS 

OBVS 

CARD 

CARD 

Field 

Column 

De  scription 

10 

45-53 

Range  rate  (KM/SEC)  with  implied  decimal  point 
between  columns  46  and  47 

or: 

maximum  frequency  shift  (cycles/sec  ) 
with  implied  decimal  point  between  columns 

52  and  53. 

11 

54 

i  0  or  A  indicates  range  rate 
l  in  field  10. 

(  1  indicates  maximum  frequenc 

quency  shift  in  field  10 

12 

55-57 

Brightness  at  ^ 

or: 

observation  time 

when  column  58  contains 
a  minus  punch,  then 

13 

58-59 

Maximum 

columns: 

brightness 

l  55-57  contain  radar  cross 

/  section  in  meters 

14 

60-6  1 

Minimum 

59-63  contain  the  fre¬ 

brightness 

quency,  with  decimal 
point  between  60  and  61. 

15 

62-63 

Time  interval  j 

16 

64-65 

Date  or  line  number  (NOTE:  Not  used  by  SPS) 

17 

66-69 

Message  number 

18 

70 

Equinox  (see  Supplement  2  -  OBVS  card,  Page 

2-129). 

19 

71-72 

Not  used 

20 

73-78 

Observation  number  (assigned  by  0RC0N) 

21 

79 

Switch  indicator  used  by  manual  system 

22 

80 

Card  type  (code  type  =  Any  numeric  between  0-9) 

identifies  an  Observation  Card.  0  =  unknown, 

1-9  coded  according  to  the  Association  Status  as 
determined  in  Report  Association. 


Figure  2-60.  Observation  Card  Format  (Continued) 


SUPPLEMENT  1  — OBSERVATION  CARD  (OBVS) 


Either  accuracy  or  signal  strength  may  be  indicated  in  column  10 
according  to  the  following: 

If  type  in  columns  4  and  5  is  31  or  greater,  column  10  contains 
signal  strength.  If  type  is  30  or  less,  column  10  contains  accuracy. 


Code 


Accuracy 


Signal 


0  Normal  observations,  Signal  strength  good, 

made  under  fair  conditions  reliable  measurement 


1 


2 


3 


Observations  slightly 
under  par  due  to 
outside  interference 
(e.g.,  some  clouds, 
reduced  visibility) 

Observations  poor  only 
due  to  outside 
inte  rf  e  rence 


Signal  fair 


Signal  weak, 
results  poor 


Only  estimates  possible  Signal  questionable 

(malfunction  of  instru¬ 
ment,  observation  time 
span  short) 
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I 


SUPPLEMENT  Z  —  OBSERVATION  CARD  (OBVS) 


Column  7  0  contains  year  of  Equinox  as  specified  by  following 


0  Year  to  date 

5  =  2000 

1  =  1900 

6  =  1850 

2  =  1920 

7  =  1855 

3  =  1950 

8  =  1875 

4  =  197  5 

9  =  I960 

2-129 


2.  1.4.2  End  of  Observation  Data 


Following  observation  cards,  it  is  mandatory  to  include  an  END0B 
card.  The  format  is  given  in  Figure  2-61. 
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ENDtf)B 


END0B 


F  ield 

Column 

Content 

i 

1-4 

Not  used 

2 

5-9 

END0B 

3 

10-80 

Not  used 

Figure  2-61.  END0B  Card  Format 
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2.  1.  5  Sensor  Data 


2.  1.  5.  1  Sensor  Cards 

A  sensor  card  must  be  provided  for  every  different  sensor  number 
which  appears  on  the  included  observation  cards.  The  presence  of  sensor 
data  is  flagged  in  column  31  of  the  JDC  card  (Section  2.  1.  2).  ESPOD 
assumes  the  new  SPADATS  sensor  card  format  (1  January  1964),  as 
shown  in  Figure  2-62. 

Sensor  cards  do  not  require  any  punching  to  identify  sigma  types. 
The  sigma  type  code  digit  for  each  sensor  is  stored  in  a  table  within 
ESPOD.  If  it  is  desired  on  any  one  run  to  vary  the  sigma  type  associated 
with  the  sensor,  an  STYPE  card  for  the  sensor  must  be  entered.  Any 
such  change  is  retained  until  the  next  cold  start. 

The  estimated  standard  deviations  for  each  sigma  type  are  stored 
in  the  SIGMA  table  within  ESPOD.  There  are  60  different  sets  of  sigmas 
identified  by  code  types  from  0 1  to  60.  If  it  is  desired  on  any  one  run  to 
change  the  four  values:  cr^ ,  o*^  (o*^),  °*e  in  any  set,  a  SIGMA 

card  respecifying  the  four  values  by  type  must  be  entered.  As  in  the 
STYPE  change,  the  respecified  values  will  be  retained  until  the  next  cold 
start. 

Credence  dependent  weights  are  obtained  by  multiplying  the  tabular 
sigmas  by  a  function  of  the  credence.  This  combined  with  the  proper 
initial  sigma  set  will  produce  closely  the  standard  deviations  tabulated  by 
credence  for  the  FPS-49  radars. 
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SENSOR 


SENSOR 


Field 

Column 

Description 

1 

1-4 

New  sensor  number 

2 

5-11 

4>°(+N)  —  latitude  (decimal  assumed  between  columns  7-8) 

3 

12-19 

\°(+W)  —  longitude  (decimal  assumed  between  columns  15-16) 

4 

20-25 

H(meters)  — altitude 

5 

26-32 

Datum  plane  r 

6 

33 

..•  4.-  Ju  =  unclassified 

Classification  (  .  . . .  , 

Ic  =  classified 

7 

34-36 

Not  used 

8 

37-50 

Location  (name) 

8 

51-54 

Temp  (if  the  sensor  is  a  temporary  one) 

9 

55 

Not  used 

10 

56-58 

Old  sensor  number 

1 1 

59 

Not  used 

12 

60-69 

R  emarks 

12 

70-71 

WSRN  or  Minitrack  station  number 

12 

75-78 

SAO  station  number 

13 

79 

^indicates  in  MS"  file 

14 

80 

nSn  to  indicate  sensor  card 

Figure  2-62.  Sensor  Card  Format 
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2.  1.  5.  2  End  of  Sensor  Data 


Following  sensor  cards,  it  is  mandatory  to  insert  an  ENDSN  card. 
The  format  is  shown  in  Figure  2-63. 
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ENDSN 


ENDSN 


F  ield 

Column 

Content 

i 

1-4 

Not  used 

2 

5-9 

ENDSN 

3 

10-80 

Not  used 

Figure  2-63.  ENDSN  Card  Format 
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Z.Z  DECK  SETUP 


This  section  explains  the  selection  and  input  order  of  input  cards 
for  ESPOD.  The  first  part  of  the  section  is  an  input  diagnostic  summary, 
which  is  essentially  a  tabulation  of  the  program  response  to  omitted  and 
included  input  cards.  This  is  followed  by  a  classification  of  input  cards 
by  type  that  are  mandatory,  conditional  mandatory,  and  optional.  Sample 
ESPOD  load  sheets  are  included  as  a  part  of  the  discussion  of  their  format 
and  content.  This  section  is  concluded  with  a  discussion  of  the  input  card 
sequence  (decks). 
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2.  2.  1  Input  Diagnostic  Summary 


Table  2 -III.  Program  Response  to  Omitted  and  Included  Input  Cards 

No  Card  Present  Card  Present 


Input 

in  Deck 

in  Deck 

70ASCHTP 

Run  will  not  start 
Mandatory  input 

Must  be  first  card  in  first 
deck 

Run  starts  normally 

Job  No. 

B2  system  exit 

Mandatory  B2  scheduled 
run  input 

Job  No.  printed  as  input 

RUN  ESPOD, 
DATA 

Run  will  not  start 
Mandatory  B2  system 
input 

Controlfrom  Flexowriter 
in  manual 

Reads  preliminary  data 

JDC 

Diagnostic  exit 

Mandatory  input 

Run  proceeds  as  indicated  by 
JDC  card  input 

REM 

Run  proceeds 

Remarks  printed  on  output 

ICTYP 

Diagnostic  exit,  if 

Run  proceeds 

IC0ND  cards  present 


IC0ND 


ITIME 


Automatic  reference  to 
SEAI  tape  for  elements 
and  epoch 

Diagnostic  exit  if 
ICTYP  =  1.  ,  2.  ,  or  3. 


Initial  conditions  taken  as 
input 

Epoch  indicated  is  input;  not 
required  for  ICTYP  =  4. 


7 -Card  Elements  Same  as  IC0ND 


Reads  into  program  and  is 
automatically  transformed 
to  Cartesian  coordinates 


DNREV 


No  update  will  be  made 
when  7 -card  elements 
are  input 


For  given  flags,  the  program 
will: 

F  lag  Option 

1.  Update  to  a  specific 

time 

2.  Updates  to  N  revolu¬ 

tions  past  epoch 

3.  Updates  to  a  specified 

revolution 

-1.  Updates  to  time  of  last 

observation,  when  using 
7-card  elements 
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Input 
CAT  1 


CAT2 


BNDS 


BISES 


NITER 


TSTEP 


DELET 


TMAX 


Table  2 -III.  Program  Response  to  Omitted  and  Included 
Input  Cards  (Continued) 

No  Card  Present  Card  Present 

in  Deck  in  Deck 


Assumes  a  nominal  set  of 
bounds  and  a  6  x  6  fit  if  CAT2 
also  absent.  If  CAT2  card 
present,  all  initial  conditions 
are  constrained  and  CAT2 
variables  are  solved. 

Assumes  nominal  bounds  and 
6x6  fit  if  CAT  1  absent  also. 
If  CAT1  card  present,  solves 
for  CAT1  variables  and  no 
CAT2  variables. 

If  C AT  1  or  CAT2  card  present, 
assumes  nominal  bounds. 

If  no  CAT  1  or  CAT2  cards, 
assumes  nominal  bounds  and 
a  6  x  6  fit. 

No  bias  estimates  input 

Program  runs  normally 


Assume  a  5  iteration  fit  is 
desired 


Assumes  an  initial  step  size 
of  numerical  integration  of 
one  minute 

No  observations  will  be  deleted 
from  either  the  solution  or  the 
residuals  calculation 


Observations  which  are  more 
than  ten  days  from  epoch  are 
rejected 


Solves  for  CAT1  variables 
flagged 


Solves  for  sensor  data  and  loca¬ 
tion  variables  indicated 


CAT  1  and  CAT2  variables 
bounded  as  input;  must  be  a 
bound  input  for  each  variable 
flagged  or  variable  will  not  be 
allowed  to  change. 

Initial  CAT2  variable  bias  esti¬ 
mate  input;  must  have  1:1 
correspondence  with  non-zero 
CAT2  variables  flagged. 

Allows  the  number  of  iterations 
input  with  this  card  as  a  maxi¬ 
mum.  The  program  may  stop 
sooner  if  convergence  occurs. 

Changes  the  initial  step  size  to 
specified  value 

The  number  of  observations  de¬ 
leted  will  be  those  from  a^ 
through  b^,  a^  through  b^ ,  etc. 

To  delete  a  single  observation, 
input  the  same  observation 
identification  number  for  both 
a  i ,  b i ,  or  a^,  b^  >  etc. 

Rejection  of  observations  in  days 
from  epoch  as  specified  on 
card 
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Input 

SIGMA 


STYPE 


DRAG 


APF10 


Z0NAL 

SECTR 

TESSR 

RADPR 


SMAT 


UPMAT 


Table  2 -III.  Program  Response  to  Omitted  and  Included 
Input  Cards  (Continued) 


No  Card  Present 
in  Deck 


Card  Present 
in  Deck 


Assumes  a  standard  set  which 
will  remain  unless  changed 


Assumes  standard  data  sig¬ 
mas,  average  surface  re- 
fractivity  for  each  station, 
and  a  nominal  gross  outliers 
criterion 

Assumes  drag  free  atmos¬ 
pheric  model 

Not  required  unless  a  dynamic 
atmosphere  is  flagged; 
zeroes  will  be  assumed 
when  card  is  absent. 


Will  assume  a  standard  set 
is  desired  as  a  nominal 
model  (Potential  Model  1). 

No  radiation  pressure 
perturbations 

No  a  priori  nSn  matrix  input; 
program  runs  normally. 

If  this  is  an  ESP0DEPH  run 
only  (no  ESP0DDC),  no 
a  priori  covariance  matrix 
update  is  performed;  pro¬ 
gram  runs  normally.  An 
update  can  be  made  after  a 
DC  run  without  this  card  if 
ESP0DEPH  and  a  covariance 
matrix  update  is  flagged  on 
the  JDC. 


Changes  a  particular  data  (R , 
A,  E,  R,  a,  6)  sigma  for  a 
particular  sigma  type  (01, 
02,  •  ■  •  ,  60) 

Changes,  for  the  given,  run, 
the  stations  data  sigma, 
surface  r efrac tivity ,  and/or 
the  gross  outliers  criterion 

Enters  atmospheric,  drag 
models,  CpA/2m,  and  K 
specified  on  card 

If  a  single  input  for  day  of 
epoch,  will  be  used  as  an 
average  value;  if  daily 
values  desired,  and  a  short 
set  is  input,  last  value  will 
apply  to  all  subsequent  days. 

Enters  parameters  as  indi¬ 
cated  by  a  flag  in  the  proper 
column  of  each  card 

Uses  radiation  pressure 
model,  with  vehicle  param¬ 
eters  as  specified 

An  a  priori  "S"  matrix  is  in¬ 
put  to  use  data  previously 
entered  into  the  solution 

If  this  is  an  ESP0DEPH  run 
only  (no  ESP0DDC),  the 
a  priori  covariance  matrix 
punched  in  UPMAT  is  input 
and  updated.  If  ESP0DDC 
is  present,  cards  are 
ignored. 
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Table  2 -III.  Program  Response  to  Omitted  and  Included 
Input  Cards  (Continued) 


Input 

DAC 


PRDCT 


DELTT 


99  CARD 


No  Card  Present 
in  Deck 


No  ESP0DEPH  run,  unless 
DELTT  or  PRDCT  card(s) 
present 

No  ESP0DEPH  run,  unless 
DELTT  or  DAC  card(s) 
present 

No  ESP0DEPH  run,  unless  a 
PRDCT  or  DAC  card(s) 
present 

No  changes  to  program 
constants 


Card  Present 
in  Deck 


Trajectory  points  predicted 
for  specific  times  on  DAC 
cards  input 

Trajectory  points  predicted 
for  specific  times  input 

Trajectory  points  printed  as 
input  for  At  increments 
over  interval  T  minutes 
from  epoch 

Changes  a  given  constant 
specified  by  a  core  location 
identification  to  the  number 
input 
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Inclusion  of  radiation 
pressure 
Constants  change 


Table  2-V.  Input  Diagnostic 


INPUTS 

MANDATORY 

OPTIONAL 

Card  Input 

ESPOD 

Nominal 

SPADATS 
Tape  Files 

Overriding 
Card  Input 

Options  Called 
by  Cards 

Optional  Input 
Data  on  Cards 

B2  SYSTEM  START 

70ASCHTP  ** 

X 

2 

JOB  NO 

X 

Remark*  Card{s) 

X 

RUN  ESPOD,  DATA 

X 

ESP0DDC  and  ESP0DEPH 

JDC 

X 

N.  A. 

Either 

1CTYP  ^ 

X 

N.  A. 

ICOND  > 

X 

N.  A. 

ITIME  ) 

X 

N.  A. 

or 

1CTYP  'l 

X 

N.  A. 

7-card  el  / 

X 

N.  A. 

DNREV  ) 

X 

N.  A. 

X 

or 

NO  ICTY  P'1 

SEA1 

DNREV  J 

Time  of 

X 

last  obs. 

DRAG 

None 

X 

APF10 

N.  A. 

X 

TSTEP 

1  min 

X 

RADPR 

None 

X 

ZONAL 

J2»  J4 

X 

SECTR 

None 

X 

TESSR 

None 

X 

ESP0DDC 

CAT  1 

ADBARV 

X 

CAT2 

None 

X 

BNDS 

Nominal 

Bounds 

X 

BISES 

N.  A. 

X 

SMAT 

N.  A. 

X 

NITER 

5  iterations 

X 

DELET 

N.  A. 

X 

TMAX 

10  days 

X 

SIGMA 

SIGMA 

TABLE 

X 

STYPE 

SENSOR 

TABLE 

X 

ESPODEPH  only 

Either 

DAC 

X 

or 

PRDCT 

X 

or 

DELTT 

X 

UPMAT 

N.  A. 

X 

GENERAL 

REM 

X 

99 

N.  A. 

X 

ENDPR 

X 

N.  A. 

ESPODDC  -  SPADATS 

File  Data  Only 

0BS  CARDS 

N.  A. 

SRADU 

ENDOB 

X 

N.  A. 

SEN  CARDS 

N.  A. 

SEAI 

ENDSN 

X 

N.  A. 

GENERAL.  FOLLOWING 

LAST  CASE 

ENDAT 

X 

N.  A. 

B2  SYSTEM  FINISH 

ENDDATA 

X 

11 

ENDOFJOB  8 

X 

2 

11 

ENDSCHED  8 

X 

2 
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2.  2.  3  Load  Sheets 


ESPOD  Load  Sheets  have  been  prepared  to  facilitate  the  task  of 
punching  input  cards.  Sample  load  sheets  with  the  variable  field  con¬ 
tent  indicated  have  been  included  in  the  section. 

The  load  sheet  has  abbreviations  for  the  flagged  columns,  i.  e.  , 
columns  30  to  56.  They  are  listed  below. 

Interpretation 


olumn 

Abbreviation 

30 

cs 

31 

PR 

32 

PS 

33 

P  0 

34 

SC 

35 

0C 

36 

PC 

37-40 

Not  used 

41 

DC 

42 

SM 

43 

PI 

44 

PS 

45 

R  <£> 

46 

NE 

47-50 

Not  used 

51 

P  0 

52 

XT 

53 

EP 

54 

Not  used 

55 

UP 

56 

PS 

Cold  start 

Preprocessing  necessary 
Print  sensors 
Print  observations 
Sensor  cards  present 
Observation  cards  present 
Print  program  constants 

Differential  correction 
S -Matrix  input 
Punch  AXA  inverse 
Punch  ATA 
Residual  observations 
Use  new  elements 

Postprocessing  required, 
i.  e.  ,  execute  ESP0DEPH 

DAC  cards  present  this  run 

Generate  ephemeris 

Perform  matrix  update 

Punch  inverse  of  update 
covariance  matrix 


The  actual  column  content  is  given  in  Section  2.  1.  2. 
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DATE-  - 

NAME 

MNOSkEM  NO. 


ESP0D  INPUT  SHEETS 


Note  to  k.ypuncfrer 

Card  having  columns 
rule  are  cards  hav: 


10  through  00  blank  should  not  be  punched.  Exceptions  to  this 
ng  ENDPR,  BND0B  or  BOAT  Id  the  first  name  field. 


MNIONlT* 
KEYFUMCMEO  SV 


'  11 1 1 

■«  1  *  1  •  1 1 

•  1  •  1101111*211  JIM  1 14' 14, 17  (14  !*«  20]21!27i21!24(24(24i77!  24i2*  1*1  J1 

17  jSS  !*4|7M*!  171  m  >*!4*l  4<J47  14  *]44)4SL4*,  *7J*t  '4r  so  %\ 

42.  47  74  44  44[47]t*  74  *0-4  1  *7 ■  4  >|  4  4(44(44 |«7L4«|Mt 7»l7*  |  7 2,  74l  T4|  74l  74l  77J74  (7»  I»0 

j,nc 

Veh. 

Vehicle  Kune 
i  i  i  i  i  1  i  i  i 

Reader 

I  L  1  1  1  1  1  1  1  1  !■ 

C 

S 

P 

R 

P  P  S  fT  P 

i  i  i 

D 

C 

,  1 1 

i 

M.  V  , . ,. 

c 

i  i  i  t  i  i  i  i 

i  i  i  i  i  i  i  i 

i  i  i  i  i  i  i 

1  1  1  1  1  L i L 

1  1 -1—1  1  1  1  1  1  1  1  1  1  1  1  1  1  1 

9eq 

Np 

It«- 

Haste 

till 

Variable 

L 1  1  1  1  1  i  1  LJ^i  i  1 

Name 

iiil 

Variable 

till  11  1  LI  1  1  1  ' 

Has*. 

III! 

Variable 

1  I  I  1  i  i  i  1  I  1  1  i  i 

Name 

Variable 

1  i  1  1  l  i  i  i  l  l  1  1  1 

R.A*  . 

7WrS  XS  sf  A 

1  1  1  1  1  1  1  1  l  1  1  1  ! 1  1  1  1  1  1  1  1  1  1  1  1  <1  1  1  1 

:<s?s?o  * 

i  i  i  i 

t  ;  i  i  i  t  i  i  i  i  i  i  i 

iiii 

i  i  i  i  i  i  i  i  i  i  i  i  i 

L 

WW 

/-  or  J2.  or  S  f  r 

1  1  1  1  1  1  1  1  1  1  1 

till 

1  1  l  1  I  I  1  1  L  1  1L  L 

iiii 

l  l  :  l  i  i  l  l  1  1  1  1  i 

i  i  i  l 

iiiiiiiiiiiti 

Ojl 

I.C,ril,D 

oc  or  ^  or  X 
lilt l  l l  1  i 1  1  1  i 

$  Or  $ 

1  1  1  1  1  1  1  1  l  l  1  l  1 

iiii 

j#  or  Z 

i  l  i  I  I  1  i  i  ,  I  1  .  l 

iiii 

/)  Or  X 

1  I  I  1  *  i  •  j  i  I  »  «  » 

°P 

1 

I,C,A«,D 

G  or  £ 

1  1  1  1  1  1  T  1  1  1  I  1  1 

i  i  i  i 

V  or  X 

l  i  i  i  l  t  1  1  l  l  l  l  l 

iiii 

*  l  1  1  1  1  l  l  i  1.  1  1 

llll 

1  1  1  1  1  1  1  1  I  1  1  1  1 

OjL 

J 

Y£*4. 

i  i  i  i  i  i  i  i  i  i  i  i  i 

lilt 

/>;  a  /y  r// 

ii(iiitiiiii 

iiii 

I  1  I  I  i  i  i  i  i  1  1  1  1 

l  i  l  i 

1  1  1  1  1  1  1  1  1  1  - 1  -l i 

| 

^1*.  w 

/ritflOTt 

l  l  1  1  l  i  l  l  l  l  l  l  i 

lilt 

ScCQ  A/D 

\  i  l  l  1  l  l  l  l  I  l  l  l 

1  I  1  1  11  1  1 1  1  l  l  l 

iiii 

1  1  1  11-1  1   1- 1   1  1  1  1  . 

i 

D.,R(  A,Q( 

A /z  At 

iiiiiiiiiiirt 

i  i  i  i 

K 

i  i  i  t  i  i  i  i  i  i  i  ii 

iiii 

£>r*n  cptn'H  /  cr  X. 

1  i  l  l  i  i  i  i  i  i  i  i  i 

iiii 

At  >***  AC.r  J  rrS  4 

T1  1  1-  l  -1-1  ill  1  1  1 

i 

Bi  "i  D.  ^ 

r/v 

A  1  A  L  1  1  1  1  1  1  1  1  1 

or 

iill 

C/frz  a*d  C/frz 

ii  i  i  i  i  i  i  i  i  i  i  i 

5Vc lc  rn 
iiii 

'/V  V /I  *>/)!?££  5 

1  1  1  1  1  1  1  1  1  1  1  1  1 

. I  t  1  1  ]  1  1  1 

°F 

i 

>l!M 

KCtc) 

1  1  1  1  1  l  1  l  1  1  1  1  1 

iiii 

(e'Zt  1 

i  i  l  i  l  l  l  l  1  1  1  l  1 

iiii 

1  1  1  I  l  I  1  1  1  1  1  1  1 

l  l  t  l 

1  1  1  1  1  f  1  1  1  1  1  1  1 

i 

c^^i, 

ft  ,  •  ,  , 

l  i  i  i 

i  l  i  i  1  1  11  J  1  I  1  1 

iiii 

i  i  i  i  i  i  i  i  i  i  1  i  1 

llll 

1-  1  I  1  1  1  1  1  1  1  1  1  1 

i 

■iiw 

H.O 

l I  1  l.LU  LI  11  l  1 

111! 

1  1  l  1  1  1  1  1  1  1  1  1  1 

Iiii 

1  1  1  1  1  1  1  1  1  1  i  i  i 

iiii 

1  1  1  1  1  1  1  1  -J  1  1  1  1 

i 

ViAV 

1  1  1  1  1  1  1  .1  1  1  1  1  1 

1  l  ll 

CA  or  <7^ 

l  l  l  i  i  l  i  l  !  1  l  i  1 

-1..1  i  i 

CTg  or 

i  i  i  i  i  i  i  r  i  j  i  i  i 

i  i  i  i 

1  1  1  1  1  1  1  1  1  L  L   1 — 1 

i 

ViW 

iiiiiiiiiiiti 

1111 

(Wc) 

i  i  i  i  i  i  i  i  i  l  i  i  l 

IIII 

iitiiit  iiiiii 

llll 

(«•  Aj 

1  1  1  1  1  1  1  1  1  1  1  1  1 

i 

R-l1!  °lMl  A 

0 etc) 

i  i  i  i  i  i  i  i  i  i  i  i  i 

1  1  J 1 

(<?  te) 

i  i  i  l  i  i  i  i  i  i  i  i  l 

iiii 

*<) 

i  i  i  i  >  i  i  i  i  ii  ii 

llll 

1  1  1  1  1  J  i  J  1  1  1  1  1 

1  1  1  1 

i  i  i  i  t  i  i  i  i  i  i  i  i 

i  i  i  i  i  i  i  i  i  i  i  l  i 

lilt 

i  i  i  i  i  i  i  i  i  i-  i  i  i 

iiii 

1  1  1  1  1  1  1  1  1  1  1  I  1 

*,1  DP,R 

1 1  * 

rrhn 

•  It  1  1  Is  I* 

74  »0 ' 4  1  1  42 1  4  J  J  J  JJl  **l4Ti  *H>»1  4 1*4» 

fnfefiiferhri 

4«i  4 J  soUrwUj'to  mUoUi 

4 2 It 4144 144 1 44  4»T40T4«  W 2 1171  IT 4 174 17»1«  177 IT*  17*  1*0 

oati  ESP0D  INPUT  SHEETS  e*oi  of 

NANf  Kot>  to  kpvnunch«r!  priori  tv 

MOlLtM 

NO,  OF  CA 

NO. 

Cards  having  columns  10  through  00  blank  should  not  be  punched.  Exceptions  to  this 
rule  are  cards  having  ENDPR*  KNDffe  car  ENDAT  In  the"  first  name  field. 

VI 

!VPUNCHBD  RV 

ROI _ _ 

[R1ME0  RY    

im  * 

4  14  !4  17 

t  [  4  11*111  |iai4|14l181Hil7 

14  114120  12 1:221 24! 24(24! 24127  2*  24' 40 

I, 

12144,44  (  4414*  ;47|»:7*(4e 

41 

«7|«9]44  dt  ;44j4  7»««{4*;SC  S*  S2.  $  ) .  S4( t* i  St  ] S  7  *  *  5»  4QC  1  rl?  L  •»i44i4t14«|«7|4C  4»  70 ,7!  |  72|  1 1|  74;  7»i  7#|  7 7l?«  {7t  |40 

AM 

Veh. 

M'l 

Vehicle  Haute 

i  1  i  1  1  i  1  1  1 

Header 

11111,11111 

C 

s 

p 

R 

p  p  sT  p 

3MC.C,  ,  ,  , 

D 

C 

mi, 

s 

if,  ,n, 

L  llll  1  1  1  1  1  1  1  1  1  1  1  1  1  1 

Ji  R 

C 

A  1 

1 

1  i  1  1  1  1  1 

iiiiii 

i  i  i  i 

1 1  1  1.1-1  1  1 

iiii 

AIL  1_ 1_ 1  L 

l  1_ 1_ 1_ 1_ 1  LI  Ii  1  II  1_ 1_ 1_ 1 

Seq 

*P 

Iter 

Np 

i  i 

Variable 

-Li  1  1  1  1  1  i  1  1  1  1  1 

Name 

Iill 

Variable 

:  i  i  i  i  i  i  i  i  I  i  i  i 

Name 

(II! 

Variable 

i  i  t  i  i  i  i  i  i  i  i  i  i 

Name 

L  L  1  1 

Variable 

1  i  i  i  i  i  i  i  i  l  l  l  i  t 

L 

Jl i 

L.l  1  1  1  1  1  1  1  1  i  1 

till 

THIS  IS  A  SAMPLE  I 

i  i  i  i  i  i  i  t  t  i  i  i  i 

-,OAD  SHI 
i  t  t  i 

EET 

l l 1  l  1  1  1  1  1  1  i  1  1 

i  1  1  1  1  1  1  1  1  1  1  1  1  I  1  -L  1  1 

i 

WW 

1 .  or  2.  or  3.  or  4, 

1  1  l  l  1  1  L 1  1  1  1  1  1 

iiii 

i  i  i  i  i  i  i  i  i  i  i  i  i 

iiii 

i  j  l  i  l  i  ii  i-i-i  i  .1 

lilt 

i  i  i  i  i  i  i  i  i  i  i  t  i 

i 

i.c,  *Ud 

a  or  X  or  x 

l  i  i  l  i  l  i  i  i  i  i  i  i 

1 l L  1 

6  or  y 

I  L  1  l  1  1  1  1  1  1  1  1  1 

iiii 

p  or  z 

1  1  I  1  1  1  !  1  1  1  1  1  l 

iiii 

A  or  x 

i  i  i  i  i  i  i  i  i  i  i  i  i 

P? 

i 

I,C,  ^H,D 

R  or  y 

i  i  i  i  i  i  i  i  i  i  i  i  i 

-1  L  L  1 

v  or  z 

I  i  l  i  i  i  i  i  i  i  i  i  i 

iiii 

i  i  i  i  i  i  i  i  i  i  I 

llll 

i  1  i  l  1  1  1  i  1  1  1  1  1 

i 

YEAR 

I  i  i  i  i  i  l  i  i  i  i  i  i 

llll 

M0NTH 

i  i  i  i  i  i  t  i  i  i  i  i  i 

iiii 

DAY 

i  I  1  1  1  1  4  1  1  1  i  1  1 

llll 

H0UR 

i  i  t  i  i  i  i  i  i  i  i  i  i 

1 

MINUTE 

i  i  i  1  1  l  i  l  l  l  l  l  l 

lilt 

SEC0ND 

i  i  i  i  1  1  1  I  1  1  1  1  1 

llll 

i  i  i  i  i  i  i  i  l  l  i  i  i 

llll 

i  i  i  i  i  i  i  i  i  i  i  i  i 

i 

d.R,  A°. 

Cj-jA/2m 

■  i  i  i  l  i  i  i  i  l  i  l  i 

llll 

K 

i  i  i  t  i  i  i  i  i  i  i  i  i 

llll 

Drag  option  1.  or  2. 

i  i  i  i  i  i  i  i  i  i  i  i  i 

11  11 

3p?5».  ?r/ 

| 

Bi  R| D,  ^ 

A  1  1  1  1  1  1  1  1  1  1  1  1 

IN  0R 

llll 

DER  0F  CAT  1  AND  CAT 
l  l  1  l  l  1  1  1  1  l  l  t  l 

:  S0LUT1 
llll 

0N  VARIABLES 

1  i  1  1  1  i  1  l  i  1  1  i  i 

1  1  i  1 

1  1  1  1  1  1  1  1  1  1  !  1  1 

°? 

1 

(etc) 

1  i  l  1  1  l  i  i  i  l  l  1  i 

iiii 

(etc) 

i  i  i  i  i  1  l  i  i  i  1  i  i 

iill 

1  l  i  i  1  i  i  l  i  l  l  l  1 

llll 

11  1  1  1  1  1  1  1  1  1  1  1 

<¥■ 

1 

aif  .  a  A,R,v  ^Kl  1  ,  .  1  i 

iiii 

i  i  i  i  i  t  i  i  i  i  i  i  i 

llll 

i  i  i  i  i  i  i  i  i  i  i  i  i 

llll 

1  1  1  1  1  1  1  1  1  1  1  1  1 

1 

N,  0 

-i  l  l  l l l  1  1  l  1  1  l  1 

III! 

1  1  1  i  i  1  1  1  1  1  1  i  1 

llll 

1  1  1  i  i  1  i  i  i  l  i  i  i 

llll 

11111,1111111 

1 

*R  or  ra 

i  i  i  l  1  l  l  i  l  i  l  l  l 

1111 

tA  °r  *6 

i  i  i  i  i  i  i  i  i  i  i  i  l 

llll 

<r  „  or  <r 

E  Q 

i  i  i  i  i  i  i  i  i  i  i  i  i 

llll 

tr  or  % 

i  i  i  i  i  i  i  t  l  i  i  i  i 

i 

(etc) 

i  i  i  i  i  i  i  i  i  i  i  i  i 

1  1  1  l 

(etc) 

i  i  i  i  i  i  i  i  i  i  i  i  i 

llll 

(etc) 

i  i  i  i  i  i  t  i  i  i  i  l  l 

till 

(etc) 

l  l  l  1  1  1  1  1  t  1  1  1  1 

1 

*A<k*i* 

(etc) 

llll  ill  1 L  1  1  1  1 

llll 

(etc) 

l  1  1  I  l  i  1  1  1  1  1  1  1 

llll 

(etc) 

i  i  i  i  i  i  i  i  i  i  i  i  i 

llll 

(etc) 

i  i  i  i  i  i  i  l  l  i  i  l  i 

iiii 

1  l  1  1  1  1  1  1  1  I  1  1  1 

l  1  1  1 

i  i  i  i  i  i  i  ■  1  i  l  ■  ■ 

l L i 1 

1  1  1  1  1  1  1  1  1  L  1  1  1 

llll 

i  i  i  i  i  i  l  l  i  i  i  i  i 

IK,  DPR 

l  1_ 1_ 1_ 1_ 1_ 1_ 1_ 1  i  1  ft  1 

1  'I  4  U 

4  It  1  7l*  la  14(  111  1  21 1  41 1  411 4  M4  ri  711  4114  1  >0' 2  l’  221 2 1  2  4 ' 7V  24  2?  24  24^  »071  n42l  1  >  4M  ,V  TtliTI  >tT»7T  404  T47  4  jT44T4sT4V  47 

4*  4*12  ,  '42,4|!4«-44'44  '*7'44'4*l46l4i  ilbllllltlltl  4  Ttttl  4  4l  TOf?  1 171T7 4 1741701 74 1711 7*174  1*0 

Figure  2-64. 


Sample  Load  Sheets 
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igure  2-65.  Sample  Load  Sheet 


2.  2.  4  Sample  Input  Decks 


The  ESPQ)D  input  deck,  Figure  2-66  illustrates  the  data  arrangement 
for  an  orbit  determination  run.  The  deck  setup  for  a  run  involving 
ESP0DEPH  only  is  the  same  except  that  observation  cards  and  sensor 
cards  are  not  present. 

The  input  deck  as  shown  is  the  order  in  which  the  various  input 
categories,  referred  to  as  subdecks,  are  arranged.  The  content  of  the 
preliminary  data  subdeck  determines  the  nature  of  the  ESP($D  run. 

Various  possible  versions  of  this  subdeck  will  now  be  considered. 

Reference  to  Sections  2.  2.  2  and  2.  2.  3  will  help  ascertain  the 
functions  of  the  subdecks.  The  content  of  the  associated  JDC  card  will 
be  given  with  each  preliminary  subdeck,  since  it  is  the  control  card  for 
the  program.  The  contents  of  the  first  29  columns  of  the  JDC  card  have 
been  omitted  from  the  card  content  description  as  it  is  not  relevant  to 
the  discussion.  The  End  of  Preliminary  Data  Card  (ENDPR)  was  excluded 
from  the  sample  decks  which  follow  for  the  same  reason. 

Case  1 

The  simplest  differential  correction  run  consists  of  a  JDC  card 
only.  File  elements  and  sensor  data  are  obtained  from  the  SEAI  tape, 
and  observations  obtained  from  the  SRADU  tape.  In  this  case,  the  nomi¬ 
nal  epoch  is  the  last  time  of  observation.  A  DNREV  card  could  override 
this  epoch  selection.  Since  no  solution  vector  is  specified,  a  nominal 
6x6  fit  (first  six  variables  of  CAT1)  with  nominal  bounds  is  solved. 

This  is  a  drag  free  model  with  a  nominal  (Model  1)  potential  model. 

There  is  a  maximum  of  five  (5)  iterations  and  observations  which  are 
more  than  ten  days  from  epoch  are  rejected.  These  nominal  values 
could  be  changed  with  a  NITER  card  and  a  TMAX  card,  respectively. 

The  output  from  this  run  (see  Section  4.  0  for  detailed  discussion 
of  output  formats)  consists  of  station  number  and  name,  observation 
time  (from  epoch),  numbered  residual,  a  summary  table  of  the  mean 
and  RMS  residuals  for  each  station,  the  delta  correction  to  each  variable, 
the  old  value  of  the  variable  (initial  estimate),  the  new  value,  the  uncer¬ 
tainty,  the  bound,  the  current  RMS,  predicted  RMS,  best  RMS  so  far, 
the  covariance  matrix,  the  correlation  matrix  and  a  message  to  indicate 
whether  the  solution  is  converging. 
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Figure  2-66.  ESP0D  Input  Deck 
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Figure  2-67.  Load  Sheet  for  Problem 


The  critical  columns  of  the  JDC  card  for  this  particular  run  are 
31,  34,  35,  and  41.  The  load  sheet  for  this  case  follows,  consisting 
of  a  JDC  card  only. 

Case  2 

The  simplest  card  input  subdeck  for  an  ESP$DEPH  run  consists 
of  two  cards.  The  run  consists  of  an  ephemeris  prediction  only,  using 
file  elements  and  epoch. 

Tapes 


1. 

JDC  Card 

1. 

SEAI 

2. 

PRDCT  Card 

2. 

SRADU 

Since  there  is  no  specification  of  initial  conditions  on  cards,  the 
initial  conditions  will  be  read  off  the  SEAI  tape,  in  the  SPADATS  mean 
elements  format.  The  epoch  revolution  number  is  specified  in  the 
SPADATS  mean  elements,  and  this  is  the  epoch  which  is  used  in  this 
case.  Note  that  the  radiation  pressure  model  and  the  atmospheric  models 
will  not  be  included,  and  only  a  nominal  potential  model  (J^>  J J^)  is 
used  in  the  ephemeris  prediction. 

A  variation  of  this  case  would  be  to  include  a  DNREV  card  for  the 
purpose  of  updating  the  file  elements  to  a  desired  epoch. 

The  ESPC^DEPH  output  for  each  time  specified  is:  vehicle  position 
and  velocity  in  polar  and  Cartesian  form,  the  classical  orbital  elements, 
vehicle  altitude,  geodetic  latitude,  longitude,  argument  of  latitude,  period, 
perigee  and  apogee  altitudes,  time  from  epoch  to  the  next  crossing  of  the 
ascending  node,  and  a  set  of  indeterminancy  free  elements.  If  column  53 
of  the  JDC  card  is  punched,  the  special  binary  ephemeris  tape  containing 
x,  y,  z,  x,  y,  z  versus  time  is  generated. 

The  load  sheet  for  this  run,  consisting  of  two  cards,  is  shown  in 
Figure  2-68. 
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Figure  2-68.  Load  Sheet  for  Problem 


Case  3 


The  input  sequence  of  this  case  is  given  below. 


1.  JDC 


2.  ICTYP 

3.  IC0ND 

4.  ITIME 


Tapes  Needed 

1)  SEAI 

2)  SRADU 


5.  DRAG 


6.  APF10 

7.  CAT2 

8.  BISES 

9.  STYPE 
10.  SIGMA 

This  is  preliminary  subdeck  setup  (except  for  the  JDC  card)  for  a 
differential  correction  run.  The  initial  conditions  are  specified  on  cards 
and  the  solution  vector  consists  of  Category  2  variables  only.  The  pre¬ 
sence  of  a  CAT2  card  without  a  CAT1  card  suppresses  the  nominal 
6x6  solution  vector.  The  BISES  card  contains  the  initial  bias  estimates 
of  the  CAT2  variables.  The  only  non-nominal  perturbation  in  this  inte¬ 
gration  model  is  the  dynamic  atmosphere,  which  requires  APF10  cards. 
The  STYPE  card  is  used  to  change  sensor  characteristics  and  the  SIGMA 
card  changes  standard  deviations  in  the  SIGMA  table.  (The  load  sheets 
for  this  case  are  shown  in  Figures  2-69  and  2-70. 


Card  1 :  JDC 


Card  2-3:  IC0ND 


Specifies  the  initial  conditions  in 
Cartesian  coordinates 


Card  4-5:  ITIME 


Defines  time  of  epoch 


Card  6:  ICTYP 


The  M3."  Specified  type  3  elements; 
i.  e. ,  Cartesian. 


Card  7:  DRAG 


The  first  two  variable  fields  contain 
the  drag  parameter  and  the  drag  vari¬ 
ation  respectively.  The  third  variable 
field  (content:  2)  calls  for  secular 
drag  variation;  and  the  fourth  variable 
field  (content:  4)  flags  the  COESA 
Dynamic  Atmosphere. 
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Card  8:  SIGMA 


Cards  9-10:  APF10 


Card  11:  CAT2 


Cards  12-14:  BISES 


Card  15:  STYPE 


This  card  is  used  to  change  the 
standard  deviations  of  sensor  type  17, 
indicated  in  the  first  two  columns  of  the 
card.  The  four  variable  fields  contain 
the  new  standard  deviations  for  this 
sensor  type  (for  this  run  only). 

These  cards  contain  the  Ap  (geomag¬ 
netic  planetary  amplitude  )  and  FlO  (10.7 
centimeter  solar  flux)  coefficients  for 
the  days  covering  the  period  of  integra¬ 
tion.  (8  days) 

The  first  four  columns  of  the  variable 
field  identify  the  sensor  (0329).  The 
,flMs  call  the  variables  to  be  solved 
for.  For  Sensor  No.  0329  they  are: 
range  bias,  azimuth  bias,  elevation 
bias,  range  rate  bias,  latitude  bias, 
and  longitude  bias.  In  the  second 
variable  field,  the  solution  variables 
for  Sensor  No.  0330  are:  azimuth 
bias,  elevation  bias,  range  rate  bias, 
and  time  bias. 

These  three  cards  contain  the  bias 
estimates  for  the  Category  2  variables. 
As  indicated  on  the  load  sheet,  the  bias 
estimates  correspond  one-to-one  with 
Category  2  variables. 

This  card  is  used  to  change  data  in 
the  STYPE  table;  that  is,  the  sigma 
type,  the  gross  outliers  criterion, 
the  refraction  flag,  and  the  mean 
surface  value  of  refractivity .  The  new 
parameters  are  valid  on  this  run  only. 


Cas  e  4 


This  case  is  a  combined  ESP0DDC  and  ESP0DEPH  run,  and  con¬ 
tains  most  of  the  options.  The  input  preliminary  cards  are: 


1. 

JDC 

13. 

CAT  1 

2. 

REM 

14. 

CAT2 

3. 

ICTYP 

15. 

BNDS 

4. 

SPADATS  -  7  element 

16. 

BISES 

5. 

DNREV 

17. 

NITER 

6. 

DRAG 

18. 

DELET 

7. 

APF10 

19. 

TMAX 

8. 

radpr 

20. 

SIGMA 

9. 

Z0NAL 

21. 

STYPE 

10. 

SECTR 

22. 

PRDCT 

11. 

TESSR 

23. 

99  cards 
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Figure  2-69.  Load  Sheet  1  for  Problem 
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Figure  2-70.  Load  Sheet  2  for  Problem 


The  initial  conditions  are  specified  on  the  SPADATS  mean  element 
cards  (7),  with  the  epoch,  as  specified  with  these  elements,  to  be  updated 
to  a  desired  epoch  with  a  DNREV  card.  As  in  the  previous  case,  a  dy¬ 
namic  atmosphere  requiring  A^  and  data  is  selected. 

The  initial  integration  time  step  is  changed  as  well  as  the  maximum 
number  of  iterations  in  the  differential  correction  process.  The  most 
sophisticated  potential  model  as  well  as  the  inclusion  of  radiation  pres¬ 
sure  complete  this  perturbation  model.  The  solution  vector  consists  of 
both  CAT1  and  CAT2  variables.  As  in  the  previous  case,  initial  esti¬ 
mates  of  the  Category  2  variables  are  included.  New  bounds  for  this 
case  have  been  introduced.  Data  editing  in  the  form  of  residual  selection 
and  rejection  of  observations  too  far  removed  in  time  from  epoch  is 
included.  Data  in  the  sigma  and  STYPE  tables  are  changed  for  this  run. 

The  ephemeris  generation  is  specified  on  a  PRDCT  card,  although 
it  could  be  specified  on  DAC  or  DELTT  cards.  The  load  sheets  for 
this  run  are  shown  in  Figures  2-71  and  2-72. 


Card  1 :  JDC 


Card  2:  REM 


Columns  41  and  51  call  for  the 
differential  correction  package  and 
ESPC^DEPH.  Column  36  indicates 
that  the  program  constants  are  to  be 
printed.  In  column  46,  the  "newest 
elements"  are  selected  to  generate  the 
updated  times  in  ESP0DEPH.  These 
are  the  elements  whose  residuals 
have  not  been  calculated.  Column  53 
calls  for  the  generation  of  the  spe¬ 
cial  binary  ephemeris  tape. 

The  remarks  card  need  not  follow 
the  JDC  card. 


Card  3:  ICTYP 


Card  4:  DRAG 


The  content  (4)  indicates  that  type  4 
elements  are  read  in  as  initial  condi¬ 
tions  (SPADATS  mean  elements). 

Same  as  Case  3 


Card  5-6:  BNDS  These  two  cards  contain  the  bounds 

for  the  Categroy  1  and  Category  2 
variables  respectively. 

Card  7:  CAT1  The  Category  1  variables  with  a  "l" 

indicate  inclusion  in  the  solution 
vector.  They  are:  right  ascension, 
declination,  flight  path  angle  ((3), 
azimuth,  C^A/2m,  and  K. 
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Card  8:  NITER 


Card  9:  SIGMA 


Card  10:  STYPE 


Cards  11-12:  APF10 


Card  13:  DNREV 


Card  14:  RADPR 


Card  15: 

Z0NAL  >. 

Card  16: 

SECTR 

Card  17: 

TESSR 

Card  18: 

99  card 

Card  19: 

CAT2 

Card  20: 

BISES 

This  permits  the  maximum  iterations 
in  the  differential  correction  process 
to  be  8.  The  nominal  value  is  5 
iterations. 

The  first  two  columns  indicate  the 
sensor  type  (17),  and  the  variable 
fields  contain  the  new  standard 
deviations  for  that  sensor  type. 

This  card  can  change  the  remaining 
sensor  characteristics  of  sensor 
type  17;  that  is,  gross  outliers  cri¬ 
terion,  refraction  flag,  and  mean 
surface  value  of  refractivity. 

These  two  cards  contain  the  Ap  and 
F  iq  coefficients  for  the  period  of 
integration. 

The  content  (-1.  0)  is  a  flag  to  update 
the  epoch  (when  using  SPADATS  mean 
elements)  to  the  last  time  of 
observation. 

This  calls  the  radiation  pressure 
option.  The  variable  fields  contain 
the  effective  area  of  the  vehicle  in 
square  meters,  and  its  mass  in 
kilograms. 


These  four  cards  specify  the  fourth 
potential  model  as  described  in 
Section  2.  1.3.  7.  3 


The  Categroy  2  variables  selected 
are  the  latitude  and  longitude  biases 
for  sensor  No.  0330.  Hence,  the 
solution  vector  of  this  run  is  an 
8x8  system. 

This  card  contains  the  bias  estimates 
of  the  Category  2  variables,  in  the 
order  given.  That  is,  the  first  vari¬ 
able  field  contains  the  latitude  bias 
and  the  second  field,  the  longitude 
bias. 
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Card  21:  DELET 


Card  22:  TMAX 


Card  23:  PRDCT 


Cards  24-30: 


2.  2.  5  Stacked  Cases 


This  card  specifies  the  first  and 
last  residual,  in  each  V-field,  of  a 
list  of  residuals  to  be  deleted.  The 
identification  numbers  are  obtained 
from  the  residuals  listing  output  after 
running  at  least  once  with  this  case. 

This  card  permits  the  rejection  of 
observations  which  are  more  than 
8  days  from  epoch.  The  nominal 
value  is  ten  days. 

This  card  specifies  the  update  times 
at  which  the  vehicle  position  and 
velocity  are  calculated  and  printed. 
These  times  could  have  been  specified 
on  DAC  or  DELTT  cards. 

SPADATS  mean  elements.  See 
Section  2.  1.  3.  1.  1  for  format. 


It  is  possible  to  stack  cases  for  a  single  run.  For  each  case  to  be 
run,  an  adequate  set  of  case  description  cards  (b,  c,  d,  e)  are  prepared. 
Cards  for  input  are  in  seven  categories  as  listed  below: 

1)  Job  start  in  B2  system  schedule  mode 

2)  ESP(#D  job  description  card,  JDC 

3)  Preliminary  data 

4)  Observation  cards,  if  any 

5)  Sensor  cards,  if  any 

6)  ESP$D  end  of  data  (ENDAT) 

7)  Job  finish  in  B2  system  schedule  mode 
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Figure  2-71,  Load  Sheet  1  for  Problem 
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igure  2-72.  Load  Sheet  2  for  Problem 


The  case  description  cards  are  arranged  in  multiple  consecutive 
sets  between  the  B2  system  Job  start  and  Job  finish  cards  (see  Section 
2.  1.  1),  For  example,  the  card  arrangement  for  three  cases  would  be: 

a 

bi  ci 

b2  c2 

b3  c3 
f 

g 

When  running  with  File  Data,  the  observation  and  sensor  data  are 
obtained  from  the  SRADU  and  SEAI  tapes,  respectively;  this  would  elimi¬ 
nate  the  d  and  e  cards  from  the  above  sequence. 
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3.  OPERATIONAL  INSTRUCTIONS 


ESPOD  may  be  executed  in  two  operational  modes  in  the  BZ  system; 
the  manual  mode  and  the  schedule  mode  as  discussed  below. 

3.  1  MANUAL  MODE 

Table  3-1  lists  the  required  and  optional  tapes  for  the  manual  mode. 


Table  3-1.  Tape  Setup  for  Manual  Mode 


Logical 

Tape  Number 

Setup 

Tape  Description 

® 

Input 

Prestored  input  data 

i 

System 

BZ  master  system  tape 

4 

SEAI 

Master  SEAI  (sensor,  elements, 
acquisition,  and  information  sites) 

6 

SRADU 

Observation  tape 

© 

Scratch 

ESPOD  writes  blocks  of  C$MM$N 
data  and  observations  on  this  tape 

<°> 

Planetary  ephemeris  tape 

Scratch 

_ ^  ^ 

Special  binary  ephemeris  tape  (R,  V 
versus  time)  generated  by  ESP0DEPH 
on  option 

© 

OUTPUT 

Off  line  output  tape 

<0> 

Prestored  octal  change  cards 
(TOGGLE  45,  ON,  if  octals  present) 

Write  ring  required 


May  be  left  out  if  proper  conditions  are  met 
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When  the  B2  System  has  been  initiated,  giving  control  to  the  flexo- 
writer,  the  following  is  typed  out: 

MANUAL  Typed  in  by  operator 

When  the  manual  mode  is  indicated  as  ready, 

ESP0D  Typed  in  by  operator 

The  computer  system  will  locate  the  ESP$D  segment  on  the  B2  sys¬ 
tem  tape  and  transfer  control  to  it.  The  following  will  then  by  typed  by  the 
System  on  the  flexowriter: 

START  ESP0D 

START  ESP0DDC  (if  selected) 

START  ESP0DEPH  (if  selected) 

START  ESP0D  (for  next  case) 


NEXT  FUNCTION 

When  NEXT  FUNCTION  is  typed  out,  the  orbit  determination  has 
been  completed  and  the  output  is  on  L0G  Tape  11.  Printed  output  is  in 
DS-0;  wrap  up  output  is  in  DS- 1 . 
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3.  2  SCHEDULE  MODE 


The  schedule  mode  requires  no  operator  intervention;  all  necessary 
inputs  for  system  operation  are  on  the  schedule  tape.  The  tape  setup  for 
the  schedule  mode  is  given  below. 


Table  3 - II.  Tape  Setup  for  Schedule  Mode 


Logical 

Tape  Number 

Setup 

Tape  Description 

Sc  ratch 

Data  is  transferred  to  this  tape 

i 

Sy  stem 

B2  master  system  tape 

2 

Schedule 

Job  tape  (input) 

<3> 

SEAI 

Backup  tape  for  logical  tape  NO.  4 

4 

SEAI 

Master  SEAI  tape.  ESP0D  uses  only 
the  sensors  and  elements  from  this  tap* 

6 

SRADU 

Observation  tape 

© 

Sc  ratch 

ESP$D  writes  blocks  of  C(#MM0N 
DATA  and  observations  on  this  tape 

<3> 

Planetary  ephemeris  tape 

O 

Scratch 

Special  binary  ephemeris  tape  (R,  V 
versus  Time)  generated  by  ESP$DEPH 
on  option 

© 

Output 

Off-line  output  tape 

<5> 

Prestored  Octal  change  cards 
(  Toggle  45,  ON) 

O 

Write  ring 

required 

O 

May  be  left 

out  if  proper  conditions  are  met 
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Toggle  Settings: 

Toggle  24,  ON 

Toggle  45,  ON,  if  octals  are  present  (Tape  15). 

The  input  deck  must  be  in  the  order  as  specified  in  Section  2.2.5 
(input  decks,  stacked  cases).  The  SCHEDULE  tape  is  interrogated  by  the 
B2  system  and  ESP0D  is  called  in  after  data  has  been  transferred  to  TapeO. 
If  it  is  a  cold-start,  then  Tape  L0G  7  is  not  read,  but  is  written  before  call¬ 
ing  the  next  section  of  the  orbit  determination  package.  If  it  is  a  condi¬ 
tioned  start,  ESP$D  reads  Tape  L$G  7,  expecting  to  find  C$MM$N  storage 
already  recorded  there. 

DS-0  prints  all  normal  ESP$DEPH  output.  DS-  1  prints  card  image 
of  card  prepared  by  ESP(#DEPH  for  GIPAR,  and  DS-2  punches  all  cards. 

3.  3  COLD  AND  CONDITIONED  START 

Column  30  of  the  JDC  card  indicates  whether  the  current  run  is 
started  from  a  cold  start,  a  conditioned  start,  or  a  conditional  start.  The 
distinction  between  the  cold  start  and  either  the  conditioned  or  conditional 
start  is  that  with  the  latter  cases  the  observations  are  recovered  from  a 
L($G  7  tape  which  was  generated  on  a  previous  cold- start  run.  This  elimi¬ 
nates  the  necessity  of  reprocessing  observation  cards  when  the  only  change 
to  the  input  deck  is  in  the  JDC  options  or  the  preliminary  data.  The  condi¬ 
tional  start  is  a  suboption  of  a  conditioned  start. 

3.  3.  1  Cold  Start 

The  following  procedure  describes  the  method  of  input  processing 
from  a  cold  start.  First,  the  JDC  card  is  read  and  interpreted  and  fol¬ 
lowed  by  the  preliminary  data.  The  program  then  interrogates  the  JDC 
card  to  determine  if  it  is  an  ESP$DEPH  run  only.  If  this  is  so,  ESPQiDEPH 
is  called  into  core  from  the  B2  Master  Tape  and  control  is  transferred  to 
it.  However,  if  ESP(jDDC  is  to  be  executed,  the  observation  data  is  read 
(from  the  SRADU  tape  or  observation  cards),  time  sorted  and  written  on  the 
data  tape.  Sensor  data  (from  the  SEAI  tape  or  on  cards),  if  any,  is  proc¬ 
essed  and  a  tape  is  generated  of  the  data  in  C0MM$N  storage.  ESPQiDDC 
is  then  called  from  the  B2  Master  Tape  and  control  is  given  to  it.  After 
execution  of  the  differential  corrector  package,  the  JDC  card  is  again 

interrogated  and  ESPCjDEPH  is  executed,  if  called. 
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3.  3.  2  Conditioned  Start 


The  conditioned  start  options  (column  30  of  JDC  card  =  1  or  2)  permit 
the  running  conditions  to  be  obtained  from  a  tape  record  generated  in  the 
last  iteration  of  some  preceding  differential  correction.  The  tape  left  on 
Tape  L(jG  7  after  a  differential  correction  contains  the  data  necessary  for 
this  type  of  continuation.  The  uses  of  the  conditioned  start  are: 

a)  Restarting  an  iteration  after  machine  failure  or  an  inter¬ 
ruption,  without  reprocessing  the  preliminary  data,  the 
observations,  and  sensor  data 

b)  Extending  the  previous  differential  correction  with  more 
iterations,  with  modifications  to  the  solution  vector  or  to 
the  force  models 

c)  Performing  a  prediction  with  the  results  of  a  previous 
differential  correction 

As  indicated  above,  it  is  possible  to  modify  some  of  the  preliminary 
data  associated  with  the  previous  differential  correction.  In  a  conditioned 
start,  a  new  JDC  card,  which  calls  for  a  conditioned  start  (column  30)  with 
new  preliminary  data,  if  any,  is  input.  The  foregoing  restrictions  apply: 

a)  The  object  must  remain  the  same 

b)  The  initial  conditions  may  be  changed  but  this  nullifies  any 
benefit  from  a  previous  differential  correction  run 

c)  The  solution  vector  may  be  changed;  if  drag  is  added,  and  a 
drag  model  was  not  previously  included,  it  must  be  provided 
with  the  new  preliminary  data 

d)  The  force  model  may  be  changed 

e)  New  observations  may  not  be  added 

f)  New  sensor  data  may  not  be  added 

g)  Sensor  type  data  and  standard  deviations  may  be  modified 
(STYPE  and  SIGMA  cards) 

3.3.3  Conditional  Start 

The  conditional  start  is  a  suboption  of  the  conditioned  start.  It  is 
provided  principally  to  prevent  running  a  prediction  using  initial  conditions 
which,  because  of  the  circumstances  of  their  formation,  may  be  worse  than 
the  best  available.  In  the  conditional  start  mode,  the  run  will  not  proceed 
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if  two  conditions  were  present  when  the  correction  terminated:  a)  converg¬ 
ence  was  not  obtained  on  the  final  iteration;  and  b)  the  next  to  last  iteration 
was  a  diverging  step.  ESP0DEPH  will  not  run  on  a  conditional  start  if  it 
follows  immediately  after  an  ESP$DDC  run  which  was  not  converging.  In 
this  case,  no  prediction  will  be  performed. 

3.4  ERROR  RETURNS  AND  RECOVERY 

Section  4.  1  of  this  document  lists  the  on-line  messages  and  Sec¬ 
tion  4.  2.  9,  the  off-line  messages.  When  an  incorrect  card  is  input,  the 
program  continues  to  read  until  it  finds  the  next  JDC  card,  or  exits  if  it 
is  the  last  case.  An  exception  to  this  rule  occurs  when  there  are  no  ob¬ 
servations  on  the  SRADU  Tape  for  the  vehicle  in  question.  The  program 
then  exercises  the  STOP-GO  option  until  the  operator  types  STOP  or  GO 
as  instructed  on  the  Flexowriter.  The  same  rule  applies  when  the  SEAI 
Tape  does  not  contain  the  initial  conditions  for  the  vehicle. 

3.5  INTERFACE  CONSIDERATIONS 
3.  5.  1  Core  Interface 

ESPOD  stores  the  predicted  Cartesian  position  coordinates  (X,  Y 
Z  in  Earth  radii)  and  the  associated  time  in  days  and  fractions  of  days  for 
the  times  specified  on  DAC  cards.  These  parameters  are  stored  in  core, 
beginning  in  location  15100.  These  parameters  are  not  left  in  core  if  the 
ESPQftDEPH  times  are  specified  on  DELTT  or  PRDCT  cards  (up  to  four 
sets  of  data  can  be  stored). 

3.5.2  Special  Binary  Ephemeris  Tape 

ESPOD  generates  an  ephemeris  tape  (Tape  L(#G  10)  of  X,  Y,  Z,  X, 

Y,  Z  versus  time  on  option  (column  53  of  JDC  card).  This  tape  may  be 
used  as  input  to  other  system  programs  (for  use  with  XYZLA).  The  Bi¬ 
nary  Ephemeris  Tape  format  is  described  in  Section  3.  7.  5  of  ESD-TDR- 
64-395  "ESPOD  Mathematical  and  Subroutine  Description." 
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4.  OUTPUT  DESCRIPTION 


ESPOD  produces  output  in  the  form  of  on-line  messages,  punched 
cards,  magnetic  tape  and  off-line  printout.  This  section  explains  the  for¬ 
mat,  the  parameters  produced,  and  in  the  case  of  optional  output,  how  to 
obtain  it.  Each  kind  of  output  will  be  discussed  in  turn,  and  samples  are 
given  to  clarify  the  description. 

4.  1  ON-LINE  MESSAGES 

On-line  messages  may  be  typed  by  the  Flexowriter  during  an  ESPOD 
run.  They  are  error  diagnostics  which  direct  the  operator  to  the  source 
of  a  program  stoppage. 

The  messages  below  are  arranged  alphabetically  to  aid  the  operator 
in  using  this  document. 

CONDITIONAL  REQUIREMENTS  NOT  MET  -  PROCEEDING  TO  NEXT  CASE 

This  message  appears  when  CLDSTR  =  2  and  CONVR  ^  0.  The 
program  then  finds  the  next  JDC  card  and  starts  to  read  more  data, 
or  exits  if  this  is  the  last  case. 

DATA  NAME  (  )  NOT  FOUND 

This  message  indicates  that  an  illegal  card  name  was  used,  and  the 
program  proceeds  reading  cards 

ERROR .  NO  OBS  ON  SRADU  FOR  SATELLITE  NO. 

TYPE  GO  TO  REREAD  TAPE,  STOP  FOR  NEXT  CASE. 

When  this  message  is  typed,  the  program  exercises  the  STOP-GO 
option  until  the  operator  types  GO  or  STOP  as  instructed  on  the 
Flexowriter. 

OBS  OVERFLOW  CORE,  NOT  IN  REVERSE  SORT,  ERROR. 

When  this  message  is  typed,  the  program  finds  the  next  JDC  card 
and  starts  to  read  more  data,  or  exits  if  this  is  the  last  case. 

SENSOR  DATA  OVERFLOWS  COMMON,  ERROR 

When  this  message  is  typed,  the  program  finds  the  next  JDC  card 
and  starts  to  read  more  data,  or  exits  if  this  is  the  last  case. 

TAPE  04  BAD- --MOUNT  BACKUP 

ELEMENTS  FOR  SATT.  NOT  ON  SEAI. 

TYPE- --GO  TO  REREAD  SEAI,  STOP  FOR  NEXT  CASE. 

When  this  message  is  typed,  the  program  exercises  the  STOP-GO 
option  until  the  operator  types  GO  or  STOP  as  instructed  on  the 
Flexowriter . 
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TAPE  04  BAD - MOUNT  BACKUP 

TYPE - GO  RETRY  TAPE,  STOP  NEXT  CASE. 

This  message  indicates  that  the  program  will  exercise  the  STOP- 
GO  option  while  the  backup  tape  is  being  mounted,  and  will  do  so 
until  the  operator  types  GO  or  STOP  as  instructed  on  the  Flexowriter. 

TAPE  ON  LOG  7  NOT  CORRECT.  I.  D.  IS 

TYPE-- -GO  TO  RETRY  TAPE,  STOP  FOR  NEXT  CASE. 

This  message  indicates  that  the  program  will  exercise  the  STOP-GO 
option  until  the  operator  changes  the  tape  and/or  types  GO  or  STOP  as 
instructed  on  the  Flexowriter. 

4.  Z  STANDARD  OUTPUT 

ESPOD's  standard  output  is  an  off-line  printout.  This  information 
contained  in  the  run  is  indexed  in  Figure  4-1,  which  directs  attention  to 
more  detailed  descriptions  of  particular  output  sections. 

4.  Z.  1  Input  Listing 

This  page  lists  the  input  cards  which  were  entered  for  the  run.  The 
columns  in  the  listing  correspond  exactly  with  the  columns  punched  on  the 
input  cards.  See  Section  Z.O  for  the  input  card  formats. 

These  data  are  a  quick  and  convenient  check  for  proper  input  to  the 
run.  Figure  4-Z  shows  a  sample  input  listing  page. 

4.  Z.  Z  Run  Header  Data 

The  Job  Description  Card  (JDC),  run  identification  data,  program 
control  and  reference  data,  and  initial  conditions  comprise  the  run  header 
data.  See  Section  Z.O  for  input  quantity  descriptions. 

These  data  identify  the  run,  labeling  the  various  input  quantities  and 
presenting  standard  coordinate  initial  conditions.  Figure  4-3  shows  a 
sample  header  data  page.  Each  line  is  described  as  follows: 

Page  Print 

Line  No.  Line  Description 

1  An  identification  obtained  from  JDC  columns  Z0-Z9;  run 

number,  programmer’s  name,  etc. 

Z  Vehicle  number  from  JDC  columns  4-7,  vehicle  name 

from  columns  8-19. 
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SAMPLE  RUN  DATA 


NO.  PAGES 


qy 

§■■■>.._ 

lltflaa..- 

isnisf 


INPUT  LISTING 

RUN  HEADER  DATA 

OBSERVATION  TYPE 

SENSOR  LOCATIONS 

RESIDUALS  PRINT 

ESTIMATES  OF  MEAN 

AND  STANDARD  DEVIATIONS 

CURVE  FIT  INTERATION 
SUMMARY 

TRAJECTORY  PRINT 


1 


1 


1  OR  MORE 


1 


1  OR  MORE 


1 


1 


1  OR  MORE 


Figure  4-1.  ESPOD  Standard  Run  Output 


SECTION 

4.2.1 

4.2.2 

4.2.3 

4.2.4 

4.2.5 

4.2.6 

4.2.7 

4.2.8 
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Figure  4-2.  Sample  ESPOD  Input  Listing  Page 
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Figure  4-3.  Sample  ESPOD  Header  Data  Page 


Page  Print 

Line  No.  Line  Description 


3 

COLD  START  or  NON-COLD  START:  By  flag  in  JDC 
column  30;  0  =  cold  start,  1  =  non-cold  start. 

4 

ALPHA  G  ZERO:  The  right  ascension  of  the  Greenwich 
meridian  at  midnight  on  the  day  of  epoch. 

5 

Numerical  value  of  Alpha  G  Zero,  in  degrees. 

6,  7,  8 

Epoch  time  of  the  run. 

9 

INITIAL  CONDITIONS:  First  estimate  of  conditions 
at  epoch. 

10,  11 

Symbols  and  numerical  values  of  initial  conditions  in 
geocentric  inertial  Cartesian  coordinates,  units  are 
kilometers  and  kilometer s  / second. 

12,  13 

Symbols  and  numerical  values  of  initial  conditions  in 
geocentric  inertial  spherical  polar  coordinates,  units 
are  degrees,  kilometers,  and  kilometer s  / second. 

The  following  are  printed  only  if  a  drag  parameter  value  has  been  input: 


14 

CDA/2M:  atmospheric  drag  parameter. 

15 

Numerical  value  input  for  C^A/2m,  units  are 
meter  s  2 /kilogram. 

4.  2.  3  Observation  Type 

This  page  lists  the  time-ordered  observations  input  for  this  run. 
Figure  4-4  shows  a  sample  observation  type  page.  The  column  symbols 
and  their  descriptions  are: 

Column 

Symbol  Description 


ID 

Observing  station!s  identification  number,  corresponds 
to  a  number  under  SENSOR  LOCATIONS. 

T-TO 

Time  of  the  observation  in  minutes  from  the  epoch  (TO). 

>> 

YR 

MN 

DAY 

HR 

MIN 

SEC 

Greenwich  time  at  which  the  observation  was  obtained: 
Year,  month,  day,  hour,  minute,  and  decimal  seconds. 

j 
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Figure  4-4.  Sample  ESPOD  Observation  Type  Pag 


Column 

Symbol  Description 

R  Range  observation,  in  kilometers. 


A 

Azimuth  observation,  in 
from  true  north. 

degrees  positive  clockwise 

EL 

Elevation  observation,  in  degrees  from  the  horizon. 

R  DOT 

Range  rate  observation, 

in  kilometers/second. 

HA 

Hour  angle  observation, 

in  degrees. 

D 

Declination  observation, 

in  degrees. 

The  numerical  values  for  the  observations  are  listed  in  the  above 
columns  on  one  or  more  pages.  All  quantities  are  positive  unless  prefixed 
by  a  minus  sign.  Columns  for  which  the  particular  sensor  does  not  report 
data  automatically  contain  zeroes. 


4.  2.  4  Sensor  Locations 

This  page  shows  the  locations  of  the  sensors  for  which  observations 
have  been  input  for  this  particular  run.  Figure  4-5  shows  a  sample  page 
print  of  the  sensor  locations.  The  following  listing  describes  the  lines  on 
this  page: 


Page  Print 
Line  No. 


3  to  Page 
End 


Line  Description 

SENSOR  LOCATIONS:  The  identification  and  geodetic 
location  of  stations  furnishing  observations,  referenced 
to  the  WGS  I960  spheroid. 

Column  identifications 

ST,  NO.  Observing  station’s  identification  number 

LAT.  North  latitude  of  observing  station,  in  degrees 

LONG.  West  longitude  of  observing  station,  in  degrees 

ALT.  Altitude  of  observing  station  above  mean  sea 

level,  in  meters 

Numerical  values  in  the  columns  indicated  in  line  17. 


Zeroes  if  the  station  reports  angular  observations  in  azimuth  and  elevation, 
and  vice  versa. 
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igure  4-5.  Sample  ESPOD  Sensor  Eocations  Page 


4.  2.  5  Residuals  Print 


The  residuals  printed  are  defined  as  the  difference  between  the 
measured  observations  and  the  computed  observations:  measured  -  com¬ 
puted  =  residual.  The  computed  observations  are  based  on  the  initial  condi 
tions  before  the  first  iteration  and  on  the  revised  initial  conditions  after 
each  successive  iteration.  The  residuals  are  time-ordered  like  the  obser¬ 
vation  type  listing  just  described,  and  correspond  exactly  in  time  with  their 
associated  measured  observations.  Figure  4-6  shows  a  sample  residuals 
print  page.  The  column  symbols  and  their  descriptions  are: 


Column 

Symbol 

Description 

ID 

Observing  station ’ s  identification  number;  corresponds 
to  a  number  under  SENSOR  LOCATIONS. 

DATE 

Y  ^ 

MM 

DD 

HH  ' 

MM 

SS.  SSj 

Time  associated  with  the  measured  observation 
differenced  to  obtain  this  residual.  Symbol  indicates 
allowable  number  of  digits  which  can  be  shown,  i.  e., 

Y  =  year  and  last  digit  only.  Successively,  others  are 
month,  day,  hour,  minute,  and  decimal  seconds. 

N 

Serial  number  assigned  to  each  residual  for  identification 
purposes.  It  is  constant  through  the  run. 

RANGE 

KM 

Range  residual  in  kilometers 

AZ,  HA 

DEG 

Azimuth  or  hour  angle  residual,  whichever  type  of 
observation  this  sensor  measures,  in  degrees. 

EL,  DC 

DEG 

Elevation  or  declination  residual,  whichever  type  of 
observation  this  sensor  measures,  in  degrees. 

RDOT 

KM/SEC 

Range  rate  residual,  in  kilometers /second. 

U 

KM 

Up  component  of  the  position  residual,  collinear  with  and 
positive  in  the  same  direction  as  the  radius  vector,  in 
kilometers . 

V 

KM 

Down  component  of  the  position  residual,  orthogonal  to 
the  radius  vector,  in  the  orbit  plane,  positive  in  the 
direction  of  motion,  in  kilometers. 

W 

KM 

Cross  component  of  the  position  residual,  normal  to  the 
orbit  plane  and  positive  in  the  direction  of  the  angular 
momentum  vector  to  complete  a  right-handed  coordinate 
system,  in  kilometers. 
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igure  4-6.  Sample  ESPOD  Residuals  Page 


Column 

Symbol  Description 


VM 

KM 

Magnitude  of  the  position  residual  or  displacement 
vector,  in  kilometers. 

DEL  T 
MIN 

Time  separation  of  computed  and  measured  positions, 
assuming  Keplerian  mean  motion,  in  minutes. 

U 

DEG 

The  argument  of  latitude  (angle  from  the  ascending  mode 
to  the  radius  vector  in  the  orbit  plane)  at  this  time,  in 
degrees. 

BETA 

DEG 

Residual  angle  between  the  measured  position  vector 
and  the  computed  orbit  plane. 

The  numerical  values  for  the  residuals  are  listed  in  the  above  col¬ 
umns  on  one  or  more  pages.  A  symbol  or  letter  following  a  numerical 
value  indicates  that  residual  has  been  deleted  from  further  calculation  on 
a  given  criterion,  which  is  one  of  the  following: 


Symbol 

Criterion 

sjc 

Observation  deleted  by  a  DELETE  input  (removed  by 
analyst’s  choice). 

G 

Observation  deleted  as  a  gross  outlier.  See 

Section  2.2. 3. 5  for  further  details. 

K 

Observation  deleted  from  this  iteration  by  KRMS  test. 

The  nominal  value  for  K  may  be  changed  as  described 
in  Section  2.2.3.  11. 

If  observations  are  input  from  a  station  for  which  the  location  has  not 
been  input,  a  note  will  appear  in  place  of  the  residuals  for  this  sensor  as 
follows  : 

STATION  (ST.  ID.  )  NOT  IN  MASTER  SENSOR  LIST 

Note  that  the  program  cannot  compute  a  position  vector  without  know¬ 
ing  the  station  location,  and  hence  cannot  calculate  the  residual. 
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4.  2.  6  Estimates  of  Mean  and  Standard  Deviations 


To  aid  the  analyst  in  determining  where  systematic  error  solution 
is  required,  or  in  deleting  data,  estimates  of  mean  and  standard  devia¬ 
tions  by  sensor  and  type  are  printed.  Figure  4-7  shows  a  sample  page 
containing  these  data.  A  given  station's  data  are  arranged  into  columns 
of  three  lines  each.  The  column  symbols  and  their  descriptions  are: 

Column 

Symbol  Description 


ST.  ID 

Observing  station's  identification  number,  corresponds  to 
a  number  under  SENSOR  LOCATIONS. 

R 

Numerical  values  associated  with  range,  in  kilometers. 

A, HA 

Numerical  values  associated  with  azimuth  or  hour  angle, 
in  degrees . 

E,  D 

Numerical  values  associated  with  elevation  or  declination, 
in  degrees . 

RDOT 

Numerical  values  associated  with  range  rate,  in  kilometers 
per  second. 

The  last  column  on  the  right  identifies  the  data  on  each  line  as 


follows : 

MEAN 

The  arithmetic  mean  or  average  values  of  the  residuals 
for  each  type  of  data  observed  by  the  given  station. 

ESTD 

The  estimated  standard  deviation  (one  sigma)  for  each  type 
of  data  observed  by  the  given  station.* 

NA/NR 

The  number  of  observations  accepted/the  number  of  obser¬ 
vations  rejected,  for  each  type  of  data  observed  by  the  given 
station. 

Note:  The  estimated  standard  deviation  is  calculated  about  the  computed 
mean  value  and  not  about  an  assumed  mean  of  zero. 
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Figure  4-7.  Sample  ESPOD  Estimates  of  Mean  and  Standard  Deviations  Page 


4.2.7  Curve  Fit  Iteration  Summary 


This  page  shows  the  results  and  convergence  status  for  a  given 
iteration.  If  the  solution  is  converging,  correlation  and  covariance  lower 
triangular  matrices  are  printed.  Figure  4-8  shows  a  sample  curve  fit 
summary  page.  The  first  line  on  the  page  identifies  the  iteration  number 
in  this  particular  run.  The  following  describes  the  tabulated  solution 
vector  data,  arranged  in  columns  as  shown  below: 


Column 

Symbol  Description 

CATEGORY  1  These  are  the  orbital  and  drag  parameters,  identified  by 
VARIABLES  a  number  and  a  name  as  follows: 


1  ALPHA 

2  DELTA 

3  BETA 

4  AZ 

5  R 


Right  ascension  (degrees) 

Declination  (degrees) 

Flight  path  angle  from  local  vertical 
(degrees) 

Azimuth  to  inertial  velocity  vector, 
clockwise  from  true  north  (degrees) 

Radius  vector  from  geocenter 
(kilometers) 


6  V 

7  CDA/2M 


Velocity  vector  magnitude 
(kilometers/ sec) 

,  2 

Drag  parameter  (meters  per 
kilograms) 


8  K  Drag  variation  (secular  option: 

meters^  per  kilogram  per  day; 
periodic  option:  meters^  per 
kilogram) 


DELTA  The  corrections  applied  by  the  program  to  each  variable. 

OLD  Numerical  values  for  the  variables  from  the  previous 

ite  ration. 


NEW  Numerical  values  for  the  variables  from  this  iteration 

(NEW  =  OLD  +  DELTA). 

SIGMA  The  uncertainty  in  each  variable,  computed  from  the 

covariance  matrix. 

BOUNDS  The  constraints  applied  to  the  changes  which  the  program 

is  allowed  make  to  the  variables. 
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Figure  4-8.  Sample  ESPOD  Curve  Fit  Summary  Page 


If  Category  2  variables  (station  and  data  biases)  are  being  solved 
for,  a  data  tabulation  much  like  that  just  described  will  be  listed  next. 
Its  columns  are  as  follows: 


Col  um  n 

Symbol  Description 


CATEGORY  2  These  variables  will  be  numbered  from  the  next  digit 
VARIABLES  following  the  last  Category  1  variable  number  to  50. 

Each  number  will  be  followed  by  a  station  identification 
number,  as  used  in  SENSOR  LOCATIONS,  and  the  name 
of  the  variable  which  can  be  any  of  the  following: 


R 

Range  (kilometers) 

A,  HA 

Azimuth  or  hour  angle  (degrees) 

E,  D 

Elevation  or  declination  (degrees) 

R  DOT 

Range  rate  (kilometers  per  second) 

LAT 

Station  north  latitude  (degrees) 

LONG 

Station  east  longitude  (degrees) 

ALT 

Station  altitude  (meters) 

T 

Time  (seconds) 

DELTA  1 

OLD 

NEW 

SIGMA 
BOUNDS  j 

[  s 

>  Same  as 

for  Category  1  variables. 

The  following  line  states  "SOLUTION  IS  (IS  NOT)  CONVERGING.  " 
The  solution  is  converging  if  the  current  RMS  of  the  residuals  is  smaller 
than  the  best  RMS  to  this  point  in  the  run. 


This  message  is  followed  immediately  by  "SOLUTION  IS  (IS  NOT) 

AFFECTED  BY  BOUNDS.  "  The  program  first  forms  an  unbounded 

f  DELTA’  ^2 

solution,  and  if  it  satisfies  the  condition  ^(b'qujsjdS-  )  ^  "SOLUTION 

IS  NOT  AFFECTED  BY  BOUNDS"  is  printed.  If  the  condition  is  not 
satisfied,  the  system  of  equations  is  solved  until  the  condition  is  satisfied, 
and  "SOLUTION  IS  AFFECTED  BY  BOUNDS"  is  printed. 
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Next  are  listed  the  current  RMS  of  the  residuals,  the  RMS  predicted 
for  the  next  iteration  (based  on  the  DELTA'S  and  (A^A)  ),  and  the  best 

RMS  so  far  in  the  curve  fit  run. 

A  correlation  matrix  in  lower  triangular  form  is  printed  next  if 
the  solution  is  converging.  Its  rows  and  columns  are  numbered  to  cor¬ 
respond  to  the  Category  1  variables  as  they  are  numbered  in  the  solution 
vector  data  tabulation  previously  described. 

T  “ 1 

Following  the  correlation  matrix  is  the  (AXA)  matrix  or  covariance 

matrix,  also  in  lower  triangular  form.  It  is  printed  only  if  the  solution 

is  converging.  If  it  were  printed  and  the  solution  were  not  converging, 

the  information  in  the  matrix  would  be  false,  because  the  solution  would 

T  - 1 

be  based  on  the  previous  (A  A)  in  the  last  converging  iteration  and 
—  ^^UNDS)  ^  the  most  rightly  constrained  solution  which  the  program 
can  produce. 

The  next  item  will  appear  on  the  final  iteration  only,  and  will  be  one 
of  the  following: 

SOLUTION  HAS  CONVERGED 

MAXIMUM  ITERATIONS  EXCEEDED 

BOUNDS/8  FAILED 

This  indicates  the  criterion  on  which  the  run  has  been  terminated. 

4,2.8  Trajectory  Printout 

The  last  page(s)  of  the  run  consist  of  blocks  of  trajectory  and  re¬ 
lated  data,  updated  to  times  from  the  epoch  which  were  specified  by  input. 
Figure  4-9  shows  a  sample  trajectory  printout  page.  Each  item  in  the 
block  is  preceded  by  an  identifying  symbol.  The  first  line  of  each  block 
contains  the  date  and  Greenwich  mean  time  for  the  data,  time  in  minutes 
from  epoch,  and  time  in  days  from  January  0  of  the  year  of  the  epoch. 

The  following  list  describes  the  parameters  associated  with  each  identi- 
fying  symbol: 
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Figure  4-9.  Sample  ESPOD  Trajectory  Printout  Page 


Symbol  Description 


x  ) 

Y 

? 

X  ' 

Y 

z 

Components  of  the  position  and  velocity  vector  in  geo¬ 
centric  inertial  Cartesian  coordinates.  It  is  a  right- 
handed  orthogonal  system  where  the  X  axis  is  in  the 
direction  of  the  vernal  equinox  and  the  Z  axis  is  in  the  di¬ 
rection  ot  true  north.  Units  are  kilometers  and  kilometers/ 
second.  Coordinates  are  true  of  ol'O  day  of  epoch. 

ALFA 

Right  ascension,  in  degrees 

DLTA 

Declination,  in  degrees 

BETA 

Flight  path  angle,  in  degrees,  positive  downward  from 
the  local  vertical. 

AZ 

Azimuth  of  the  velocity  vector,  in  degrees,  positive 
clockwise  from  true  north. 

R 

Radius  vector  from  the  geocenter,  in  kilometers. 

V 

Magnitude  of  the  velocity  vector,  in  kilometers  per 
second. 

ALT 

Height  of  the  satellite  above  mean  sea  level,  in  kilo¬ 
meters. 

LAT 

Geodetic  north  latitude  of  the  satellite,  in  degrees. 

LONE 

East  longitude  of  the  satellite,  in  degrees. 

SMA 

Semimajor  axis,  in  kilometers. 

ECC 

Eccentricity  of  the  orbit. 

INC 

Inclination  of  the  orbit  plane  to  the  equator,  in  degrees, 
positive  counterclockwise  from  the  equatorial  to  the  orbit 
plane  at  the  ascending  node. 

NODE 

Right  ascension  of  the  ascending  node,  in  degrees. 

OMG 

Argument  of  perigee,  in  degrees,  positive  in  the  direction 
of  motion  from  the  ascending  node. 

M 

Mean  anomaly,  in  degrees,  positive  in  the  direction  of 
motion  from  perigee. 

UX  \ 

UY  ) 

UZ 

Direction  cosines  of  the  position  in  Cartesian  coordinates, 
with  axes  directed  as  in  the  XYZ  system  previously 
described. 
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Symbol 

RPVX  I 
RPVY  ) 
RPVZ  J 

ALAT 


TAU 


PRD 

1/A 


D 

APOG 


Description 

Components  in  Cartesian  coordinates  of  a  vector  in  the 
orbit  plane  which  is  orthogonal  to  the  position  (r)  and 
angular  momentum  (h)  vectors. 

Argument  of  latitude,  in  degrees,  equals  the  sum  of  the 
argument  of  perigee  and  the  true  anomaly. 

Time  until  the  next  ascending  nodal  crossing,  in  minutes 
from  epoch. 


Osculating  period  of  the  orbit, 

Inverse  of  the  semimajor  axis 
(indeterminacy-free  element). 

Indeterminacy-free  element  = 


in  minute  s . 
in  earth  radii 


R  •  R 


in  (earth  radii) 


1/2 


Altitude  of  apogee  above  a  mean  equator,  in  kilometers. 


PRG 


Altitude  of  perigee  above  a  mean  equator,  in  kilometers. 


ELLIPSE 

HYPERBOLA 


Prints  one  or  the  other  to  describe  the  orbit's  conic  form. 
If  HYPERBOLA,  many  of  the  preceding  values  will  be 
omitted. 


XVM 
YVM 
ZVM  l 
XDVM  f 
YDVM 
ZDVM  t 


Vehicle  position  and  velocity  relative  to  the  moon,  in  a 
coordinate  system  defined  like  the  geocentric  inertial 
Cartesian  system  in  this  section,  but  translated  to  the 
selenocenter .  These  coordinates  are  true  of  0.0  day  of 
epoch. 


XVS 

YVS 

ZVS 

XDVS 

YD  VS 

ZDVS 


Vehicle  position  and  velocity  relative  to  the  sun,  in  a 
coordinate  system  defined  like  the  geocentric  inertial 
Cartesian  system  in  this  section,  but  translated  to  the 
heliocenter.  These  coordinates  are  true  of  0^0  day  of 
epoch. 


DVFM 


Distance  from  the  vehicle  to  the  selenocenter,  in 
kilometers . 


DVFS 


Distance  from  the  vehicle  to  the  heliocenter,  in  kilometers 


The  above  parameters  are  repeated  for  as  many  updates  as  were  re¬ 
quested  in  the  input.  These  data  complete  the  standard  output  from  ESPOD 
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4.2.9  Off-Line  Messages 


Off-line  messages  are  given  in  the  printed  output  to  aid  the  analyst 
in  diagnosing  stoppages  which  occurred  during  a  run.  In  many  cases,  the 
program  will  bypass  the  case  in  which  an  error  halt  occurred  and  go  on  to 
the  next  one.  When  the  printed  output  is  reviewed,  the  program  will  have 
indicated  the  case  in  which  the  error  halt  occurred  with  an  off-line  me ssage . 

To  aid  the  operator  in  quickly  finding  the  explanations  for  these 
messages,  they  are  arranged  alphabetically  in  the  following  list. 

COLUMN  68  in  CARD  7  IS  ILLEGAL.  CARD  IMAGE  BELOW 

(Followed  by  the  data  actually  punched  on  card  7) 

This  message  is  printed  when  column  68  in  element  card  7  is  not 
a  0,  1,  or  2.  The  program  then  finds  the  next  JDC  card  and 
continues  reading  input.  If  this  is  the  last  case,  the  program  exits. 

CONDITIONAL  REQUIREMENTS  NOT  MET - PROCEEDING  TO  NEXT  CASE. 

This  message  indicates  that  the  cold  start  flag  CLDSTR=2  and 
CONVR  /  0.  The  program  finds  the  next  JDC  card  and  continues 
reading  input.  If  this  is  the  last  case,  the  program  exits. 

DATA  NAME)*  •  •)  NOT  FOUND. 

An  illegal  card  name  has  been  used  in  place  of  ( •  •  •)•  The  program 
finds  the  next  JDC  card  and  continues  reading  input,  or  exits,  if 
this  is  the  last  case. 

ELEMENTS  FOR  SATT. _ NOT  ON  SEAL 

This  message  is  printed  off-line  for  the  record  and  on-line 
for  action.  The  accompanying  on-line  message  is: 

Tape  04  bad  -  mount  backup 

Elements  for  satt.  not  in  SEAI 

Type  —  GO  to  reread  SEAI,  STOP  for  next  case. 

❖  *  ^END  OF  FILE  ENCOUNTERED  READING  THE  EPHEMERIS  TAPE 

This  message  records  the  reason  for  the  program’s  halt  on  this 
case.  It  then  finds  the  next  JDC  card,  reads  the  input  and 
continues  or  exits  if  this  is  the  last  case. 
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*  *  *  EPHEMERIS  TAPE  ARGUMENT  TOO  LARGE  •  •  • 

T  =  . _ SECONDS  FROM  1950.0 

This  message  indicates  that  the  time  for  which  data  are  requested 
is  after  19  December  1969,  which  is  the  end  of  the  ephemeris 
tape's  data.  The  program  proceeds  to  the  next  case  or  exits  if 
this  is  the  last  case. 

*  *  *  EPHEMERIS  TAPE  ARGUMENT  TOO  SMALL  •  *  * 

T  =  .  SECONDS  FROM  1950.0 


This  message  indicates  that  the  time  for  which  data  are 
requested  is  before  28  August  I960,  which  is  the  beginning 
of  the  ephemeris  tape's  data.  The  program  proceeds  to  the 
next  case  or  exits  if  this  is  the  last  one. 

ERROR.  NO  OBS  ON  SRADU  FOR  SATELLITE  NO. _ 

This  message  records  the  reason  for  proceeding  to  the  next  case 
or  exiting  if  this  is  the  last  one.  The  accompanying  on-line 
comment  for  action  is: 

Error.  No  obs  on  SRADU  for  satellite  no. _ 

Type  GO  to  reread  tape,  STOP  for  next  case. 

ICOND"  ^ 

ICTYP  >  CARD  REQUIRED  FOR  THIS  RUN 
ITIME  J 

This  message  indicates  that  an  ICOND,  ICTYP,  or  ITIME  card 
was  missing  from  the  input  deck.  The  program  proceeds  to  the 
next  case  or  exits  if  this  is  the  last  one. 

ILLEGAL  TESSERAL  NM  REQUESTED,  THE  PROGRAM  IS  IGNORING 
IT  AND  PROCEEDING 

This  message  is  printed  if  the  value  of  a  requested  tesseral 
harmonic  is  greater  than  or  equal  to  the  degree.  The  program 
proceeds  with  the  case  as  though  the  illegal  input  were  never 
there . 

JDC  CARD  NOT  FOUND.  ID  WHICH  TERMINATED  RUN  IS  _ 

While  searching  the  input  for  a  JDC  card,  the  program  found  an 
identification  name  to  be  ENDAT,  ZZZZZ,  or  EEEEE.  The 
program  exits. 


One  of  these  will  appear  in  the  message. 
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*  *  ❖MAJOR  PROGRAM  ERROR  •  *  •  POSSIBLE  INPUT  AND/OR 
MACHINE  ERROR 

This  message  indicates  that  the  total  RMS  on  the  first  iteration 
has  exceeded  the  maximum  floating  point  number  which  the 
computer  can  handle.  The  program  returns  and  reads  input 
for  the  next  case  or  exits  if  this  is  the  last  one. 

NO  JDC  IN  400  CARDS 

This  message  records  the  fact  stated.  The  program  exits  when 
this  message  has  been  written  for  printing. 

OBS  OVERFLOW  CORE,  NOT  IN  REVERSE  SORT,  ERROR. 

This  message  records  the  facts  stated  to  indicate  why  the  program 
proceeded  to  the  next  case  or  exited.  Observations  may  be  added 
until  they  overflow  core  if  and  only  if  they  are  in  reverse  sort. 

This  message  is  typed  on-line  at  the  same  time  it  is  written  for 
off-line  printing. 

PARTIAL  (J)  WITH  RESPECT  (I)  ASKED  FOR 

This  message  indicates  that  a  non-orbital  parameter  partial 
derivative  has  been  requested  which  the  program  does  not 
recognize,  because  the  allowable  numbers  (J  =  6,  1=10)  have 
been  exceeded.  The  program  sets  this  illegal  partial  derivative 
equal  to  zero  and  proceeds. 

REQUIRED  SENSORS  TABLE  FULL.  ITEM  NOT  SAVED  IS _ ’ 

This  message  is  written  for  printing  when  the  stated  fact  occurs. 

SENSOR _ NOT  IN  MASTER  SENSOR  LIST 

This  message  is  printed  when  the  program  cannot  find  this  sensor 
identification  in  the  master  list,  and  so  cannot  calculate  a  residual 
because  it  does  not  know  the  sensor  location.  The  program  continues 
through  the  observations  until  a  different  sensor  identification  is 
found  which  is  the  master  list. 

SENSOR  DATA  OVERFLOWS  COMMON,  ERROR 

This  message  indicates  that  too  large  a  system  of  sensors  and  non¬ 
orbital  parameters  was  requested,  and  variable  storage  capacity 
has  been  exceeded.  The  program  reads  input  for  the  next  case  or 
exits  if  this  is  the  last  case. 


Will  be  a  sensor  identification  number.  If  the  table  overflows  by  more 
than  one  sensor,  a  separate  message  is  printed  for  each.  The  program 
continues  the  computation  as  though  the  overflowing  sensors  were  never 
input. 
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SUBERR  EXIT  OCCURRED,  SUBERR  J  LOC,  PROGRAM  JMP  LOC 
A  REGISTER,  Q  REGISTER 


This  message  records  that  the  exit  to  the  next  case  was  instigated 
by  the  ERR$R  subroutine.  The  contents  of  the  A,  Q,  and  JA 
registers  are  printed  to  aid  in  diagnosing  the  error.  SUBERR 
J  LOC  is  the  address  from  which  the  ERROR  subroutine  was 
called,  and  the  PROGRAM  JMP  LOC  is  the  subroutine  which 
called  the  error-incurring  subroutine. 

TAPE  ON  LOG  7  NOT  CORRECT.  I.D.  IS _ . 

This  message  records  the  reason  for  exit  to  the  next  case,  and 
is  written  for  off-line  printing  at  the  same  time  as  this  on-line 
mes  s age : 

Tape  on  LOG  7  not  correct.  I.D.  is 

Type  -  GO  to  retry  tape,  STOP  for  next  case. 

THE  FOLLOWING  CARD(S)  COULD  NOT  BE  CONVERTED  ERR 

LOCATION 

This  message  is  printed  when  subroutine  BCDC^BS  cannot  read 
an  observation  card  for  any  reason.  The  core  address  is  given 
in  the  message  where  the  subroutine  tried  to  read  the  observation 
card,  and  the  card  is  imaged  as  read  on  the  line  following  the 
message.  The  program  proceeds  as  though  the  illegal  card  were 
never  present. 

4.3  OPTIONAL  OUTPUT 

This  section  describes  optional  outputs  and  how  to  obtain  them. 
These  options  are  described  in  the  same  order  in  which  they  would  be 
present  on  the  printout.  Figure  4-10  shows  the  ordering  of  the  optional 
print  and  the  section  in  which  it  is  subsequently  described. 

4.3.1  Input  Listing 


Any  flags  input  to  produce  optional  output  will  be  present  in  the 
input  listing.  Since  it  is  simply  an  interpretation  of  the  punched  input 
cards  or  listing  from  the  magnetic  tape,  there  are  no  optional  outputs  for 
this  page  of  the  printout. 

4.3.2  Run  Header  Data 


There  are  no  optional  outputs  for  the  run  header  data. 
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SECTION 


SAMPLE  RUN 


DATA 


INPUT  LISTING 


HEADER  DATA 


PROGRAM  CONSTANTS* 


SENSOR  INFORMATION* 


OBSERVATION  TYPE 


SENSOR  LOCATIONS 


RESIDUALS  PRINT* 


ESTIMATES  OF  MEAN 
AND  STANDARD  DEVIATIONS 

CURVE  FIT  INTERATION 
SUMMARY 


1>._ _ 


TRAJECTORY  PRINT* 


♦OPTIONAL  PRINT  AVAILABLE 


NO.  PAGES 


1 


1 


1 


1 


1  OR  MORE 


1 


1  OR  MORE 


1 


1 


1  OR  MORE 


4.3.1 


4.3.2 


4.3.3 


4.3.4 


4.3.5 


4.3.6 


4.3.7 


4.3.8 


4.3.9 


4.3.10 


Figure  4-10.  ESPOD  Standard  and  Optional  Run  Output 
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4.  3.  3  Program  Constants 


The  program  constants  may  be  printed  by  placing  a  1  flag  in  JDC 
card  column  36.  These  constants  may  be  changed  by  making  an  appropriate 
input  on  a  99-card,  as  described  in  Section  5.0. 

For  a  listing  of  the  constants  shown  by  this  print  option,  see 
Section  5.0.  Figure  4-11  shows  a  sample  program  constants  print  page . 

4.3.4  Sensor  Information 


The  flag  which  causes  the  program  constants  to  be  printed  also  re¬ 
sults  in  the  sensor  information,  which  the  program  has  in  its  master  list, 
being  printed.  Figure  4-12  shows  a  sample  sensor  information  print  page. 
The  column  symbols  and  the  information  pertaining  to  each  sensor  are 
described  in  the  following  list: 


Column  Symbol 


SENSOR 

NO. 


SIGMA 

TYPE 


Description 

The  sensor  identification  number  which  is 
used  to  reference  a  given  sensor  location 
or  properties,  when  they  are  desired  for 
use  in  the  program. 

The  identification  for  the  group  of  standard 
deviations  associated  with  certain  stations1 
data. 


RANGE/RA(FR)  The  standard  deviation  in  range,  or,  if  the 

sensor  is  a  Baker-Nunn  camera,  the  standard 
deviation  in  right  ascension  when  it  is  a 
field- reduced  observation. 


AZ/DEC(FR) 


EL/RA(PR) 


RDT  /DEC(PR) 


GSUBS 


The  standard  deviation  in  azimuth,  or  in  the 
case  of  a  Baker-Nunn  camera  sensor,  the 
standard  deviation  in  declination  when  the 
observation  is  field  reduced. 

The  standard  deviation  in  elevation,  or  in 
the  case  of  a  Baker-Nunn  camera  sensor, 
the  standard  deviation  in  right  ascension 
when  the  observation  is  precision  reduced. 

The  standard  deviation  in  range  rate,  or  in 
the  case  of  a  Baker-Nunn  camera  sensor, 
the  standard  deviation  in  declination  when 
the  observation  is  precision  reduced. 

The  gross  outlier  editing  criterion  for  the 
observations  in  terms  of  N  sigma. 
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Figure  4-11.  Sample  ESPOD  Program  Constants  Page 
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Figure  4-1Z.  Sample  ESPOD  Sensor  Information  Page 


Column  Symbol 


Description 


REFR.  Applies  the  program's  refraction  correction 

FLAG  to  all  observations  from  a  given  sensor  if 

a  "1"  is  input.  If  a  "0"  is  input,  no  refraction 
correction  is  performed  by  the  program. 

NSUBS  The  local  index  of  refraction  for  the  given 

station  and  the  time  of  interest.  Used  in  the 
refraction  correction  routine. 


4.  3. 5  Observation  Type 

There  are  no  optional  outputs  for  the  observation  type,  since  it  is 
simply  a  time-ordered  listing  of  the  observations  input  to  the  program. 

4.3.6  Sensor  Locations 

There  are  no  print  options  for  the  sensor  locations.  It  is  simply  a 
listing  of  input  data  for  this  run  from  the  master  sensor  data  table  con¬ 
tained  in  the  program. 

4. 3.  7  Residuals  Print 

Two  options  are  available  on  the  residuals  print  which  replace  the 
columns  which  are  headed  with  U,  V  and  W  on  the  standard  output.  By 
setting  the  DCFLG(5)  =  1,  the  components  of  a  vector  describing  the  dis¬ 
placement  of  the  observed  with  respect  to  the  computed  position  are  printed 
for  each  observation,  in  a  normal,  tangential  and  cross  direction.  Fig¬ 
ure  4-13  shows  a  residuals  page  print  with  the  STW  option.  The  column 
symbols  are  described  as  follows: 

Column  Symbol  Description 

S  Vector  component  orthogonal  to  the  velocity 

KM  vector  and  in  the  orbit  plane,  forming  a 

right-handed  system  with  T  and  W. 

T  Vector  component  collinear  with  and  in  the 

KM  same  direction  as  the  velocity  vector. 

W  Vector  component  normal  to  the  orbit  plane 

KM  and  in  the  same  direction  as  the  angular 

momentum  vector. 
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Figure  4-13.  Sample  ESPOD  Residuals  (STW  Option)  Page 


By  setting  DCFLG(5)  =  2,  the  components  of  a  vector  describing  the 
displacement  of  the  observing  station  which  would  be  required  to  reduce 
the  residual  to  zero  are  computed.  Figure  4-14  shows  a  residuals  page 
print  with  this  LLH  option.  The  column  symbols  and  their  explanations 
are  as  follows : 


Column  Symbol 


Desc  r  iption 


ST.  L AT 
DEG 

ST.  LONE 
DEG 

ST.  HT 
KM 


Displacement  component  in  geocentric  latitude, 
in  degrees,  positive  northward. 

Displacement  component  in  longitude,  in 
degrees,  positive  eastward. 

Displacement  component  in  height,  in  kilometers 
above  mean  equatorial  sea  level. 


4.  3.  8  Estimates  of  Means  and  Standard  Deviations 


There  are  no  output  print  options  for  the  page  containing  the  estimates 
of  mean  and  standard  deviations. 


4.  3.  9  Curve  Fit  Iteration  Summary 

There  are  no  output  print  options  for  the  page  containing  the  curve 
fit  iteration  summary. 

4.  3.  10  Trajectory  Printout 

In  addition  to  the  block  of  data  printed  for  each  trajectory  update 
time,  as  described  in  Section  4.  2.  7,  updated  error  matrices  are  available 
in  a  variety  of  coordinate  systems.  The  variations  are: 

Orbit  Plane  Coordinates  (UVW) 

Spherical  Polar  Coordinates  (ADBARV) 

Cartesian  Coordinates  (XYZ) 

This  print  option  is  obtained  by  placing  a  11 1  11  flag  in  column  55  of 
the  JDC  card.  Figure  4-15  shows  a  sample  page  print  of  the  updated 
error  matrices  option. 

The  optional  updated  matrices  are  Msigma  and  rho"  symmetric 
matrices  given  in  lower  triangular  form.  The  diagonal  terms  are  the 
estimated  standard  deviation  (a)  and  the  off  diagonal  terms  are  the  corre¬ 
lation  coefficients  (p).  The  following  list  describes  the  parameters  given 
in  the  columns  of  each  optional  updated  matrix: 
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Figure  4-14.  Sample  ESPOD  Residuals  (LLH  Option)  Page 


4-34 


Figure  4-15.  Sample  ESPOD  Trajectory  (Updated  Error  Matrices  Option)  Printout  Page 


ORBIT  PLANE  SIGMA  AND  RHO  MATRIX 


Column  Symbol 

U 

1 


V 

2 


W 

3 


UDOT 

4 

VDOT 

5 

WDOT 

6 

CDA/2M 

7 

or 

K 

7 


Description 

Position  vector  (up)  component  in  the  direction 
of  the  radius  vector  and  collinear  with  it,  in 
kilometer  s. 

Position  vector  (down)  component  in  the 
direction  of  motion  and  orthogonal  to  the  radius 
vector  in  the  orbit  plane,  in  kilometers. 

Position  Vector  (cross)  component  normal  to 
the  orbit  plane  in  the  direction  of  the  angular 
momentum  vector  to  complete  a  right-handed 
coordinate  system,  in  kilometers. 

Velocity  vector  (up)  component,  defined  like 

U,  in  kilometers  per  second. 

Velocity  vector  (down)  component,  defined  like 

V,  in  kilometers  per  second. 

Velocity  vector  (cross)  component,  defined  like 

W,  in  kilometers  per  second. 

2 

Drag  parameter,  in  meters  per  kilogram. 


Drag  variation  parameter,  in  meters  per 
kilogram  per  day  for  the  secular  option  and 
in  meters^  per  kilogram  for  the  periodic  option. 


Note:  Any  further  headings  and  numbers  will  apply  to  nonorbital  param¬ 
eters  which  have  been  entered  into  the  solution. 

The  following  parameters  provide  information  concerning  the  size, 
shape  and  orientation  of  the  error  ellipsoid  about  the  position  of  the  vehicle 
for  a  particular  update  time. 

EIGENVECTORS  OF  U,  V,  W,  COVARIANCE  MATRIX 

These  are  the  normalized  eigenvectors  of  the  preceding  U,  V,  W 
covariance  matrix.  They  are  direction  cosine  components  with 
respect  to  the  U,  V,  W  axes  for  the  orthogonal  axes  of  the  error 
ellipsoid  defined  by  the  covariance  matrix. 
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SQUARE  ROOTS  OF  THE  EIGENVALUES 

The  square  roots  of  the  eigenvalues  represent  the  lengths  of  the 
orthogonal  semi-principal  axes  of  the  error  ellipsoid.  They 
correspond  respectively  to  the  normalized  eigenvector  columns 
immediately  above.  They  are  rotated  from  the  U,  V,  W  coordinate 
axes  as  described  next. 

TO  ALIGN  U,  V,  W  WITH  THE  PRINCIPAL  AXES 

YAW  RIGHT  (*...)  DEG 

PITCH  DOWN  (*...)  DEG 

ROLL  CLOCKWISE  (*  .  .  .  )  DEG 

These  are  a  series  of  ordered  rotations  which  will  reposition  an 
observer,  facing  initially  in  the  positive  V  direction,  to  a  new  orientation 
facing  along  the  error  ellipsoid's  nearest  principal  axis. 

(*...)  indicates  a  numerical  value  for  each  direction  of  rotation. 

POLAR  SIGMA  AND  RHO  MATRIX 


Column  Symbol 


De  sc  r  lpt  ion 


ALPHA 

1 


Right  ascension,  in  degrees 


DELTA  Declination,  in  degrees 

2 

BETA  Flight  path  angle,  in  degrees,  positive 

3  downward  from  the  local  vertical. 


AZ  Azimuth  of  the  velocity  vector,  in  degrees, 

4  positive  clockwise  from  true  north. 

R  Radius  vector  from  the  geocenter,  in 

5  kilometers. 


V  Magnitude  of  the  velocity  vector,  in  kilometers 

6  per  second. 


Note:  Any  further  headings  will  apply  to  the  drag  parameter,  the  drag 
variation,  and/or  any  nonorbital  parameters  which  were  entered  into 
the  solution. 


CARTESIAN  SIGMA  AND  RHO  MATRIX 

Column  Symbol  De  sc  r  iption 

X  Position  vector  component  from  the  geocenter 

1  in  the  direction  of  the  vernal  equinox,  in 

kilometers. 
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CARTESIAN  SIGMA  AND  RHO  MATRIX  (Continued) 


Column  Symbol 

Y 

2 


Z 

3 

XDOT  ^ 

4 

YDOT 

5 

ZDOT 

6 


Desc  ription 

Position  vector  component  from  the  geocenter 
orthogonal  to  X  and  Z  to  complete  a  right- 
handed  coordinate  system,  in  kilometers. 

Position  vector  component  from  the  geocenter 
in  the  direction  of  true  north,  in  kilometers. 

Velocity  vector  components  from  the  position 
defined  by  X,  Y,  Z  and  in  the  directions 
defined  by  X,  Y,  Z,  respectively.  In  kilometers 
per  second. 


Note:  Any  further  headings  will  apply  to  the  drag  parameter,  the  drag 
variation,  and/or  any  non-orbital  parameters  which  were  entered  into 
the  solution. 
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4.  4  PUNCHED  CARD  OUTPUT 


ESPOD  will  produce  punched  card  output  in  three  formats.  One 
format  is  used  for  the  7-card  element  set,  another  for  ESPOD  solution 
parameters,  and  a  third  for  DAC  data. 

4.4.  1  7 -Card  Element  Set  Output 

The  7-card  element  set  is  in  the  standard  SPS  element  format, 
and  will  always  be  osculating  elements.  At  the  end  of  each  differential 
correction  run  with  ESPOD,  the  resulting  best  elements  which  have 
been  tested  in  the  program  (labeled  OLD  on  the  printout)  may  be  obtained 
by  operating  DATA  SELECT  2  on  the  Universal  Buffer  Controller  (UBC). 

If  elements  are  wanted  which  have  the  last  changes  made  by  ESPOD,  a 
1  flag  in  JDC  card  Column  46  will  cause  the  latest  elements  from  the 
best  iteration  (labeled  NEW  on  the  printout)  to  be  punched  by  operating 
DATA  SELECT  2  on  the  UBC. 

The  card  field  assignments  for  the  7-card  element  set  produced  by 
ESPOD  are  given  in  Figures  4-16  through  4-22.  All  values  are  floating 
point,  and  although  the  decimal  is  indicated  to  be  in  a  given  column,  it  may 
be  anywhere  in  the  field  without  invalidating  the  format.  The  decimals  are 
shown  in  the  columns  where  they  will  appear  in  the  ESPOD  output. 

Note  that  in  all  subsequent  sample  cards  shown,  areas  which  should 
not  contain  punches  are  cross-hatched.  These  samples  may  be  used  to 
test  questionable  cards  quickly  by  placing  them  over  the  full-sized  sample 
and  noting  if  any  punches  fall  in  cross-hatched  areas. 
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Figure  4-16.  7-Card  Elements,  Card  1 
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Figure  4-17.  7-Card  Elements,  Card  2 


4-40 


7  CARD  ELEMENTS  NO.  3 


7  CARD  ELEMENTS  NO.  3 


Card  Field 

Columns 

Data 

Remarks 

3-1 

1-3 

Satellite  number  | 

Right  adjusted 

4-6 

Elements  number  J 

to  Column  6 

7 

Not  used 

3-2 

8 

Card  number  of  set 

Always  3 

9-22 

Not  used 

23-26 

Right  Ascension  of  Node  at  tQ 

Whole  degrees 

27 

Decimal 

28-36 

Right  Ascension  of  Node  at  tQ 

Decimal  degrees 

37 

Not  used 

3-3 

38-40 

Argument  of  Perigee  at  tQ 

Whole  degrees 

41 

Decimal 

42-50 

Argument  of  Perigee  at  tG 

Decimal  degrees 

51-53 

Not  used 

3-4 

54 

Perigee  distance 

Whole  earth  radii 

55 

Decimal 

56-64 

Perigee  distance 

Decimal  earth  radii 

65-79 

Not  used 

3-5 

80 

Always  E 

Figure  4-18.  7-Card  Elements,  Card  3 
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Figure  4-19.  7-Card  Elements,  Card  4 
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Figure  4-20.  7-Card  Elements,  Card  5 
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Figure  4-21.  7-Card  Elements,  Card  6 
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Figure  4-22.  7-Card  Elements,  Card  7 
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4.  4.  Z  ESPOD  Solution  Parameters 


A  second  category  of  punched  card  output  available  from  ESPOD 
serves  to  record  data  from  a  run  for  convenient  resubmission  to  the 
computer  at  a  later  time.  Since  these  parameters  correspond  exactly 
to  input  cards  described  in  Section  Z,  only  a  description  of  the  card  fields 
is  given  here.  Parameters  output  in  this  format  include  IC0ND,  DRAG, 
BISES,  BNDS,  ICTYP,  and  ITIME. 

This  punched  output  is  identified  by  a  header  card  for  each  itera¬ 
tion's  output,  so  that  the  results  of  a  given  iteration  may  be  selected  for 
r  esubmis  s  ion. 

Card  fields,  and  the  field  assignments  are  shown  in  Figure  4-Z3. 
Sample  cards  are  shown  in  Figures  4-Z5,  4-Z5,  and  4-Z6. 
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Figure  4-23.  Solution  Parameter  Punched  Card  Output 
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Figure  4-24.  Sample  Solution  Parameter  Cards,  ICOND,  DRAG 
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Figure  4-25.  Sample  Solution  Parameter  Cards,  BNDS,  BISES 
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Figure  4-26.  Sample  Solution  Parameter  Cards,  UPMAT,  SMAT 
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4.4.3  DAC  Output 


DAC  parameters  are  a  special  output  of  position  and  time  data  de¬ 
signed  to  interface  with  other  SPADATS  programs.  These  parameter s  are 
punched  on  cards ,  printed,  and  stored  in  specific  computer  core  locations  . 

DAC  parameters  are  output  automatically  when  DAC  cards  are  input 
to  supply  the  desired  update  times.  Core  storage  results  automatically, 
and  cards  are  punched  from  logical  tape  1 1  in  DATA  SELECT  2. 

4.  4.  3.  1  DAC  Punched  Cards 

ESPOD  will  punch  the  second  and  third  card  of  a  3-card  set.  The 
first  card  is  a  header  card  which  must  be  keypunched  or  otherwise  obtained 
before  using  the  3-card  set  as  input  to  another  program.  All  three 
cards  are  described  in  Figures  4-27,  4-28,  and  4-29. 
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Figure  4-27.  DAC  Card  Set,  Card  1  (Header) 
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Figure  4-28.  DAC  Card  Set,  Card  2 
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Figure  4-29.  DAC  Card  Set,  Card  3 
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4.  4.  3. 2  DAC  Core  Storage 


The  data  on  the  3-card  set,  when  the  DAC  data  option  is  elected, 
will  be  stored  starting  at  the  octal  locations  for  the  first  case  as  shown  in 
Table  4-1.  A  word  of  ZZZZZZZZ  indicates  the  end  of  the  last  set  of  DAC 
data. 

Table  4-1.  DAC  Core  Storage 


Core 

Location 


Content 


Description 


15100 


(floating  point  number) 


X  position  component  in 
earth  radii 


15101 


(floating  point  number) 


Y  position  component  in 
earth  radii 


15102 


(floating  point  number) 


Z  position  component  in 
earth  radii 


15103 

15104 


(floating  point  number)  Time  from  January  0  — 

whole  days 

(floating  point  number)  Time  from  January  0  — 

fractional  days 


Ranges  from  ZZZZZZZZ  Words  of  Z's  which 

15105  to  15124  indicates  that  this  is 

the  end  of  the  last  set. 

4 . 4  .  3  .  3  DAC  Printed  Output 

While  the  punched  cards  are  the  form  in  which  the  DAC  output  is 
used,  the  data  are  also  presented  in  the  trajectory  printout.  Figure  4-30 
shows  DAC  data  printed  output  following  a  block  of  trajectory  printout. 

At  the  bottom  of  Figure  4-30  is  the  comparison  point  from  a  subsequent 
ESPOD  run  which  had  the  data  immediately  above  as  DAC  input.  The 
black  background  for  these  subsequently  obtained  data  is  used  to  emphasize 
that  they  will  never  appear  consecutively  in  this  manner,  but  are  shown 
here  for  ease  of  illustration.  The  DAC  data  printed  are  described  as 
follows : 
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COMPARISON  POINT  1 

Column 

Symbol 

Description 

X(ER) 

Y(ER) 

Z(ER) 

The  X,  Y,  Z  position  components  in  the  geocentric  inertial 
Cartesian  coordinates  used  in  ESPOD.  They  are  true  of 

0^0  day  of  epoch.  Units  are  earth  radii. 

R(ER) 

Radius  vector  magnitude  to  the  vehicle  position,  in  earth 
radii  from  the  geocenter. 

R(KM) 

Radius  vector  magnitude  to  the  vehicle  position,  in  kilometers 
from  the  geocenter. 

Line 

Symbol 

Description 

REFR 

The  "reference"  or  value  input  to  ESPOD  or  DAC  cards. 

These  can  be  all  zeroes  when  DAC  output  is  being  obtained 
initially  for  new  update  times. 

ESPOD 

The  values  for  the  preceding  parameters  which  result  from 
this  iteration,  taken  at  the  same  update  time. 

DIFF 

The  "reference"  value  minus  the  ESPOD  value.  Indicates 
the  change  in  the  position  prediction  resulting  from  this 
iteration. 

PUNCH  CARD  OUTPUT 

Line 

Description 

1 

The  data  punched  on  Card  2  of  a  3  Card  Set  of  DAC  data, 
in  an  easily  interpreted  form. 

2 

The  card  image,  exactly  as  it  is  punched,  of  Card  2  of  the 

3  Card  Set. 

3 

The  data  punched  on  Card  3  of  a  3  Card  Set  of  DAC  data, 
in  an  easily  interpreted  form. 

4 

The  card  image,  exactly  as  it  is  punched,  of  Card  3  of  the 

3  Card  Set. 

Note:  The  parameters  punched  on  these  cards  are  described  in 
Section  4.4.  3.  1  of  this  document. 
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4.  5  MAGNETIC  TAPE  OUTPUT 


A  magnetic  tape  containing  satellite  ephemeris  data  may  be  written 
by  placing  a  "l"  flag  in  Column  53  of  the  JDC  card.  This  tape  format  is 
compatible  with  this  SPADATS  program-generated  ephemeris  tapes  of  a 
similar  nature,  and  may  be  processed  in  a  like  manner. 

This  magnetic  tape  output  format  is  described  in  detail  in  the 
"ESPOD  Mathematical  and  Subroutine  Des  cription,  M  ESD - TDR -64  - 395 . 
The  description  has  not  been  given  here  because  it  is  beyond  the  scope 
of  an  operating  manual. 
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5.  ESPOD  CONSTANTS 


This  section  lists  the  program  constants  which  can  be  changed  with 
a  99-card.  They  fall  into  the  general  classifications  listed  below. 

Astronomical  constants 
Earth-potential  constants 
Geophysical  constants 
Miscellaneous  constants 

The  utilization  of  these  constants  is  discussed  in  detail  in  ESD-TDR- 
395  ESPOD,  "Mathematical  and  Subroutine  Description,"  Section  5. 

Core 


Location 

Symbol 

Value 

Description 

00141 

CO 

E 

0.437526906  x  10-2 

Earth  rotational  rate  (rad/min) 

00142 

f 

0. 335232986  x  10"2 

Ellipticity  -  1/298.  3 

00143 

g 

32. 087598 

g  -  ft/ sec^ 

00144 

* 

K 

e 

0. 55303934  x  10'2 

Earth  gravity  constant 
(e .  r.^/min^) 

00145 

1. 0 

Earth-Earth  mass  ratio 

00146 

81 . 3015 

Earth-Moon  mass  ratio 

00147 

332951. 3 

Sun-Earth  mass  ratio 

00148 

0. 81476905 

Venus-Earth  mass  ratio 

00149 

0.  107821 

Mars-Earth  mass  ratio 

00150 

94.952 

Jupiter- Earth  mass  ratio 

00151 

317. 88665 

Saturn-Earth  mass  ratio 

It  should  be  noted  that  the  program  uses  a  value  of  Ke  given  in  SPADATS/ 
SPACETRACK  Center  Constants  (March  1964).  If  elements  derived  from 
the  obsolete  value  of  Ke  are  used,  there  may  be  a  slight  inconsistency  be¬ 
tween  the  update  as  given  by  ESPOD  and  the  update  as  given  by  SGP. 
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Core 

Location 

Flag 

Description 

00152 

0 

Earth- selector  flag  (automatic) 

00153 

1 

Moon  selector 

00154 

1 

Sun  selector 

00155 

0 

Venus  selector 

00156 

0 

Mars  selector 

00157 

0 

Jupitor  selector 

00158 

0 

Saturn  selector 

Core 

Location 

Symbol 

Value 

Description 

00159 

0. 10823  X  10"2 

See  Section  5.  8 

5*  for  interpre- 

tation  of  Earth  potential 

00160 

J3 

-0. 230  X  10 

coefficients 

00161 

J4 

-0. 180  x  10"5 

00162 

J5 

-0. 64  x  10-7 

00163 

-6 

J6 

0.  390  x  10 

00164 

-6 

J7 

-0.470  x  10 

00165 

00 

-0. 20  x  10" ' 

00166 

J9 

0. 117  x  10"6 

00167 

J10 

0.  0 

00168 

J11 

0.  0 

00169 

J12 

0.  0 

00170 

Jl,l 

0.  0 

00171 

J,  ! 

0.  0 

2,  1 

00172 

J3,  1 

0. 395  x  10"5 

00173 

J4(l 

0.  264  x  10"5 

1 

\ 

*ESD-TDR' 

-64-395, 

"Mathematical  and  Subroutine  Description." 
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Core 

Location 

Symbol 

Value 

00174 

J5,l 

0.  0 

00175 

J6.1 

0.  0 

00177 

J2,  2 

0. 232  X  10'5 

00178 

J3,2 

0.41  x  10'6 

00179 

J4,2 

0. 167  x  10'5 

00180 

J5,2 

0.  0 

00181 

J6,2 

0.  0 

00182 

X3,  1 

22.  0 

00183 

2 

31 . 0 

00184 

J3.  3 

0. 191  x  10"5 

00185 

J4,3 

0.  46  x  10"6 

00186 

J5,  3 

0.  0 

00187 

J6,  3 

0.  0 

00188 

X4.  1 

163.  5 

00189 

X4,2 

54.  0 

00190 

X4,  3 

-13.0 

00191 

J.  . 

4,4 

0.  56  x  10'b 

00192 

J5, 4 

0.  0 

00193 

J6,4 

0.  0 

00194 

X5,l 

0.0 

00195 

X5,2 

0.  0 

00196 

X5,  3 

0.  0 

00197 

X5.4 

0.  0 

00198 

J5,5 

0.  0 

Description 

See  Section  5.8  for  interpre¬ 
tation  of  Earth  potential 
coefficients 


*ESD-TDR-64-395, 


"Mathematical  and  Subroutine  Description." 
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Core 

Location 

Symbol 

Value 

00199 

J6,  5 

0.  0 

00200 

X6,  1 

0.  0 

00201 

X6,Z 

0.0 

00202 

X6,  3 

0.  0 

00203 

X6,4 

0.  0 

00204 

X6,5 

0.  0 

00205 

J6,6 

0.  0 

00206 

X2,  2 

-37.  5 

00207 

X3,  3 

51. 3 

00208 

X4,4 

50.  3 

00209 

X5,  5 

0.  0 

00210 

X6,  6 

0.  0 

00215 

re 

0. 20925738  x  10 

00216 

0. 3048  x  10~3 

00217 

6378. 165 

00219 

6378165. 

00220 

23454. 865 

00221 

57. 29577951 

00222 

6076. 1152 

00223 

0.  0 

00452 

KRMS 

1. 5 

00454 

€ 

0. 1  x  10"2 

00479 

At 

128. 

Description 

$ 

See  Section  5.  8  for  interpre¬ 
tation  of  Earth  potential 
coefficients 


Feet  per  Earth  radius 

Kilometers  per  foot 

Kilometers  per  Earth  radius 

Meters  per  Earth  radius 

Earth  radii  per  a.  u. 

Degrees  per  radian 

Feet  per  nautical  mile 

Nautical  miles  per  Earth  radius 

rejection  criterion 

Convergence  criterion 

Maximum  allowable  step  size, 
in  minutes 


❖ 


ESD-TDR-64-395, 


"Mathematical  and  Subroutine  Description." 
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Core 

Location 

Symbol 

Value 

Description 

00480 

At 

0.  0 

Minimum  allowable  step  size 
in  minutes 

00481 

y 

7  min 

0.  1 

Minimum  y  for  variable  step 
test 

00482 

103-5 

0.  1  > 

00 

1 

o 

Error  control  tolerance 

00483 

R 

o 

00 

Ratio  of  Runge-Kutta  steps 
per  Cowell 
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